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He R. d’Adhémar te un mémoire intitulé : Théorie du mouvement gyros- 
des projectiles. Les états @instabilité virtuelle. Les indélerminations. Les 
erturbations initiales. 


4 M. O. Godart présente un travail intitulé : Contribution à la théorie des effets de lati- 
tude et d'asymétrie des rayons cosmiques. V. Détermination des exposants caracté- . 
ristiques des orbites périodiques. Te 
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PE Ces mémoires paraitront dans les Annales. 


M. G. Lemaitre fait un court exposé d’un mémoire intitulé : Généralisation de la 
_ méthode de Hill, qu'il a publié ailleurs en collaboration avec M. O. Godart. 
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4 : | rie 53 | ail a 
Les observations d’équinoxes et de sol ces - 
dans le chapitre 1 du livre 3 du Commentaire — 


— 


sur l'Almageste par Théon d'Alexandrie 
Note de M. l'abbé A. Rome. 


(Seconde partie) 


Résumé. — Dans une première partie (Ann. Soc. scient, 57 (1937) 1, p. 219-236), on a 
décrit ies observations directes d’équinoxes et de solstices, et la discussion que les anciens | 
font des erreurs de ces observations. On examine maintenant les difficultés que créaient 

à Hipparque certaines mesures de la distance entre l'Épi et le point automnal. Enfin, de — 
l'ensemble des faits rencontrés dans les deux parties , on tire des conclusions sur le but 
d’Hipparque et les procédés de Ptolémée. 


11. Mesure indirecte de la longueur de l’année. — Hipparque avait 
à s'occuper d’une autre catégorie d'observations semblant indiquer une 
irrégularité dans la longueur de l’année (°) : en mesurant la longitude de 
’Epi de la Vierge, la fameuse étoile qui lui sert à exposer, si pas à calculer, 
la précession des équinoxes, on l'avait trouvée, en — 145, à 6° 1/9 à l'Ouest 
du point automnal et en — 134 à 5° 1/4 du même point. La précession ne 
suffisait pas à expliquer cela. Fallait-il introduire une nouvelle hypothèse, 
par exemple, une variation dans la position des points équinoxiaux entrai- 
nant une variation correspondante dans la longueur de l’année tropique ? 
ou bien, y avait-il erreur, et en quoi s’était-on trompé ? Hipparque admet- 
tait finalement qu’on s'était trompé. 

Théon (°”) explique comme suit la façon de mesurer la distance entre 
l’'Epi et le point automnal : 

Soit ABTA l’écliptique, la direction ABT représentant le sens direct. A 
est le point vernal, T le point automnal. Supposons une éclipse de lune se 


(©) Cfr. CLauDIt ProLEMAEI Opera quae exstant omnia. vol L pars I ed. J.-L. Heiberg. 
Leipzig 1898, p. 199,1-201,8. Nous citons en abrégé : Alm., et le numéro de la page de la 
pars 1; lorsqu'il s’agit de la pars 2 nous le mentionnons explicitement. — THÉON 
D’ALEXANDRIE, Commentaire sur le 3° livre de VAlmageste, inédit. En attendant qu’il ait 
paru nous citons les folios du Medicaeus 28.18. Le passage que nous étudions ici se trouve 
fol. 1318 sqq. 

(*) Cfr. Med. 28.18, fol. 131Y. La figure est une simple circonférence. Nous donnons la 
substance de l’explication, sans traduire mot pour mot le texte de Théon. 


Seely BS 


produisant aux environs de l’équinoxe de printemps. Au milieu de l’éclipse, 
le soleil sera par exemple en Z, un peu avant le point vernal A. La lune 
est en E, diamétralement opposée. Une étoile fixe se trouve près de la 
lune, en H (un peu plus loin, Théon suppose que cette étoile est occultée 
par la lune au moment du premier contact). 

On note l’heure du milieu de léclipse, et l’heure de Péquinoxe. Le 
temps qui s’écoule entre les deux, multiplié par le mouvement horaire 
moyen du soleil, donne la longueur de l’arc ZA = Er. 

€ Puisqu’au moment du premier contact, la lune se trouvait sur l'étoile 
» fixe H, comme le dit Hipparque dans le traité de la précession (:}) », 
d’après le temps qui s’est écoulé entre le premier contact et le milieu, on 
calcule la longueur de Parc HE, entre la lune et l’astre-repère H. 

On mesure directement la distance HO entre l’astre-repère H et l’Épi de 
ta Vierge ©. 

Une simple addition va ensuite donner @r. 

C’est cet arc OF qui a présenté des variations insolites : En (?) l’an 32 
de la 3° période de Callippe, Hipparque a trouvé |’Epi de la Vierge 6° 4/2 
en avant (c'est-à-dire à l'Ouest) du point automnal ; après 11 ans de la 
même période, il l’a trouvé toujours en avant du point automnal, mais de 
5° 1/4 seulement. 

Théon présente cette méthode comme étant celle d’Hipparque, et il cite, 
en le raccourcissant un peu, le titre du traité sur la précession, à l’appui 
du détail qu’il ajoute à l’Almageste : |’étoile-repére est occultée par la lune 
au moment du 1% contact (*). Il n’y a, en soi, auctine raison de penser que 
Théon n’a pas lu Hipparque. Lorsqu’on lit d’un trait le passage corres- 
pondant de l’Almageste, on a bien le sentiment que tout, y compris les 
deux dates d’observation, — 145 et — 134, vient du traité de la préces- 
sion. I] faut pour cela prendre |’Almageste à partir de la p. 194,16. Le 
traité de la précession y est nommé, avec son titre complet (*), puis vient 
une citation littérale du traité, p. 194,23-195, 9. La liste d’équinoxes qui suit 
(automne, p. 195,10 à 196,4 ; printemps, p. 196,5 à 196,21) est présentée 
explicitement comme provenant du même ouvrage. Ensuite vient une 
discussion par Ptolémée des résultats d’Hipparque (p. 196,21 à 197,23). 
Mais Hipparque, dit Ptolémée (p. 198,1 sqq.), n’a pas cru devoir tirer de 
ces observations des conclusions sur Virrégularité de l’année ; certaines 
éclipses de lune au contraire lui font trouver des variations dans la posi- 
tion de l’Épi. L’enchainement de ces quelques pages semble bien montrer 


(1) Cfr. Med. 28.18, fol. 1328. cita, émAoyioduevoc Tocov  GEAyvn <dmexer> Tod H, êk 
TOU AMO Tic EuTTTWOEWS Xpovou, ÊTei KATA THY GPX THs EuTTTWOEWS KATO Too GOTEPOG 
€O0ewpeito, Ka0wWe Kai 6 “Immapxog ev TH TeEpi Thc UEeTATTWOEWS TWY TPOTIKWY. 

(2) Cfr. Med. 28.18, fol. 131v. 

(*) Mais voir plus loin, p. 9; n. 1. , 

(4) mepi ueranTboews Tv Tpomkdv Kai iONMEPIVÈV ONMEtWY. 


. 


que tout provient du traité de la précession, dont le titre se retrour # 
abrégé, dans Théon. Ce dernier semble donc encore avoir lu le trai & 
d’Hipparque, et nous apporter un détail, nouveau pour nous, quil ya 
trouvé. 10. 


12. Données numériques concernant les deux observations. — 
Il est facile de réunir un certain nombre de données concrètes qui servi- 


‘à ront de point de départ à la discussion. 
| La première observation date donc de l’an 32 de la 3° période de Callippe, 
F soit —145. Nous avons déjà parlé dans la 4° partie de notre travail de 


l’équinoxe de printemps observé par Hipparque en cette année-là. Il s’est 
passé, dit Hipparque, le 27 Mechir au matin, mais le cercle d'Alexandrie 
a indiqué l’équinoxe à 5 h. (= 11 h. du matin). La date retenue par Hip- 
parque est le 27 Mechir au matin. D’après la table de Schram ('), nous 
avons trouvé que l’équinoxe avait dû avoir lieu le 24 mars —145 a 
16 h. 12 m. 8,64s. temps vrai de Greenwich, soit (*) à Rhodes, 18 h. 4 m., 
à Nicée à 18 h. 11 m., et par 30° de longitude à 18 h. 12 m. 

L’éclipse a eu lieu le 2 avril—145. D’après la table de Schoch (*) nous 
trouvons que le milieu se place à 21 h. 12,4 m. temps vrai Greenwich 
(21 h. 27,6 m. temps moyen). Demi durée de la totalité : 51,4 minutes ; 
demi-durée de l’éclipse : 111,8 minutes. Grandeur de Péclipse : 20,7 doigts. 

Au milieu de lPéclipse : 

longitude de la lune : 207°,382 ; latitude (*) : + 0°,09. | 
Au 1* contact extérieur : 

longitude de la lune : 206°,34 ; latitude : 0,00. 
rayon apparent de la lune : 0,26 (fraction décimale de degré) 
Parallaxe lunaire à Rhodes : en longitude + 0,687 

en latitude — 0,45 

à Nicée : en longitude + 0,656 

en latitude — 0,51 

Il est inutile de tenir compte de la réfraction, puisque les étoiles en sont 
affectées aussi. 


() R. Scuram, Kalendariographische und Chronologische Tafeln, Leipzig 1908, p. 358 
sqq- Nous comptons les jours à partir de minuit. 

(?) Nous adoptons pour Rhodes les coordonnées du mont Atabyrion où s’est trouvé un 
temple de Zeus, et qui est le point culminant de l’île (1215 m.), soit 27059 (=4 h. 51 m. 28s.) 
Est Greenwich, et 36°12' lat. nord. Pour Nicée, 29,8° (ou 1 h. 59 m.) Est Greenwich et 40,5° 
lat. nord (en degrés et fractions décimales de degrés). 

(*) G. Scuoca, Die Neubearbeitung der Syzygientafeln von Oppolzer dans Mitteilungen 
des astronomischen Rechen Instituts Berlin-Dahlem 2,2 (1928). 

(*) Latitudes de la lune, course de la lune pendant la demi-durée de l’éclipse, et paral- 
laxes calculées d’après les tables de P. V. NEUGEBAUER, Tafeln zur astronomischen 


Chronologie, Leipzig 1922 à 1915 et P. V. NEUGEBAUER, Astronomische Chronologie, 
Berlin 1929. 
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“4 Nous ne trouvons pas de planète (*), ni d’étoile mentionnée par l’Alma- 
… geste, qui sort occultée ou presque, au moment du 1* contact. « de la 

_ Balance (long. 207,916 ; lat. + 0,283) se trouve sur le chemin de la lune 
entre le milieu et la fin de l’éclipse. Mais la parallaxe à Rhodes et à Nicée fait 
apercevoir la lune beaucoup plus au Sud. N’empéche que cette étoile est 
_la plus proche de la lune au moment du premier contact. À moins qu’on 
_ ne suppose qu’une étoile non mentionnée par l’Almageste et visible à l'œil 
nu (il n’y ena pas beaucoup) ait été occultée. 

Finalement, l’Epi (a Vierge) en —145 a une longitude de 174,06 et une 
latitude de —1,916. 


La seconde observation est de l’an 43 de la 3° période de Callippe, soit 


(9 G. Perers et E. KNogez, Ptolemy’s catalogue of stars, Carnegie Institute Publica- 

- -tions n° 86, Washington 1915, donnent, pour chaque étoile mentionnée dans l’Almageste 
sa longitude et sa latitude, pour les années +100 et —129. Admettant pour la valeur 

» annuelle de la précession, la formule 50”,2564 +-0",0222 T, où T est le nombre de siècles 

.. compté à partir de 1900,0, on trouve que 49,8" est une bonne valeur de la précession 

> en un an pendant le 2° siècle av. J.-C. Il suffit d'introduire cette valeur dans le catalogue 
de Peters et Knobel pour —129, et l’on trouve la position des étoiles pour l'instant de 
éclipse, avec une approximation suffisante parce que la date de —145 est très proche 
de —129. 

Voici dans ces conditions, les étoiles mentionnées par l’Almageste, dont la longitude 

- était peu différente de celle de la lune au moment de l’éclipse de —145 : 

€ Serpent long. 204,349 lat. +24,149 

« Centaure long. 205,049 lat. —23,800 

8 Loup long. 205,283 lat. — 24,183 

y Balance long. 205,300 lat. + 4,600 

« Hercule long. 205,749 lat. +37,449 

u Serpent long. 206,116 Jat. 416,483 


Tt Loup long. 207,883 lat. —28,183 
« Balance long. 207,916 lat. + 0,283 
À Loup . long. 207,949 lat. —26,300 


Arg 14782 long. 298,300 lat. — 1,183 
n Hercule long. 208,600 lat. +60,549 
39 Balance long. 208,816 lat. — 8,266 
à Loup long. 208,883 lat. —21,200 
y Hercule long. 209,250 lat. +40,217 
t Balance long. 209,549 lat. — 9,749 


x Loup | long. 209,716 lat. — 29,400 
6 Balance long. 210,000 lat. + 3,600 
e Loup long. 210,350 lat. —95,000 
4 Balance long. 210,550 lat. + 6,316 
u Loup long. 210,616 lat. —28,250 


À Balance long. 210,650 lat. + 0,316 
Pour les planètes, la table de P. V. NEUGEBAUER dans Astronomische Nachrichten 248 
(1933) n° 5937 nous a donné des positions approximatives tellement distantes de celles 
de la lune, qu'il était inutile de faire des calculs plus exacts. De même pour Péclipse 
de —134. 


pr Abe | le 24 mars à 8 h. 30 m. 41 s. temps vrai. 
40 h. 2m. à Rhodes, 40 h. 9 m. à Nicée, et 10 bape i 
longitude. 


~Véclipse a eu lieu le 21 mars à 0b. 24,0 m. temps vrai de Gr iwich 

_ (0 h. 42,1 m. temps moyen). + d'a PE 
Derni-d'uréë de Ia totalité “50.7 m. ; demi-durée de té nées 116, 6m.; 

grandeur : 18,8 doigts. Late 


. Au milieu es Péclipse : 
Longitude de la lune : 176, 929 ; varias —0,13. 
A 1° contact extérieur : 
Longitude de la lune : 176,00 ; latitude —0,05. 
Rayon apparent de la lune : 0,24 degré 
Parallaxe de la lune à Rhodes : en longitude : + 0,13 
en latitude :— 0,53 
Parallaxe de la lune à Nicée : en longitude : + 0, 12 
M en latitude : —0, 760. 
Aucune planète n’est aux environs de la lune. h de la Vierge vient de 
sortir du disque de la lune au moment du 1* contact : en effet, cette étoile 
a, en —134, une longitude de 175,600 et une latitude de —0,300. Gomme 
sa grandeur est de 5,4 (6 dans l’Almageste) il faut qu’elle soit fort loin du 
disque lunaire pour qu’on la voie s’en détacher (!). 


rs 
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() Voici, obtenu de la même façon que la liste précédente, un tableau des étoiles, 


PL EM par lPAlmageste, qui avaient en —134 une longitude comprise entre 173° 
e 


& Vierge long. 172,583 lat. +8,783 
03 Vierge long. 173,016 lat. —7,666 


PA Vierge long. 173,916 lat. +3,233 
(l'Épi) a Vierge long. 174,200 lat. — 1,916 
i Vierge long. 175,133 lat ne Tes 


h Vierge long. 175,600 lat. — 0,300 
61 Vierge long. 175,683 lat. — 8,449 
m Vierge long. 177,100 lat. 1,866 
p Vierge long. 177,500 lat. +9,750 


Nous avons contrôlé nos calculs des éclipses en consultant les tables de Ginzel et de - 
von Oppolzer. D'autre part, nous les avons refaits en nous servant uniquement de l’Alma- | 
geste : L’éclipse de —145 tombe le 25 Pharmouthi 602, ère de Nabonassar, à 11 h. 44 m.” 
à partir de midi, temps vrai d'Alexandrie ; la longitude vraie de la lune au moment de la” 
syzygie vraie est 207037 et la grandeur re léclipse est de 20 doigts. L'éclipse de —134 
est arrivée le 27 Mechir de l'an 615, ère de Nabonassar, à 14 h. 24 m. à partir de midi, 
temps vrai d'Alexandrie ; la longitude vraie de Ia lune à la syzygie vraie est 176216 et la 
grandeur de l’éclipse est Pr 19 doigts. Somme toute ce résultat, comparé à ceux des tables : 
modernes, n’est pas mauvais du tout. Mais il est probable qu'Hipparque ne disposait pas 


Sy oe : : 


13. Détails de l'observation. — La position de la lune étant obtenue | 
en calculant le mouvement moyen du soleil entre l’équinoxe et le milieu 
de l’éclipse, la première chose qu’il fallait savoir était l'instant de Péqui- 
noxe. Nous en avons parlé assez dans la première partie de cet article. 
A fallait aussi noter exactement Vinstant du milieu de Péclipse. On 
___ disposait pour cela de plusieurs moyens. L’astrolabe plan, par exemple, 
; parait avoir été connu d’Hipparque. Mais Fotheringham (!) remarque qu'il 
: ne doit pas avoir été considéré comme un instrument de précision. 
a L’astrolabe (sphérique) on le météoroscope, ou toute autre variété de 
A Sphère armillaire permettait aussi de trouver l'heure la nuit pour l’obser- 
_ vation de hauteurs d’étoiles. Nous avons décrit le procédé ici-même (?), 
7 et nous n’y reviendrons plus. F 
5 que Pogo signale (*) que les Egyptiens avaient des listes de levers 
4 d’étoiles qui permettaient de trouver l’heure la nuit si la date est donnée. 
À Tous ces moyens de connaître l’heure la nuit venaient à tout le moins 
* controler et rectifier le procédé que l’on croit avoir été le plus souvent 
é employé, la clepsydre. Fotheringham et Schoch (+) se sont occupés de 
- l'erreur commise par la clepsydre. Schoch pense que cette erreur devait 
s'accumuler : si une clepsydre à marche régulière avançait de 5 minutes 
à la 1° heure, elle avancait de 10 minutes à la 2° et ainsi de suite; en 
sorte qu’au bout de la nuit elle aurait avancé de toute une heure. Aussi 
ces deux auteurs considérent-ils comme seules bien notées les heures des 
éclipses de lune qui se passent peu après le coucher du soleil. A première 
vue, cette remarque paraît évidente. Seulement, l’observation des étoiles 
fournissait un moyen facile de contrôler. Et même une erreur s’accumu- 
lant au point où Schoch le dit se contrôlerait toute seule : la clepsydre 
avançant d’une heure au bout de la nuit, on ne pouvait manquer de 
constater que le soleil n’était pas levé, et compter de combien le lever 
suivait la fin de la 12° heure marquée par la clepsydre. Une réduction 
: proportionnelle aurait alors donné une meilleure approximation de tous 
les temps observés pendant cette nuit-là. D'ailleurs Pappus () (et aussi le 


de tables aussi bonnes : en tous cas, Ptolémée dit avoir dû modifier certaines données 
d'Hipparque, relatives au mouvement de la lune. Cfr. Alm. livre 4, chap. 2 et 11; lintro- 
duction, au 5° livre, d’une seconde anomalie n’affecte pas les calculs d’éclipses. 

(‘) Cfr. Nature (Londres) 131 (1933) p. 821. 

(2) Cfr. Astrolabe et le Météorascope, dans Annales de la Société scientifique de 
Bruxelles, 47 (1927) série A, p. 88, et THÉON, Commentaire sur les livres 1 et 2 de l’Alma- 
geste, Studi e Testi 12 (1936) p. 419, n. 1. 

(3) A. Poco, Three unpublished calenders from Asyut dans Osiris 1 (1936) p. 501. 

(4) J. K. Foraertncuam, Two babylonian eclipses dans Monthly notices 95 (1955) 
p. 719. — C. Scuocu. Die säkulare Acceleration des Mondes und der Sonne dans Astro- 
nomische Abhandlungen, Ergänzungshefte zu den Astronomischen Nachrichten 8,2 
(1931) p. B2. 

(5) Pappus, Commentaire sur les livres 5 et 6 de l’Almageste. Studi e Testi 54 (1931) 


papyrus astronomique d’Oslo) indique une autre façon d'observer le 
à la clepsydre : c’est de comparer, à l’aide de l’eau écoulée, le temps _ 
à mesurer avec la longueur du jour entier. De cette façon, les défauts de _ 
graduation n’entraient plus si fort en jeu : il suffisait d’avoir un écoule- 
ment constant. Théon admet (*) comme une chose entendue que se 
tromper de deux heures dans la prédiction d’une éclipse de lune est une 
erreur qui ne peut échapper à l’attention. D'autre part Ptolémée (ie 
considère une erreur_de '/,, d'heure, dans des calculs d’éclipse de lune, 
comme tellement minime qu’ergoter là-dessus est plutôt digne d’un 
faiseur d’embarras que d’un savant cherchant la vérité. Pappus trouve 
tout naturel qu’on se trompe souvent de !/, d’heure en observant la lune. 
Les anciens pensaient donc savoir l’heure la nuit à dix minutes près. 

| Ce qui ne veut pas dire qu’ils ne se faisaient pas illusion. : 
D’après Théon, il semblerait qu’on a pris pour repère une étoile 
occultée par la lune au moment du premier contact. Pour l’éclipse de 
—134,h de la Vierge parait répondre assez bien à ce signalement. Mais pour 
Péclipse de —145 nous n’avons pas pu découvrir une étoile mentionnée 
par l’Almageste, occultée au moment du premier contact. On pourrait 
souger à expliquer autrement le texte emprunté par Théon à Hipparque : 
l'étoile aurait même longitude que la lune. La liste de la note 1, p. 9 
montre qu’en ce cas, la lune s’encadre entre u Serpent 206.116 et x Loup 
207.883; mais comme son disque ne s'étend, y compris la parallaxe 
à Rhodes, que de 206.77 à 207.29, toutes deux sont fort loin et x Balance 
(long. 207.916 lat. 0,283) est à tout prendre, encore plus satisfaisante. 
(1 faudrait donc admettre que l’étoile occultée est une étoile non cataloguée 
dans l’Almageste, ou bien remplacer l'étoile occultée par celle qui se 
trouve le plus près de la lune. 

On est tenté d'admettre cette seconde explication parce que Théon, 
analysant les causes d'erreurs dans un passage que nous étudierons plus 
loin (*), nomme en première ligne l’estimation approximative de la 
distance entre la lune et Pastre repère. Ge détail semble même exclure 
l’occultation. Il n’en est rien : décider à l’œil nu qu’une étoile entre dans 
le disque de la lune ou en sort, est réellement quelque chose d’approxi- 
matif : Schoch s’est beaucoup occupé de cette question, après Fother- 


p. 87. — S. Etrrem, Fragment of astronomy on an Oslo Papyrus dans Aegyptus 13 (1933), 
p. 479-486. — M. Pieper, Zum asironomischen Papyrus von Oslo dans Symbolae Osloenses 
13 (1934), p. 147. — O. NEUGERAUER, Zum astronomischen Papyrusfragment Pap. Osl. 73 
dans Symbolae Osloenses 17 (1937), p. 49-53. 

(‘) Toon, Clr. Med. 28,18, fol. 1338. Du reste, Ptolémée dit exactement la même 
chose au passage correspondant, Alm, 201.5. 


(?) PTOLEMBE, Cfr. Alm., 505 25. — PAppus, Commentaire sur les livres 5 et 6 de 
l’'Almageste, p. 58.11-14. 
(5) Cfr. infra p. 16, et Med. 22.18, fol. 1328 


= rh 


ingham. Mais tandis que Fotheringham admettait que l’étoile disparait 


lorsqu’elle se trouve à 3° du disque lunaire, Schoch trouve que ce chiffre 


ne vaut que pour les étoiles de 1" grandeur — il nomme l’Épi —, tandis 
que pour les autres la distance de l'étoile au bord du disque au moment 
de l'immersion augmente lorsque la grandeur diminue. Pour 19 du 
Taureau (5° grandeur) l’occultation paraît se produire lorsque l’astre est 
à 8’ de la lune. Cet arcus visionis » n’est d’ailleurs pas le même sous 
tous les climats. A propos de la question analogue du lever héliaque des 
étoiles, M. P. V. Neugebauer note que les valeurs que Schoch lui assigne 
sont bonnes pour Babylone mais pas pour d’autres régions, et que les 
expériences de Schoch montrent tout au plus qu’on peut parfois par 
hasard avoir sous nos climats un ciel aussi pur que celui de la Mésopota- 
mie (*). Pour les occultations par la lune, M. P. V. Neugebauer (*) paraît 
n’admettre comme observables à l’œil nu que celles où entrent en jeu une 
planète ou Aldébaran, Regulus, l'Épi ou Antarés, en quoi il semble 
contredit par la pratique des astronomes anciens, mieux au courant, de 
par la force des choses, des possibilités d'observation à l'œil nu. 

La distance en longitude entre l’Épi et l’astre repère se lit directement 
sur un astrolabe sphérique. Nous avons expliqué ici-même (*) le procédé 
suivi : l’on doit toujours partir d’une étoile repère dont la longitude est 
connue d'avance. C’est ce qui explique qu’on ait eu recours, dans le cas 
dont nous nous occupons aujourd’hui, à tout cet échafaudage de longitudes 
qui finissent par être déduites d’une observation d’équinoxe. 

Somme toute, en mettant ensemble tous les détails que nous pouvons 
recueillir dans nos diverses sources, il semble que le procédé d’observa- 
tion peut être reconstitué avec assez de précision. 


14. Critiques de ces observations par les auteurs modernes. — 
On pourra trouver étrange que nous ouvrions la liste des auteurs modernes 
par Delambre, l’un des créateurs, il y a 150 ans, du système métrique 
décimal. Celui qui le négligerait risquerait de faire dans Théon d’Alexan- 
drie des découvertes qui ne seraient pas très neuves : il est un des rares 
historiens qui aient bien dépouillé Théon. Le troisième livre lui a pourtant 
échappé, parce qu’il ne se trouve dans aucun manuscrit parisien. Il a posé 
toute une série de questions, au sujet des observations d éclipses dont 
nous parlons en ce moment (*). Ne possédant aucune édition de l’Alma- 
geste à laquelle on pouvait se fier, il lui arrive d’être trop sceptique, au 


(2) P. V. NEUGEBAUER, Astronomische Chronologie, p. 154. 
(2) P. V. NEUGEBAUER, Astronomische Chronologie, p. 88. a oly 
(3) Cfr. L'astrolabe et le météoroscope dans Ann. Soc. scient. Brux., 47 (1927) série A, 


. 88. A 83 0e 
} (4) DezamBre, Histoire de l’astronomie ancienne. Paris, 1817, tome 2, p. 104 sqq. ; 


voir aussi p. 207. 
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point par exemple de n'avoir pas remarqué que l’on pouvait aisément 
dater les deux éclipses, d’après ce que Ptolémée en dit. Nous ne voyons 
pas non plus qu’Hipparque connaissait mieux que Ptolémée la parallaxe 
de la lune. Mais dans la liste de questions qu’il pose, il s’en trouve encore 
que nous n’avons pas su résoudre : les deux éclipses ont-elles été observées 
dans le même lieu, ou dans des lieux dont la différence de longitude fût 
suffisamment connue? Comment avait-on déterminé l’heure du phéno- 
mène ? Il termine par une suggestion intéressante : la variation de la 
distance entre l’Epi et le point automnal ne serait-elle pas due au fait 
qu’on a négligé la parallaxe de la lune ? « 38’ de parallaxe en longitude, 
» en un sens, dans une observation, et 37° dans un autre sens dans l’autre 
» observation », suffiraient à expliquer cette variation de 1° ‘/,. 

Sans savoir mieux que Delambre en quel endroit les éclipses ont été 
observées, nous pouvons nous donner une idée de la parallaxe. On voit, 
pp. 8 et 13 ci-dessus, qu’à Rhodes comme à Nicée la parallaxe de la lune 
doit être ajoutée à la longitude, lors des deux éclipses. Il faudrait supposer, 
pour avoir une parallaxe négative lors du premier contact, que l’observa- 
tion de —134 a été faite à Babylone : la parallaxe en longitude est en effet 
de — 0,06 en cet endroit. Mais pour que la parallaxe soit à la fois négative 
et assez grande pour que les choses se soient passées comme Delambre le 
suggère, il faudrait aller beaucoup plus loin à l'Est, et sortir complétement 
des régions qui se trouvaient en relations régulières avec le monde grec. 
D'ailleurs les deux équinoxes sont très probablement observés par 
Hipparque lui-même, et les deux éclipses aussi. 

Nous n’arrivons pas non plus à expliquer par un effet de parallaxe la 
variation de 1°'/, si nous prenons, contre l’avis de Théon, le milieu de 
Péclipse : en —145, à Rhodes, au milieu de l’éclipse, la parallaxe en 
longitude est de +9,33 environ, et en —134 elle est de — 0,16 dans les 
mêmes conditions. La suggestion de Delambre ne doit pourtant pas être 
rejetée sans plus, car elle pourrait tout au moins rendre compte d’une 
partie de l'erreur. Nous y reviendrons plus loin, p. 48. 

Fotheringham (') s’est occupé ex professo des occultations d'étoiles par 
la lune mentionnées dans l’Almageste. Il écarte ces deux observations-ci, 
qui sont visiblement inexactes et mentionnées comme inexactes par 
Ptolémée. I] se borne à celles du livre 7 chapitre 3. Mais nos deux obser- 
vations d’éclipses sont examinées en passant dans le mémoire que nous 
avons eu fréquemment à citer dans la première partie du présent article (?). 
Il considère ces deux observations comme deux cas extrêmes de toute une 


(1) J. K. FOrTHERINGHAM et G. Loncporrom, The secular accelerations of the moon’s 
mean motions as determined from the occultations in the Almagest dans Monthly notices 
of the royal astron. soc , 75 (1915), p. 377-394. | 

(®) J. K. ForHeRINGHAM, The secular acceleration of the sun as determined from 
Hipparchus’ equinox observations, dans Monthly notices, 78 (1918), p.-417-419. 
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série (*) (cfr. Alm., 128.13-15). Ses données numériques viennent de von 


- Oppolzer et de Newcomb, et diffèrent peu de celles que nous avons 
obtenues. On sait d’ailleurs que Schoch a été le collaborateur de Fother- 
ingham. 


Nous pensons qu’Hipparque a rejeté l'hypothèse d’une oscillation des 


“points équinoxiaux, ce qui ressemblerait à la théorie connue au moyen âge 


sous le nom de trépidation, et dont il y a des traces dans l’antiquité (chez 
Théon d'Alexandrie). On peut difficilement dire qu’Hipparque a été amené 


à “favour a false hypothesis ». Il a discuté ces observations, comme il le 


devait pour que la question fût traitée complètement. Mais il ne s’est pas 
laissé entrainer par elles. 

D’après Ptolémée (*), qui cite le texte même d’Hipparque, la théorie de 
la précession des équinoxes était exposée par Hipparque à partir des 


variations survenues dans la distance entre l’Épi et le point automnal. 


Ces variations étaient décelées par l'observation d’éclipses de lune au 


temps de Timocharis (—294 environ) comparées à des observations du : 
temps d’Hipparque. D’autres étoiles que l’Épi ont été observées, Ptolémée 
le dit explicitement. La longitude de l’Épi avait passé, d’après l'exemple 
que donne Ptolémée, de 172° à 174. Cela ferait environ 0°0’45” par an, 
excellente valeur, en tous cas meilleure que celle donnée par l’Alma- 
geste (*) d’après le traité d’Hipparque sur la longueur de l’année, 1 degré 
au moins pour cent ans. Fotheringham croit que Hipparque n’a rien osé 
tirer de plus précis de ses longitudes de lPÉpi, vu les variations qu’il 
constatait dans ses propres mesures. Nous n’examinerons pas pour le 
moment cette question. Contentons-nous de signaler que les longueurs 
de 8° et de 6° assignées à l’arc d’écliptique qui sépare Epi du point 
automnal sont données en chiffres ronds. Ainsi, des observations d’occul- 
tations de l’Épi par la lune exécutées par Timocharis dans la nuit du 
9 au 10 Mars — 293 et du 8 au 9 Novembre — 282, Ptolémée déduit que 
l'Épi était respectivement à 82°20’ et à 82°30’ du solstice d’été en ces 
années-là, soit donc a 7°40’ et 7°30’ de l’équinoxe d’automne. 


15. Discussion de ces observations par les anciens. — Ptolémée 
fait remarquer que l’équinoxe de lan 32 de la 3° période (—145), point 


(: Cette série ne peut pas être bien longue, si elle se compose d'observations faites 
par Hipparque lui-même : on ne trouve pas un grand nombre d’éclipses de lune, visibles 
dans [a région de Rhodes, et arrivées pas très loin de l’équinoxe de printem ps, par 
exemple en février, mars et avril. Un simple coup d’ceil à von Oppolzer ou à Ginzel 
suffit à le faire voir : il n'y en a qu’une seule entre celle de —145 et celle de — 184, c'est 
celle du 1er avril —135. Encore n'a-t-elle que 8,8 doigts. Après, il faut aller jusqu’à celle 
du 29 février — 124, de 9,8 doigts. 

(2) Almageste, pars 2, pp. 12-13 (ed. Heiberg). 

(3) Almageste, pars 2, pp. 28 sqq. 


de départ de la mesure de longitude qui a placé Epi à 6°'/, du point 
automnal, est considéré par Hipparque comme bien observé (). De même, 
VPéquinoxe de l’an 43 de la 3° période (—134) qui a donné 5 pour la 
distance de Epi au point automnal, est compté par Hipparque comme 
bien observé (*). En doutant de la longueur de l’année a cause des 
variations de longitude, Hipparque admettait à la fois que les équinoxes 
étaient bien et mal observés (*). Aussi Ptolémée se hâte-t-il d’ajouter qu’il 
n’y avait chez Hipparque qu'un doute systématique. RER 

Il y avait autre chose à incriminer, dit Ptolémée (*), que les équinoxes. 
«Il est permis de trouver beaucoup plus probable, ou bien que la distance. 
» de la lune éclipsée aux étoiles les plus proches a été déterminée d’une 
» façon trop approximative; ou bien que les calculs, soit ceux de la 
» parallaxe lunaire permettant de trouver la position apparente, soit ceux 
» du mouvement solaire entre l'instant de l’équinoxe et celui du milieu de 
» Péclipse, n’ont pas été justes, ou poussés assez loin. » > | 

Théon (*) explique ce passage concrètement, avec les lettres de la figure 
- tracée au commencement ; il n’ajoute aux causes d’erreur énumérées par 
Ptolémée, qu’un seul détail, celui des erreurs possibles dans la mesure de 
la différence en longitude entre l’astre repère H et |’Epi, 9. Il rappelle 
qu'en calculant la distance à laquelle le soleil se trouve du point vernal, 
on néglige l’anomalie, et l’on ne tient compte que de la vitesse moyenne, 
qui serait elle-même calculée, si l’on doit prendre au pied de la lettre ce « 
que Théon dit un peu plus loin (fol. 432"), en admettant que le soleil 
parcourt 360° en 365 '/, jours. 

Essayons de voir quelle part chacune des causes d’erreurs détaillées par 
Ptolémée et Théon peut avoir eue dans l’erreur totale. 

En —145 Hipparque trouvait Epi à 6°,5 du point automnal, alors que 
la longitude de cette étoile était d’environ 174,06 : il exagère donc la 
distance de 0°,56 (ou 0°34’) environ. 

En —134 Hipparque trouvait l’Epi à 5°,25 du point automnal, alors que 
la longitude de cette étoile était d’environ 174°,90 : il réduit donc la 
distance de 0°,55 (ou 0°33’) environ. 

Suivons pour chacune des deux éclipses la liste des causes d’erreurs 
telle qu’elle figure dans Théon : 5 

1° Estimation trop approximative de la distance entre la lune au 1** 
contact et l’astre repère : On ne parle plus, ici, d’astre occulté. On ne 
Pexclut pas non plus. En fait, en —145 on n’a probablement pas eu d’astre 
occulté a sa disposition ; en ce cas, il fallait faire la mesure au moyen d’une 
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() Cfr. Alm., 199.5 et 196.6. 
(2) Cfr. Alm., 199.4 et 196.3. 
(3) Alm., 200.8. 

(4) Alm., 200.6-16. 

(5) Cfr. Med. 28 18, fol. 1328. 


à 
“4 
3 

F: 


sphère armillaire, et l'observation était faite à une dizaine de minutes près: 
si l’on suppose un astrolabe sphérique gradué en !/, de degré, et parfaite- 
ment construit et installé. En —134, on a probablement tiré parti de 
Voccultation de k de la Vierge. Mais alors on était victime de l'irradiation, 
qui faisait voir l’émersion seulement lorsque cette étoile de grandeur 5,4 
était à plus de 8’ du disque. 

I] s'agissait aussi d’estimer le diamètre de la lune. Hipparque le faisait 
au moyen de la dioptre (’), instrument que Pappus trouve assez peu précis, 
et dont les renseignements, pour donner le diamètre, devaient être traités 
par une méthode graphique, si l’on peut en croire Pappus. Si l’on faisait 
appel à la table des cordes, Pappus fait remarquer qu’une interpolation 
proportionnelle est aussi cause d'erreur. Mais au degré d’approximation 
auquel les anciens observaient, cette dernière remarque est purement 
théorique. Ptolémée a préféré se ser vir d’éclipses de lune pour déterminer 
le diamètre. Hultsch estime qu’Hipparque assigne à la lune un diamètre 
moyen de 3314". Cela ne nous apprend pas, même si c'était certain, quel 
diamètre Hipparque assignait à la lune au moment de l’éclipse de —134. 
1 était probablement trop fort. Nous avons vu p. 10 que le rayon de la 
lune devait être de 0°,24 au moment de l’éclipse de —134. On peut donc 
estimer que tout compris, Hipparque se trompait d’une dizaine de minutes 
dans la distance de la lune à l’étoile repère. Ici l’erreur ne pouvait jouer 
que dans un seul sens : écarter l’étoile-repère du point automnal, tandis 
qu’en —145 l'erreur de la sphère armillaire pouvait augmenter ou 
diminuer cette distance. 

2° Théon ne parle pas de l’erreur qu’on pouvait commettre en estimant 
la course de la lune entre le 1° contact et le milieu de léclipse. Cette 
course n’est cependant pas oubliée : il en parle dans le texte que nous 
avons résumé ci-dessus pp. 6-7. Mais il ne dit pas comment on la cal- 
cule : on voit seulement qu’elle est déduite de l'heure observée du 4° 
contact. Utilisait-on, comme on le faisait pour la course du soleil, la 
vitesse moyenne, ou bien procédait-on, comme l’Almageste le demande, 
à un calcul de la vitesse horaire de la lune,au moment de l’éclipse ? Par 
exemple, alors que la vitesse moyenne de la lune en longitude admise 
par Ptolémée est de U°32’56” à l'heure, au moment de l’éclipse de —134, 
on lui assigne d’après les règles de ’Almageste une vitesse de 0°30'23”, 
et au moment de l’éclipse de —145, une vitesse de 0°33'7". 

La vitesse moyenne en longitude admise par Ptolémée est la même que 
celle d’Hipparque : cela se remarque rien qu’en juxtaposant le chapitre 1 
du livre 3 et les chapitres 2 et 3 du livre 4 de l’Almageste. Si donc on a 


() Cfr. Pappus, Commentaire sur l’Almageste, p. 9. — HULTSCH, Winkelmessungen 
durch die hipparchische Dioptra, dans Abhandlungen zur Geschichte der Mathematik, 
9 (1899) pp. 193-209. — W, GUNDEL, Art. Mond, dans Realencyclopädie, t. XVI, col. 86. 
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pris la vitesse moyenne, la course de la lune en —145 a dû étre esti 
à 1°41’, ce qui est à peu près l’are que font trouver les tables 


milieu. On pouvait facilement se tromper de 10 minutes temps en 


- évaluant cet intervalle (qui est de près de deux heures dans les deux 


cas) : une erreur de 10 minutes temps entrainait une erreur de 0°5’ arc. 

3° La deuxième erreur possible signalée par Théon, git dans la mesure. 
de la distance en longitude entre l’astre repère et ’Epi. Nous venons de 
parler de cette cause d’erreur il y a un instant. Elle agit ici une seconde 
fois, qui peut évidemment tout aussi bien contrebalancer la première, que 
l’exagérer. 

4° On peut encore se tromper, dit Théon, dans la parallaxe de la lune. 
Delambre (') connaissait déjà le texte de Pappus contenant des fragments 
d’Hipparque sur la distance de la lune et du soleil. Pappus, dirait-on, 
représente Hipparque indécis sur la question des parallaxes : au 1* livre 
«des grandeurs et des distances », il suppose que la terre est un point par 
rapport à la sphère du soleil ; donc, que le soleil n’a pas de parallaxe. 
Au % livre du même ouvrage, Hipparque se risque à estimer la distance 
du soleil à 2499 rayons terrestres (*). De même pour l’éloignement de la 
lune, au 1° livre « des grandeurs et des distances » il admet un minimum 
de 71 rayons terrestres, un maximum de 83 et une moyenne de 77 ; mais 
à la fin du 1° livre € des grandeurs et des distances » il ajoute : « Dans ce 
» livre, nous sommes donc arrivés jusqu’à ce point-ci dans nos explica- 
» tions. Mais ne va pas conclure (lecteur) que nous avons déjà résolu en 
» tous points ce qui peut se dire au sujet de la distance de la lune : car 
» il reste encore en cetle matière certaine considération qui va nous faire 
» trouver une distance de la lune moindre que celle qui vient d’étre 
» calculée. » Et finalement, dans le deuxiéme livre « des grandeurs et des 
distances », Hipparque estime la distance de la lune à 62 rayons terrestres 
au minimum et 67 '/, en moyenne. Pappus ajoute : et 72 ?/, au maximum. 

Si l’on peut en juger par ces extraits, Hipparque n’est pas nécessaire- 
ment hésitant ; il pourrait fort bien avoir établi une première approxima- 
tion au {* livre, et continué son calcul au 2°. Mais nous ne devons pas 
oublier que Pappus, auteur sérieux, a eu, ou tout au moins semble avoir 


() DELAMBRE, op. cit., t. 2, p. 580. — Cfr. PAPPUS, op. cit., p. 68, et F. HuLTSCH, 
Hipparchos über die Grésse und Entfernung der Sonne, dans Berichte über die Verhand- 
lungen der k. Sächsischen Gesellsch. der Wiss. Phil.-Hist. KL., 52 (1900) p. 169-200. 

(?) 2490 est une conjecture de Hultsch. Une partie des manuscrits n’ont conservé que 
90, les meilleurs ont encore 490. — Cfr. Huttscu, Hipparchos über die Grosse, 
p. 190 et PAppus, op. cit., p. 68, n. 1 (p. 69, li. 2 des notes lire 27 au lieu de 12), 


de P. V. Neu- — 
_gebauer : 4°,04 ; et en —134, à 1°3’ ce qui est un peu trop fort, d’après | à 
les tables de P. V. Neugebauer, qui donnent : 0°,92, à condition que Pon — 
ait évalué exactement le temps qui s’est écoulé entre le 1° contact et le — 
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eu, l'original-d’Hipparque en mains. S’il le dit hésitant, il doit avoir des 

» motifs de le dire. 

Si Hipparque a hésité, a-t-il établi une table de la parallaxe, permettant 
de faire pratiquement le calcul, ou s’est-il obligé à recommencer chaque 

- fois tout le calcul trigonométrique, à partir des données qui lui semblaient 

les plus vraisemblables au moment ou il en avait besoin? Le calcul de la 

parallaxe verticale, étant donné la distance en rayons terrestres de la lune 

à la terre, et la distance zénithale de la lune, est fort simple : on n’a qu’à 

suivre le bref raisonnement du livre 5 chap. 46 de l’Almageste, qui 

comporte seulement deux résolutions de triangles rectangles rectilignes. 

Seulement, lorsqu’Hipparque a rédigé son traité des parallaxes, savait-il 

résoudre un triangle rectiligne ? Chose paradoxale, il semble avoir su 

résoudre des triangles sphériques, à l’aide sans doute du théorème dit de 

“Ménélas, lorsqu'il se trouvait encore désarmé devant un simple triangle 

À rectiligne, dont une table des cordes lui donnait la solution par une 

- réduction proportionnelle qui saute aux yeux. On sera plus disposé 

à admettre ce paradoxe si l’on considère que les Grecs n’ont jamais vu 

_ que cette réduction proportionnelle permet de résoudre un triangle quel- 

conque : ils n’ont jamais abordé que des triangles rectilignes rectangles, 
pas les autres. 

Nous connaissons par Ptolémée et par Pappus (') quelques autres détails 
sur le traité des parallaxes d’Hipparque : ainsi, il cherchait la distance 
zénithale de la lune elle-même, et pas, comme Ptolémée, celle du point 
de l’écliptique ayant même longitude que la lune. Mais il s'était, paraît-il, 
quelque peu perdu dans ses raisonnements. Et, en outre, Pappus rapporte 
des détails qui semblent provenir du traité des parallaxes : ce sont eux qui 
nous ont fait demander si Hipparque savait déjà résoudre un triangle 
rectiligne lorsqu'il écrivait son traité des parallaxes. 

2 Et comme d’ailleurs la table des cordes et le traité de trigonométrie 
d’Hipparque ont disparu, nous ne sommes, pas capables de dire quelle 
valeur Hipparque avait assignée à la parallaxe de la lune au moment du : 
4% contact des deux éclipses. 

Il est très facile de trouver quelle valeur lAlmageste assignerait à ces 
parallaxes (?), mais cela ne nous apprend pas grand chose. Il est probable, 


(2) Protimis, cfr. Alm, p. 450.11. — PAPPUS, op. cit., p. 150.20. — et Premiers essais 
de trigonométrie rectiligne chez les Grecs, dans l'Antiquité classique, 2 (1935) p. 177-192 
et en abrégé dans Annales dela Soc. scient. de Bruxelles, 52 (1932) série A, Comptes 
rendus des séances, p. 272-274. 

(2) Par exemple, à Rhodes, au moment de l’éclipse de — 145 : au 1° contact, parallaxe 
en longitude 0°43’23” qui augmente la longitude; au milieu 0°23'16", qui augmente la 
longitude ; au moment de l’éclipse de — 134, au 1°" contact, parallaxe en longitude 03127 
qui augmente la longitude ; au milieu, 0°10'40” qui diminue la longitude. Voir ci-dessus, 
pp. 8 et 10, les valeurs obtenues au moyen de tables modernes. 


i oe 


mais cela n’est somme toute qu’une impression, que la principale cause 
derreur se trouve dans la parallaxe. | 

On doit même se demander si on en a tenu compte (*) : dans la mesure 
de la distance entre la lune et le soleil, puis entre la lune et une étoile, 
servant à mettre en évidence la précession des équinoxes (*), la parallaxe 
est négligée complètement. 

5° La dernière des causes d'erreurs énumérées par Théon réside dans 
l'estimation de la distance entre la lune et le{point automnal, obtenue en 
calculant le chemin parcouru par le soleil entre le moment de l’équinoxe 
de printemps et le milieu de l’éclipse. 

Cette erreur comprend d’abord l'erreur dans l’observation de l’équinoxe. 
Nous avons vu précédemment (*) que l’équinoxe de printemps de —145 
était noté 12 heures trop tôt par Hipparque, tandis que celui de —134 
l'avait été 10 heures trop tôt. 

Ensuite, pour calculer le chemin du soleil, on partait de la vitesse 
moyenne du soleil, et pas de la vitesse qu'il possédait au moment de 
éclipse. 

Et finalement, dans Cette erreur intervient aussi l’erreur dans l’estima- 
tion de l'heure à laquelle s’est produite l’éclipse. 

Le mouvement du soleil étant lent, toutes ces causes d’erreurs agissaient 
peu, et cela pouvait passer pour un bon subterfuge de calcul, d'utiliser le 
mouvement du soleil pour trouver la position de la lune. 

Pourtant, le total de cette erreur semble avoir été assez considérable : 
supposant l’heure de l’éclipse bien notée, et comptant 0°59’ pour le dépla- 
cement moyen journalier du soleil, en —145, on trouve pour la course du 
soleil entre le 24 mars vers 6 heures du matin, date de l’équinoxe d’après 
Hipparque, et le 21 avril vers 23 heures, date de l’éclipse à Alexandrie ou 
à Rhodes d’après nos tables, un arc de 27°15'17’, ce qui met la lune sur 
207°15’ au lieu de 207 382 (= 207°23’) que nous avait fait trouver la table 
de Schoch. En utilisant simplement les tables d’éclipses de l’Almageste, 
on trouverait que la lune est sur 207°37' au moment de la syzygie vraie. 
Dans les mémes conditions, en —434, on trouve pour la course du soleil 
jusque minuit de la nuit du 23-24 mars, où se place l’équinoxe d’après 
Hipparque, depuis le 21 mars vers 2 heures du matin, date de l’éclipse en 
Proche-Orient d’après nos tables, un are de 1°59’, ce qui met la lune au 
moment de l’éclipse sur 178°8’ au lieu de 176.992 (= 176°55’) que font 
trouver les tables modernes. En se fiant uniquement aux tables de 
l’Almageste, on trouverait 176°16. 


(') Toutefois Ptolémée et Théon parlent*de l'erreur commise dans la parallaxe et pas 
de l’erreur commise en négligeant la parallaxe. 

(*) Alm. livre 7 chap. 2. Voir ci-après, p. 21. 

(°) Voir ci-dessus, {re partie du présent article, dans Annales de la Soc. scient., 57 
(1937) série I, p. 215, n. 3. 
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€ pote est aadente v af fe tables, et non ue d'après une | | 

_ observation d’équinoxe, combinée avec la course moyenne ; la longil 

de la lune (au 4* quartier) est évaluée en mesurant directement à l’astr 

_ labe la distance lune-soleil ; et finalement, l’on mesure la distance entre _ 
la lune et une étoile. Le ne comme nous l’avons signalé il y a un 
instant, sans tenir compte de la parallaxe ni de la réfraction. 

__- En mettant ensemble tous ces éléments, nous n’arrivons pas à retrouver 

_ tout le détail des erreurs : certains De restent problématiques. 

Il n’est pourtant pas superflu de grouper tous nos résultats : 


Pour Péclipse de —445 : 
__ 4° Distance de l’équinoxe à la lune au milieu de l’éclipse : 


2° Course de la lune entre le 4° contact et le milieu : 
semble avoir été calculée assez exactement ? 0°0/ 
3° Parallaxe de la lune: 
nous ne savons pas quelle valeur lui a assigné Hipparque ? 
4° Distance de la lune à Pétoile repère non occultée : 
Perreur ordinaire des arcs mesurés à l’astrolabe 
soit au moins + 0°10 
5° Distance de l’étoile repère à Epi: 
même erreur qu'au précédent soit au moins + 0°40’ 
En tout, on a placé l’Épi sur 173°30/ au lieu de 1744", soit 0°34’ trop peu. 
On voit que l’erreur s’est produite en faisant les HiReses observations 
et les calculs d’une façon pas trop mauvaise. 


(1) Alm., , éd. Heiberg, vol. 2 p. 14.1. — Voir aussi MANITIUS, Ptolemiius Handbuch der 
Astronomie, Leipzig 1913, vol. 2, p. 13. 

FOTHERINGHAM, The secular acceleration of the sun as determined from Hipparchus’ equi- 
nox observations. Monthly notices of the R. astr. soc., 78 (1918) p 419, trouve que Ja 
facon habituelle de raconter la découverte de la précession ne doit pas étre authentique : 
Hipparque aurait, dit-on, admis Ja précession pour expliquer que la longitude de PEpi 
- avait changé depuis le temps de Timocharis. Mais Hipparque ne pouvait pas accorder sa 
-_ confiance à une méthode conduisant, d’après ses propres déclarations (qu’il fait dans le 
traité de la précession) à des résultats si discordants. 

Il nous semble qu’on doit distinguer l'exposé qu’Hipparque fait de sa théorie, des 
recherches qui l’y ont conduit : il est très probable qu'il ne s’est pas contenté de mesurer 
une seule longitude d'étoile. En tous cas Ptolémée, à l'endroit que nous venons de citer, 

parle explicitement de mesurer la longitude de plusieurs étoiles. D'autre part, la méthode 
de mesure au chapitre de la précession, n’est pas identique à celle qui est suivie ici : on 
- élimine la lune éclipsée, l'étoile occultée, et la position du soleil comptée à partir d’une 
seule observation d’équinoxe. 
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on met la lune sur 207°15 au lieu de 207°23’ soit 0°8’ trop peu 


semble avoir été pi prise un ‘peu Pt ongue, xempl 
)°5' de trop, ce qui diminue de 0°5' la longitude « delay res Sant 
une, | soit CL 0) 
3° Parallaxe de la lune : ' aie i 
4° Distance du centre de la lune à l’étoile repère : a ds! 
l’étoile repère vient d’être occultée ; en y comprenant ath oe 
une légére erreur sur le diamétre de la lune, on la 
compte comme sortant du disque lorsqu’elle est déja 
à 10° au moins du bord, : soit donc (0°10 de tro 
5° Distance de l'étoile repère à Epi: sé. 
l'erreur ordinaire des observations à l’astrolabe, au 
moins + 0°10". 
En tout, on a mis l’Epi sur 474945! au lieu de 174919’, soit 0°33’ de trop. 
Les erreurs que nous pouvons chilfrer, soit 1°43" + 9°40" — 0°5’, ont 
donc été compensées en partie par des erreurs d'observation et de calcul 
de parallaxe. L’observation de —134 est done beaucoup plus mauvaise, 
dans ses détails, que celle de —145, quoiqu’au bout du compte, elle arrive 
à n’être pas plus distante, par excès, du résultat exact, que la première 
ne l’était par défaut. 


16. But des observations d’Hipparque. — Enumérons maintenant 
quelques conclusions qui semblent se dégager de ce qui précède. 

Hipparque admet que la durée de Pannée tropique est constante. Cela 
ressort de l'ensemble du chapitre que nous venons d’étudier dans Ptolé- 
mée et dans Théon. Ptolémée fait remarquer, entre autres choses (*), que 
s’il avait admis une irrégularité, il aurait, conformément à l'usage, 
introduit une seconde anomalie du soleil, destinée à mettre la longitude 
calculée d'accord avec la longitude observée. Par parenthèse, puisque 
cette idée se dégage nettement de l’ensemble du chapitre, on ne peut pas 
dire (*) que Ptolémée prête à Hipparque des opinions qu’il n’a jamais 
soutenues, pour se donner le plaisir de les réfuter. Ptolémée pense 
qu'Hipparque a soulevé ces questions par honnêteté scientifique, pour 
mentionner aussi les faits qui semblaient donner en démenti à ses théories. 

Ce qu’un ancien aurait appelé une irrégularité dans la durée de l’année, | 
ce serait, par exemple, une variation de deux heures au moins : c’est ce 


(1) Alm., 200.17. 


(?) Cfr. TANN ER Recherches sur l’histoire de l’astronomie ancienne, Paris 189%, 
pp. 148 et 150. 


3 chiffre que prend Ptolémée lorsqu'il veut donner un échantillon @irrégu- 
_larité qui sefait observable ('). Savoir sil y avait des irrégularités plus 


petites était complètement hors de leur horizon, puisqu'ils se rendaient 
bien compte qu’ils ne savaient pas de combien l’année était plus courte 
que 309 !/, jours. . 

Hipparque donc ne se demandait pas quelle mécanique il allait adopter 
pour rendre compte de ces variations ; ou tout au moins, s’il a pu (disons 
même s’il a dû) par moment, lorsqu'il étudiait la question, songer à l’un 
ou l’autre arrangement, il ne s’y est pas arrêté, et c’est dépasser les textes 
pour tomber dans la conjecture, que de nous montrer un Hipparque 
pensant admettre une oscillation des points équinoxiaux (?). 

Tannery (*) a fait une hypothèse ingénieuse sur ces considérations 
d’Hipparque : On a pu constater, dans les pages qui précédent, qu’une 


bonne partie de ces détails nous viennent du traité de la précession. 


Hipparque, suggère Tannery, aurait:voulu montrer que même en mettant 
les choses au pire, il restait encore dans la longitude des étoiles une 
variation qui ne pouvait s'expliquer qu’en admettant la précession. L’on 
comprendrait en ce cas qu’Hipparque ait pris des résultats d'observations 
complétement aberrants et qu’il rejetait comme invraisemblables. Cette 
hypothèse ressemble à celle de Ptolémée que nous venons d’indiquer. Pour 
Tannery, le but d’Hipparque (*) aurait été de déterminer la limite infé- 
rieure de la fraction à retrancher à 365 ‘/, jours pour trouver la longueur 
de l’année tropique. Nous ne voyons pas cela, ni dans Ptolémée, ni dans 
Théon, nos deux sources. |] nous semble qu'Hipparque a cherché à trouver 
le plus exactement possible la longueur de l’année tropique. Il n’y est pas 
arrivé de facon satisfaisante, ni Ptolémée non plus. Ce dernier le constate 
avec beaucoup de franchise, et sans chercher à jeter de la poudre aux 
yeux. C’est même celte façon d'exposer honnêtement les résultats auxquéls 
on est parvenu, et d'espérer qu'un successeur fera mieux, qui fait le 
charme de cet intéressant chapitre. 

Nos bases ne sont pas assez larges, dit-il en substance, pour arriver à des 


résultats satisfaisants. Mais nous essaierons une synthèse, et nous ne nous: 


désintéresserons pas de la question. « Quant aux affirmations valables 
» pour toute l'éternité, ou bien même valables pour un temps considéra- 
» blement plus grand que celui sur lequel portent les observations, 


1) Alm., 201.1-8. ; 

(2) FOTHERINGHAM, The secular acceleration of the sun, dans Monthly notices of 
R. A. S., 78 (1918) p. 407. L’oscillation en question ne serait pas Ja trépidation, mouve- 
ment a longue période, rejeté par Hipparque et Ptolémée, qui lui ont préféré la théorie 
de la précession. Fotheringham a repris ce nom pour désigner des fluctuations dans le 
mouvement du soleil. 

(3) TANNERY, Astronomie ancienne, p. 150. 

(4) Cfr. Alm., p. 207. — TANNERY, Op. cif 19: 
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. mais pas pour extrapoler. Si on la fait commencer à l’an 1 de Nabonassar, la table des 
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€ M. Marcot pense, et nous avions eu cette idée avant de lire ant 
» ques, dit Delambre, que ces observations sont supposées ou du moi 
» qu’elles ont été calculées tout exprès de manière à donner la urée € 
» l’année telle qu’on la connaissait. Il était bien maladroit d’altérer 
» intervalle si facile à vérifier ; il est plus naturel que Ptolémée ait cal 
» le temps de ses équinoxes par l'intervalle qu'il supposait; alors il devai 

» toujours retrouver la même année; et sil n'avait pas commis der ou 
» sur Vintervalle, la supercherie n’eût jamais été découverte. Il s’est 


» trompé. La fraude est reconnue... » (*) e ; CRAUE e 
Manitius (*) soupçonne Ptolémée de ne pas avoir observé les solstices, h 
sans quoi il aurait vu qu’il devait admettre un mouvement de la ligne des 
apsides de son excentrique. C’est, nous semble-t-il, exiger trop de Ptolémée. 
On dit souvent que le catalogue d'étoiles de ’Almageste est copié dans 


Hipparque, et que Ptolémée s’est contenté d’ajouter aux longitudes une | 


(1) Alm., 203.5-8. Ceci semble vouloir dire que les tables sont bonnes pour interpoler 


mouvements du soleil s'arrête à 63 ans après J.-C., longtemps avant Ptolémée. Serait-elle 
reprise à Hipparque ? Mais en réalité cette table n’a pas de commencement défini : on 
doit ajouter, à la fin de l’addition, la position du soleil au 1 Thot de l’an 1 de Nabonassar. 
On pourrait ajouter la position du soleil à n’importe quelle époque. Toutes les tables 
sont écrites en colonnes de 45 lignes : c’est cela qui détermine sans doute le nombre de 
périodes de {8 ans qu’on peut y inscrire. La table des syzygies, qui procède par périodes 
de 25 ans, comme tes « Tables faciles », s’arrête à 355 après J.-C. à 

(?) DELAMBRE, Histoire de l’astronomie ancienne, t. 2, p. 110. — Cfr. ci-dessus, n° 9 
(dans la 1"° partie du présent article, p. 226). 

(*) MANitius. Ptolemäus Handbuch der Astronomie. Leipzig 1912, vol. 1, p. 428, n. 23. 
Nous aurons probablement à revenir sur cette question.Mais rien qu’avec les éléments 
réuuis dans le présent article, on entrevoit deja la justification de Ptolémée. Manitius 
pense que Ptolémée ne pouvait manquer de voir que depuis Hipparque le printemps 
s'était raccourci et l'été allongé, ce quil’aurait amené à admettre un mouvement de l’ex- 
centrique du soleil. Admettons les données de Manitius : Ptolémée place l’apogée du soleil 
sur Gémeaux 5°30’, alors qu'il aurait dû trouver Gémeaux 10°22’. Il suffit de faire en sens 
inverse, et en utilisant les tables de l’Almageste, les calculs d’Alm., p. 235 (= Manitius, 
p- 168 sqq.) pour trouver que cette position de l'apogée correspond à un printemps durant 
environ 94 j. 8 h. et un été de 93 j, alors que Ptolémée admet que depuis Hipparque les 
saisons n'ont pas changé, et que le printemps a toujours 94 1/2 jours et l’été 92 1/2 jours. 
Le lecteur remarquera immédiatement qu'Hipparque ayant donné les dates d’équinoxe 
et de solstices à 6 h. près (cfr. {re partie de cet article p. 217), Ptolémée donne ici, fort 
logiquement, les longueurs idles saisons à 1/2 jour près. Mais comment veut on qu'il ait vu 
que le printemps s'était raccourci de 4 heures depuis Hipparque, s’il en connaissait la 
longueur, du temps d'Hipparque, à 12 heures près seulement ? Ceci en supposant que 


Ptolémée était capable d'observer les longueurs des saisons avec une approximation 
suffisante. 


| 
{ 
| 


iors. 


- constante pour tenir compte de la précession. Fotheringham (!) signale 
que cette opinion est déjà réfutée par Laplace depuis 1824; on retrouve, 
sous-jacente aux données du catalogue d’étoiles, une erreur dfobtervition 
- d’équinoxes correspondant à celle que Ptolémée a commise en +132. 
Évidemment ceci ne veut pas dire que Ptolémée a fait son catalogue 
d'étoiles sans tenir compte des travaux antérieurs. Mais cela semble 
indiquer qu’il ne faut pas trop se hater de déclarer que Ptolémée ment 
lorsqu'il dit avoir fait des observations. 

Il semble que la clef du procédé est fournie par une petite remarque 
faite en passant (*), à la-fin du chapitre 1% du 5° livre, tellement acciden- 
telle que Pappus n’a pas cru devoir en dire un mot dans son commentaire : 
» « Des nombreuses parallaxes que nous avons observées, (la lune étant) 
_ » dans des positions du genre (de celles qui ont été expliquées dans le 
» » chapitre 12), de nouveau, nous allons en détailler une, au moyen de 
1 » laquelle nous donnerons à la fois un échantillon de la facon de calculer, 
” >» et nous établirons la preuve de ce qui suit dans l’ordre naturel des idées.» 

Pour l’exposé, on choisit donc une observation; mais on en a fait 
plusieurs. Le résultat dépend donc de l’ensemble des observations. La 
question est de savoir comment il en dépend, et quelle observation est 
choisie pour exposer la théorie. 

Dans le cas présent, le résultat de Ptolémée, déduit de l'ensemble de sa 
> documentation, est somme toute négatif : le degré d’approximation auquel 
~ ila pu arriver ne lui permet pas de dire de combien l’année tropique est 

plus courte que 365 ‘/, jours (*). Il constate seulement qu'elle est plus 
- courte. Dans ces conditions, il pense sans doute qu'il n’y a aucun avantage 
» à changer les données d’Hipparque : l’approximation qu’il proposerait 
- serait tout aussi mauvaise que l’autre, et après tout, l’on pourrait aussi 

bien écrire n'importe quelle autre fraction Bape entre les limites des 


résultats d’observation (4). Hipparque prend 30 de jour (°). Ptolémée 


adopte done ces chiffres-là. 

(1) ForuerRINGHAM,The Secular acceleration of the sun, p. 422. Cfr. un reproche ana- 
logue émanant de Newcomb réfuté par FOTHERINGHAM et LONGBOTTOM, The Secular 
acceleration of the moon’s mean motion dans Monthly notices 75 (1915), p. 378. 

(2) Alm., 408, 5-9. 

(3) Alm., 201,32-202,9 4 

(+) Il est curieux de mettre à côté de ces constatations la phrase suivante de 
P. GoEDSEELS, Théorie des erreurs d'observation et éléments du calcul des probabilités. 
Paris 1914, p. 148 : « Le livre précité de M. Colardeau contient de nombreux passages 
» confirmant nos idées et celles de M. Estienne. Nous lisons notamment, à la p. 176, à 
» propos des cinq résultats suivants : 0,05739 ; 0,05691 ; 0,05685 ; 0,05679 ; 0,05712, que 
» les deux derniers chiffres décimaux, ne correspondent à rien de réel et que des 
» chiffres quelconques, mis à leur place seraient tout aussi bons. » 

(5) Atm., p. 207. 
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. West choisir arbitrairement, si l’on veut. Ma 
_inventer des observations qu’on n’a pas faites, et cela ne peut pas 

une supercherie : c’est application logique d’un procédé dont |. 
se trouve clairement énoncé ailleurs, et qui n’avait probable 
besoin d’être énoncé, parce qu’il était entendu que pour expose I 
ment ses résultats, il fallait élaguer le détail des nombreuses obsenvalag ns 
et de tous les tatonnements dont on était parti. "re 
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et d'asymétrie des rayons cosmiques 


” V. — Détermination des exposants caractéristiques 
des trajectoires périodiques 


PAR 


ODON GODART 


INTRODUCTION 


Parmi les trajectoires de particules chargées déviées par le champ 
magnétique terrestre, les trajectoires asymptotiques aux orbites pério- 
diques forment, au point de rencontre avec la terre, le cône appelé de 
première espèce dans les contributions précédentes. Ce cône sépare les 


orbites venant de l'infini de celles venant de la terre. La connaissance des 
orbites asymptotiques détermine donc la répartition du rayonnement 
| cosmique. 


Dans la contribution IV (), MM. les professeurs Lemaitre et Vallarta 
ont mis au point une méthode d'intégration numérique des trajectoires 
asymptotiques, lorsqu'on les connaît déjà dans les environs immédiats de 
lorbite périodique. Cette détermination préliminaire est le sujet de la 


présente contribution. 
Nous avons dans une première partie calculé neuf orbites périodiques, 


puis déterminé les exposants caractéristiques et les solutions des équations . 


aux variations de six d’entre elles. Cette première partie fut présentée 
comme dissertation de licence en sciences mathématiques en 1935. La 
méthode employée fut reprise pour le calcul d’une famille de trajectoires 
asymptotiques par MM. Lemaitre et Vallarta (’). 

La deuxiéme partie comprend la réduction des équations aux variations 
à une seule équation à coefficient périodique, équation que Hill a étudiée 
dans son mémoire sur le mouvement du périgée de la lune (*). Nous 


() Ann. Soc. Sc., A, 56, 1936, p. 102. 
(2) Acta mathematicu, 8, 1886, p.1 et WuitTTAKER et WATSON, Course of Modern 


Analysis, p. 414. 
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Les équations auxquelles obéissent les s particules cosmiques, s 
un plan méridien accompagnant la particule 


; Hé. ab han at 
(1) dai de wee aaa Ok 
dx\’ (dd nits NT 
: se) ay <5 (à) =P = 167 a+ dé ben cos} ” 


où y, est le paramètre de Stormer donnant la valeur du moment de la 
quantité de mouvement à linfini ; o une variable indépendante liée au 
temps ; À la latitude et où 1 


PETER 


r étant la distance de la particule au centre de la terre. ; 
Connaissant x et \, il nous est possible de déterminer la longitude @ par 


@) SE grt te 
do cos*h 

M. Lemaitre a exposé dans la contribution III (‘) et repris dans la 
contribution IV, une méthode permettant de calculer les solutions pério- 
diques en se donnant la valeur d’un paramètre n. Lorsque n augmente, 
Ÿ, diminue jusqu’à un minimum pour augmenter de nouveau. 

Les orbites périodiques pour les valeurs de n inférieures à celle corres- 
pondante à la valeur limite de 7, sont instables. Les autres orbites sont 
évidemment stables et n’admettent done pas de trajectoires asymptotiques. 

Nous avons donné à n les valeurs : 


0,01; °0,02; 0,08: 0,04; 006:/006-2707 


et étudié ainsi des orbites périodiques choisies dans toute la zone d’insta- 
bilité. 


(1) Ann. Soc. Sc., A, 54, 1934, p. 194. 


. po APE Aa Pe : 
B= By + 4,0, + 2,0, 
ARS + mS, + oS; + + 


CG = cosk (wo + 9), S, = sink (wo + 9), 


 w étant une constante différente pour chaque orbite périodique, et @ une 
constante arbitraire. 
__ La méthode de M. Lemaitre s’est révélée très précise pour le calcul de À 
_ déterminé toujours avec 5 décimales exactes. Le calcul de x a nécessité 
pour les plus grandes valeurs de n des corrections effectuées comme il est 
- indiqué dans la contribution IV, corrections atteignant la cinquième 
~ décimale. 
Dans le tableau des résultats figurent, en outre des constantes déjà 
_ définies, r, valeur de r pour 6 —0, z, et wy, constantes telles que 


sin À = Vz, sin [wo + y, sin 2w0]. 


TABLEAU I 
n+.0,04 0,02 0,03 0,04 0,05 0,06 0,07 
| | 0,964.053| 0,929.898 | 0,895.050 | 0,860,846 | 0,828.245 0,799.98 | 0,792.071 
= ;,| 0,973.854 | 0,945.009 | 0,912.620 | 0.875.311 | 0,820.479 | 0,772.002| 0.670.358 
= | 0,879.678 | 0,892.772 | 0,904 745 | 0,914 183) 0.918.451 | 0,910.632 | 0,857.387 
> | — 2993} — 5.000) — 8.470} — 12880) — 18.948 — 28151) — 47.814 
| 46.528} 92712| 0,128.364| 0.183.159) 0,226.391 | 0,266.149 | 0,292.240 
© | 0,632.910| 0.567.971 | 0.497.734 | 0,420.653 | 0,834.620| 0,234.97 | 0,102.44 
D «| —1933| — 4.063] — 6808) — 10.501) — 15.672) — 93.74 | — 42.04 
| 0.216.729 | 0.307.351 | 0,377.196| 0.435.595 | 0.485.726 | 0,527.13 | 0,548.56 
= |” — 677; —2018| —3.933| —6.a78| — 9.814) —14.31 | — 24.45 
E 3 24 72 164 319 57 1.02 
ep, _9 eH, th te —7B3 


La détermination de l’orbite périodique limite exposée dans la troisième 


partie nous a amené à calculer deux orbites : 


Pune est près de l'orbite 


limite; l’autre où n= 0,067.551 est la véritable orbite limite. Nous 
publions ces deux résultats dans un second tableau entre deux approxima- 
tions de l'orbite limite calculées par M. Lemaitre. 
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| 0,989 248 | 0289278 | 
_0,142.284 | 0142103 | O,141 
— 34.970 | — 34.987 | = 3507 | ps ME 

PB Eee CAEN EE rag “br 
0,547.790 | 0,547.805 | 0,547.847 | 0.547.991 Al : 

— 19.118 | — 19455 | —419474 | — 419979 
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§ 2. — Forme de la solution des équations aux variations 


Remarquons que les équations (1) peuvent se mettre sous la forme 
canonique en considérant l’hamiltonien : 


H=P— 51) )+(3 yi} 


D’après la théorie des exposants caractéristiques ('), les solutions 


périodiques des équations (1) admettent des exposants caractéristiques — 


égaux et de signe contraire ; de plus, deux d’entre eux sont nuls, car il 
existe l’intégrale H = 0. 
Si nous posons 
e=EF+e, cals eal i) 
où & et n sont les solutions périodiques, =, et n, des quantités petites telles 
que nous pouvons négliger le carré. Les équations aux variations s’écrivent: 


dE [OP ap 
do — soy, ee en 
ayy dP 

dP hie an) + + (oar) 1 


ae. ae dyn dn oP oP 
5 no see dll! Al ee Rd =a 
(9) do do ' do do \am jot où 
SE D038 10 2080 MAINS TON SNS 
() PoiNcaRÉ, Méthodes nouvelles de mécanique céleste, tome I, p. 162. 
(2) n sera dorénavant employé dans ce sens. 


(4) 


ee 


a 
"et nous en ron 


que. Nous alveriept active ainsi 
“2 &) = My + al + UG, = 


(eax) 


Ci + €,53-+ 6553 + + 


»P 
Car)= Bat ds + b,0, + + 


3 Outre la solution banale | 
LE dE dn = 


Ass Vly gh es 


les HA TAN (3) admettent des solutions de la forme 
3 (7) E = e2° E[ySx + 0x] , & pair ; 
a n, = e29 Z[uSx + vex) , k impair. 

§ 3. — Calcul de Vexposant caractéristique 


Bini gauiabsis les expressions (6) et (7) dans les équations du mouvement; 
effectuons les dérivations, les produits de sinus et de cosinus. Comme cette 
relation doit être vérifiée pour toute valeur de 6, les coefficients des S, et 


, doivent être égaux dans les deux membres. Nous avons multiplié les 
- équations"par ? pour simplifier les écritures. 


D, 2(A—0%?) z, 0, u,+D, Zs sos 
À (A, —22) v,—B,Q u, CA 
e 2C, 2,+B,2 v+(A,—222) u+C, 2 = 
e 2D, Z, +C LE (A—20) 2 00000 
a +C, v, + 9B, art (eB)  » =i 
E Non!’ | ae ne eee en 
4 

ou 

DA =a, Clee cy =o, B, = 4w 
A, = 26, + b, + Qw? C= '-2Gikle RS 0, 

DA,— 9, —b, +2w? - (= ¢+¢ D,=6,—6, 

» À, — 2a, + a, + 8w’ 

» A, = 2a, — a, + BW” 

- La condition de compatibilité de ces équations linéaires homogènes se 


met sous la forme d’un déterminant infini égalé à 0. 


st)" tee il de 
4 


» 


@) Ay—- = (CA, + DID, — 20, DA) 


’ 


où les mineurs A sont des fonctions de 2. | " 
= Une première approximation se détermine par itération au moyen de la 
formule ci-dessus. On prend comme première valeur de S°, A, et on cal 
cule les valeurs des A; l'équation (9) donne une deuxième valeur au moyen 
de laquelle on calcule les A, et ainsi de suite. | 

Ce procédé converge rapidement. 

Les mineurs de ce déterminant par rapport aux éléments de la première 
ligne nous donnent les valeurs approchées de 2, My, v,, Ye, 22. Portons ces 
valeurs dans les deux équations suivantes arrêtées aux termes en yu, et V5. 
Nous pourrons résoudre par rapport à v, et u,. En introduisant ces valeurs 
dans la première équation, nous obtiendrons une nouvelle approximation“ 
de Q?. Nous pourrons calculer de nouvelles valeurs de u,, V,, Ys, 223 Ma, Va. 
Lorsque ces valeurs auront convergé, le même procédé se continuera avec 
les deux équations suivantes, et ainsi de suite. | 

Ce calcul a été effectué pour les 6 premières orbites du tableau |. Les“ 
résultats sont publiés dans les tableaux III et 1V à la fin de la contribution. 
L’exposant caractéristique obtenu par cette méthode est écrit Q,. Nous 
avons toujours posé 2, = 0,5. | 

Cette méthode devient extrêmement pénible quand on s'approche de 
l’orbite limite ; les résultats deviennent incertains, car il est nécessaire de 
calculer des termes de plus en plus loin dans le développement de 2, n,. 
Aussi nous avons été amené à chercher une vérification de nos calculs. 
Cette vérification fut trouvée dans l’application des méthodes de Hill à 
notre problème. 


DEUXIÈME PARTIE | 


APPLICATION DE LA MÉTHODE DE HILL 


$ 4. — Déplacement normal à l'orbite périodique 


La théorie de Hill donne pour un degré de liberté une détermination 
de l’exposant caractéristique à la fois plus rigoureuse et plus simple que 
la méthode d’approximation précédente. La généralisation de cette 
méthode à un système de n équations linéaires homogènes du second 


re 


(0) 


q 
1 


; f mo) 25) Sy m.n 


pré tée à Pacademie Royale 
+ Nous ¥ pouvions donc appliquer le 


ill aux nations ae 
Mais il a paru préférable, pour simplifier les calculs numériques, de 


à ramener les équations aux variations à l’équation originelle de Hill. 


Posons — 
+n 
+ Has P ip {oa 
do ries] 
. Les équations aux variations deviennent 
Ey PTE Pan 


ni tt us pen, 


avec la condition d’énergie 


> (41) 22) + n'n1 = Profi + Pom 


L'écart £,, n, à l’orbite périodique peut se décomposer en déplacement 
tangentiel « et en déplacement normal v, par les formules 


(19) EAs 
n,= nu + Ev 
La condition d’énergie donne 
(13) (+ n°) w = 2 (EP, —n'Pi) v 


Multiplions les équations aux variations (4 0) respectivement par n’ et £’ 
ajoutons-les, et remplaçons £,, n,, wu’ par leurs expressions tirées de (12) 
_et (13) nous obtenons 

En!—n'?" 


w ES) v'+9 2 (nn"+£ ae v'+[4P,,£n AR ++ n*)(Poo—Po2)]v=0 


Pour faire disparaitre la dérivée première posons 
1 
2 
yas we (B® ne) 


Ce qui revient a faire la transformation 
| Pos Eu — Bw 
( n=nu+aw 


où a et B sont les cosinus directeurs de l’orbite périodique. 


(1) Bull. Ac. roy. Belg. (Classe des Sciences) 1, 1938, p. 19. 
LVI, I 


7 


a% 


© Nous représentons par R le rayon de courbure de | 


valeur finie 


7 1» = ut a 
ak asl ten pp ier 
punks 4 \H a? yy 


_Remarquant que | & 
ë Re ee 
Péquation du déplacement normal devient | 0 Oe 
(15) w" + [3(aB' — Ba’)? + 2aB P,, — a’P,, — B’P,,] w= 0 TE 
Introduisons langle x que fait la tangente à l’orbite périodique avec 
axe des À, l'équation (15) peut s’écrire : .. (Ur 


w"+ [3x + sin 2xP,, — sin’xP,, — cos’xP,,] w= 0 


L’équation (14) n'a pas d’infini au point wo = à pour lequel = et n° 
s’annullent. En effet, la règle de l’Hospital donne en ce point une vraie 


ae )s a "Nabi BP, ‘à La : 


se peut se calculer de 15 en 15° et analyse harmonique en donne la 


valeur en série de Fourier. 1 = 
Appelons 9 le coefficient de w dans l’équation du déplacement normal. 
à 2) . fe a 
Pour le déterminer calculons x = arc tg At Sime wpleas x de 15° en 15° ; 
dx 


P.o, Pog, P, ont été déterminés dans la première partie ; aie s’obtient par 


dérivation de l’analyse harmonique de x. 
Les coefficients de © développé en série de Fourier sous la forme 


(16) © — 0, + 96,0, + 20,0, + + 
sont publiés dans le tableau IL. 


§ 5. — Détermination de Vexposant caractéristique 


L’équation (15) étant une équation linéaire à coefficient périodique, elle 
admet des solutions de la forme 
w= 6% (w, + Du, 0, + 2w,C, + +.) 


u déduit de (14) par quadrature est de la même forme et, par conséquent, 
&, et n,. Donc & est l’exposant caractéristique de l'orbite périodique 
correspondante. La solution banale =, = =’, n, = n' est donnée par w = 0 


} 
il 
| 


pad 


On e2n two 


=) Wen el2n+v)iwe , 


——_- 
_ En introduisant dans les a ee on a pour chaque e@”+v)iwo 


+o 


+ 2n)* Wan — be Oon— 21 War = 0. 


[=— 


L’élimination des w donne un déterminant infini A(v)=0, dont les 
termes sont ; 
Ann (04 9) — 0, Am.n=— Onn (ma ni). 


En divisant chaque ligne du déterminant précédent par le terme de la 


_ diagonale principale, on construit un nouveau déterminant A a) dont 


les termes sont 


we ome à Om—n 
Daun le D — (v + In)? — 6 (m An). 


Ce déterminant converge, car le produit des termes de la diagonale 


principale est 1, et la somme des termes hors de la diagonale principale 


converge absolument. 


On voit aisément que À, (v) est une fonction périodique paire de période : 


égale à 2, qu’elle est une fonction analytique de v excepté à ces pôles 


v=— 7 + V6, et qu’elle tend vers l’unité quand v tend vers linfini par - 


valeurs imaginaires. 
Appliquant le théorème de Liouville, on a 


4 + k[cotg m (6: + v) + cotg x (86, — v)] = A,(v) = 0 
Si nous connaissons A, pour une valeur quelconque attribuée à v 


k peut être déterminé. 
Posons v = 0, nous obtenons l’équation A, (v) — 0 sous la forme 


sin? i 
9 


sin? = V0. 


= A, (a). 


re rx an a 
17 _ Soient n, n',n”… des entiers écrits en grandeur Aer ) 
et écrivons, pour simplifier [n] = (2n)? — 8 | 


' 


ea Ho 4 1—= 5-85 


> Tn] [n] [n’] | a | Re. 
gy Vin nn On Se Ont. ES On" On —n" } 
ra [n] [n°1 [n°1] sh + Gyn Opt mOn On —n'r S 
onan % aa fs Mes AFIN AU. 
À Ca] Le] n° Tu Ge] De] Ca] D — 
42 “ay Onn Op" nn ®,1v—n’ | 


| ae comets Up" Sniv—n" Ov" 
Ca] [7'] [n”| [n°1 En] 
Qn —n Op Faces 


> fn] Ca] [n°1 [n°1 [n°] [nY] 


où les 2 s'étendent pour tous les n de — 2 à + 2. 
Il est possible d'écrire ces sommations sous forme finie. 


Isolons par exemple le terme 6%, il a pour coefficient 


+e 1 
9%, Ca] [n+ k] 


Posons 9 = 46° et définissons les constantes A, B, C, D, telles que 


| aes B C D 
2 [n] Fan Ly ora C= du 4.0. Sce | 


Mais nous savons que 


| 
Dors De Dee eso, | 


d’où 


C3 ui DA 
ÿ 1 nl T cotg + V 8 


“ [n] (n EX] + 16 À + B+C<+D)r cotg T0 = ie Re: NEC) 


Z = OL 3 7 . be sf + A . , | 7 a … 
4 | a ras En 
we uvons FER de même toutes les sommations simples. L in | 
Pour les sémmations doubles, par exemple, — 
‘ + i@ 


Ti a = k] a Gr = ke 


; bétons an’, et ‘3 et procédons comme précédemment en considé: 
_rant d’abord 0 comme constante, | 


a B C 
î pie qua; ig Cr” ai he Sia ae | 
4 F G 

sal wary <5 a aco + O—n—{4" É On+k"+k ny ae oh ss 


Les coefficients A, B, C, D, E, F, G, H contiendront k’’. Nous devrons 
_-faire la somme de ces coefficients pour k’’ variant de 0 à + exceplen: 
ik’ =0,k"’ =k. Si nous faisons cette somme de 0 à +20, nous obtenons 
une expression en cotg x V6, de laquelle nous devons soustraire les ter- | 

mes à négliger. ER 
Cette méthode peut s'étendre à une sommation de n’importe quel ordre 
_ de multiplicité. 
Nous nous sommes arrêtés dans l’expression du développement aux 
termes non négligeables à la cinquième décimale dans le cas le plus 
- défavorable. À 


; 03 04 
| 1— 6, On 


3 03 04 7 (7—305) 82 04 06 A oie 
9(1-—0.) (4—8) | (1—80) Te re 0)(9—80) 


meotg (m V6) 41 2 


3 FRET re VE, Oe so ary 2 
one V6 dr Le 0) ne ) “ 

Tt cotg (n V/ 80) 62 04 
rcolg G V6 ) LM eV = at 4 = ages “a a i e 4(1—80)(4—860) 
- AS — 2 of 
5 t ey) 1 HU 3 Os 94 
q LEONE (eve 6 10 tea, aE aera 
: 1 9 0 nr cotg (7 V6) 

EL + + on V8 6 


fl 2 4 4) 2 DAS pm ane Pry 
get 921 1-0 Te) BG.) 4-8.) 9-8, 36, 
Les résultats de ces calculs sont publiés sous la notation 9, dans le 
tableau IV. 


v= =P" re 


nl qe, 4 
Ss pepo fea teceme + * 
s ’ "pa CRE ts t "i : 


ou Temailre revisant la Ener = Se oot Se ee m aon 
laquelle existe encore une orbite périodique, avait ¢ bt 
= = = 0, 788.545. Les éléments de cette orbite sont publiés dans la premi 
colonne du tableau II. ‘ 2: 

L’orbite périodique limite a son exposant caractéristique nul. L’équa- 
tion en w de l'orbite limite admet donc une solution périodique. Grac 
à cette propriété, nous avons déterminé Pexpression de Porbite périodi 
limite, exacte a six décimales. 


Supposons le terme non périodique dea ae 0, nous pouvons chose 


w de telle façon que ce terme soit égal à one 
La solution des équations aux variations devient 


{=.= 02+ 2,4+2,0,+ 2,0, --- 
(17) INS 
(n= on + 4,5, + LS, + uss 
| La variation de y, pourra être considérée du second ordre par rapport 
à la variation des autres éléments de l'orbite périodique, lorsque l’on 
passe d’une orbite à l’autre dans le voisinage de lorbite limite. Au 
premier ordre £(w + Aw) —£(w), où w est écrit pour la valeur limite, 
est solution des équations aux varialions 
E(w + Aw) — E(w) = [x, + Ar, + (x, + Ax,) ae Au) o...] 
—[x, + a, cos2wo... if 
=~ o Aw?’ + Ax, + Az,C, + Ax,C, + --- 
On voit par là que 


d dx, 
(18) = hey, 2, = — = 


de même pour les termes de À. 


‘ hed, Qw 
Si le terme non périodique de Dia 0, w est extrémé à l'orbite limite. 


Ce n’est pas le cas dans ce problème. 
Si nous soulignons d’une barre les éléments d’une orbite approchée de 


lorbite limite, la véritable orbite périodique limite s’écrit à des erreurs 
du second ordre près 


n= n + Aun, 


eT 


a—#" ae pes — 208 - 2h SO Pos) +w Pb Pa Bat. Pat =e LP, Fou 


LOT 


i. 
En nous bornant à une précision moindre, nous pourrons trouver une ae 
: notion périodique de ie 


at) Tach 


Le déplacement normal limite pourra s’écrire 
W=w+eC, + EC, … 


où €, €, --- sont des petites corrections qui peuvent être déterminées en 
même temps que Aw en introduisant l’expression de w et de © dans 
w"+ 0 (E&,n)w—0 
= Nous devons ensuite corriger l’orbite 2, n en tant qu’orbite périodique 
— de période w + Aw. Pour ce, nous corrigeons les coefficients du dévelop- 


pement et déterminons la variation de a= A6 que nous avions négli- 
i 1 


_ gée au premier ordre. + 
Ecrivons 
AS 2A, xe” — P,, = AP, 

ee — zAu,S;, NP, = APy 
; a? + \? __9Pp—9AP 
E Ces différences satisfont les équations 
j a ” 
4 os 4 Aa = P,Az+ PA + AP + Te o 
| Er 12 d 
a ‘uw? do Goo Ah = PAT AP AN Al 
uh 
¥ 


L’identification des coefficients des C;, S, donne les corrections. 


ve 


ns 


Ve eae de autres par er aux Wag charade aoa SW, 


première relation donne un résidu qui rend compte de la précision 
À l'orbite limite génératrice. 


_ L’orbite de paramètre y, = 0,788.541 donne un résidu égal à - 


sixième décimale, c’est-à-dire de l’ordre des erreurs sur la lek second 

de l’orbite périodique £ et n. Or, © est déterminé avec la même précision — 
que les dérivées troisièmes de = et n. Le résidu étant connu en fonction | 
des coefficients de ©, la limite des orbites périodiques est connue avec LE 


même précision que Tes orbites périodiques elles-mêmes. 


V1 | 964.953 
(ly 491.329 
fly —16.200 
4 9 
Ug 

by | —829.900 
bp | 145.644 
Un —2.349 
bg 62 
bg —2 
bio 

bie |. 

Cy 119.152 
C3 —812 
C5 3 
Cy 

Co 

8) | —637.536 
LE 13.440 
04 102 
0 

Qs 

IT] 


929.898 


~ 
474.535 


— 36.678 
27 


—918.031 
251.744 
—10.381 
417 

—19 


—187.400 
—2.622 
21 


— 600.939 
29.623 
496 
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TABLEAU Ill 
895.050 860.846 
446.169 399.371 

—63.677 — 99.674 
80 176 
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11 
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dw d 
Qw, 4.278 He | 805 Lae 
2w, «3119 i 10 di 
QWs 25 aes 
10 2 di 
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‘Théorie a mouvement gyroscopique des proje ¢ 


= Les états d’instabilité virtuelle Bers 
Les indéterminations 
Les perturbations initiales 
PAR à 


M. ROBERT D’ADHÉMAR 


Priirace. — J’ai constamment essayé de faire usage des travaux 
antérieurs ; j’ai tenté de justifier certaines hypothèses faites antérieure- 
ment. Mais, plus d’une fois, j’ai dû, réflexion faite, repousser des concep- 
tions généralement admises et des modes d’approximation qui paraissaient 
solides. 

Le développement de cette théorie est une suite de tâtonnements et cela 
n’a rien de surprenant, car le projectile se meut dans un milieu fluide “ 
que nous connaissons mal. 

Dans ce Mémoire, je rappellerai rapidement ce qui a été fait antérieure- 
ment. J’étudierai un minimum, voisin de zéro, de l’écart, dans un 
mouvement de révolution, pour compléter et coufirmer mes études. 
précédentes sur les valeurs nulles et presque nulles de lécart. 

Je reprendrai l’étude des perturbations initiales, parce que je me 
demande si certaines questions ont été bien posées. J’essaierai d’apporter 
plus de précision sur des questions dont mon exposé antérieur était 
défectueux. 

Le mouvement gyroscopique du projectile est difficile à saisir et je me 
demande si le vocabulaire de cette science est bien fixé, et si le même" 
terme n’est pas employé, parfois, avec des sens bien diferent: 7 


1. INTRODUCTION. 


Je rappelle le sens des notations. Nous désignons par G le centre de 
gravité du projectile. L’axe Gz coincide avec l'axe de figure, du côté du 
culot. L’axe Gx est perpendiculaire sur le plan de résistance, et l’axe Gy 
est situé dans le plan de résistance. Je suppose la rayure à gauche et le 


à rr \ 
> . rn 
¢ , 


contraire de celui 


composantes du moment résultant des forces extérieures sont 


M, M,, M. 


d'inertie transversal. Posons : v= A: B. 

Ce rapport, pour les projectiles allongés, peut être compris entre 
+. 
Pa 10 , | 

__ Soit Q la vitesse de rotation imprimée par les rayures. Nous avons 
constamment à faire usage de la vitesse réduite : 


0 w = vQ. 


Désignons par a le calibre, en mètres, et par v la vitesse, sur la trajec- — 


_ toire, en mètres par seconde. 
; Représentons par 7 la tangente de l’angle d’inclinaison des rayures, et 
- par Vo la vitesse initiale, à la bouche. Nous aurons : 
Vo 
a 


(2) Q = 9j 


Je désigne par 8 l’angle d’inclinaison de la tangente de la lrajectoire 
vraie, et par T l’angle d’inclinaison de la tangente de la trajectoire plane 
correspondante, qui nous sert d’approximation. Je désigne par n langle 
de dérivation. Je représente toujours par w’ la dérivée de « par rapport 
au temps ¢. Nous avons : 

dt ; COST 


8) ts 


Je représente souvent | t’| par J. 
Le coefficient balistique est désigné par c, et la masse est désignée par m. 
Le poids du projectile sera : K = mg. 
- Le poids de l’obus est toujours exprimé en kilogrammes. 

Nous avons, à Paris : g = 9,81. 

J’emploie souvent les notations suivantes : Soient deux nombres 
positifs x et y. Si y est pelit par rapport à x, j’écrirai : 


y = B(x). 


Il en résulte : 
TEYTZX. 


Si æ et y sont du même ordre, j'écrirai : 
a 00) ou bien ga Ota) 


a Les composantes de la rotation instantanée du’ projectile sont : p, q, p. — | 
Je désigne par A le moment d'inertie axial, et par B le moment 


i 
a 


Cette poussée es appliquée au point C, sur l'axe de figur 
Posons : GG = 1. Le point C est le centre de poussée. 


mr 


Breet: 


Nous pouvons dire, approximativement, que le point C correspond au 
milieu de l’ogive. Au point de vue théorique, il est inutile de préciser 
davantage. 

La poussée R est inclinée sur l’axe de figure ; soit y l’angle correspon- 
dant (Fig. 1). 

On appelle écart langle de la tangente de la trajectoire et de l’axe de 
figure, que je désigne par d. Pour conserver toujours le même sens 
d'orientation, je suppose essentiellement : 


(6) Rare 


4 
7 
} i 
; 
if 
“4 


po 


ù n'est pas une variable algébrique. J’ai donné des explications, à ce 
sujet, dans Mon Mémoire précédent. 

Le premier moment Mz est le moment de R, par rapport à G, et ce - 
moment définit le demi-axe Gz. Nous avons : 


(7) M.= Risiny. 
Nous posons : 
(8) sin y = kd; 
k est le coefficient d’obliquité. Si l’angle y est assez petit, nous avons : 
EM Tey i Lap 
(9) : Sin OLD 


Le paramétre k n’est pas parfaitement connu, et il varie probablement 
sur la trajectoire, surtout dans le cas du tir courbe. Lorsque le tir est de 
plein fouet, nous pouvons supposer que le coefficient de k est à peu près 
constant. 

Les paramètres k et / ne sont pas hien connus, et il faut, dans les études 
théoriques, leur donner des valeurs numériques vraisemblables ('). 

D’après ce qui précède, nous écrirons : 


(10) M.=Rklo>0. 
Rkl ners 4e: 
Posons : P = AQ? cette fonction joue un rôle important. 
Posons : Q — 7 c’est la Fonction de Mayevski, dont le rôle est 
fondamental. 
Il est souvent commode d’écrire la formule (10) sous la forme suivante : 
(11) M: = A2Pd= BwPd. 


On fait, par la formule (4), une séparation trés brutale des variables. 
Il est clair que, les Tables numériques de F(v) ayant été établies pour 
certains types de projectiles, ne seraient pas valables pour des projectiles 
trop nettement différents de ceux qui ont été utilisés pour les mesures. ~ 

Nous avons aussi des Tables numériques très précieuses, qui donnent 
le degré de résistance n, défini par la relation : 

4 dF 
(2) F dv 

Notons bien les faits. Le coefficient balistique n’est pas parfaitement 
connu. La trajectoire vraie n’est pas parfaitement connue, et nous lui 
Re TE Pe oe al PEI RE 


(1) J'ai essayé de donner quelques indications, à ce sujet, dans mon Mémoire : Note 
sur les trois couples fondamentaux et sur l’approximation. Ann. de la Société scient. de 


Bruxelles, 1935. 
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substituons la trajectoire plane. Nous ne connaissons pas bien 
variations des paramètres k et L. l'a Ftimss 


* 


Que faut-il conclure ? C’est que la théorie du mouvement gyroscopique — 


des projectiles tournants peut être faite dans ses grandes lignes, que l’on 
peut saisir l’aspect général du phénomène, mais que nous ne pouvons 
préciser le détail. 

Cherchons à voir l’ensemble ; c’est déjà assez difficile. 

Nous avons défini le premier moment Que dire du deuxième moment 
M,, qui correspond à des forces dont la résultante est perpendiculaire 
sur le plan de résistance ? 

Ces forces correspondent, d’après M. Esclangon, aux frottements 
latéraux ('). Autrefois, on admettait que le deuxième moment provient 
de l’effet Magnus, mais c’est assez vague. 

On admettait aussi que M, a un signe constant et contient 0 en facteur. 
C'est une hypothèse faite à priori, qui n’est pas toujours vérifiée. 

Posons : 0 — M, : Mz. Le deuxième moment a souvent été négligé, 
sans explication, ou bien l’on affirmait que le rapport © est nécessairement 
petit par rapport à l’unité. C’est une hypothèse sans fondement. 

La notion des frottements latéraux de M. Esclangon est la plus naturelle 
et la plus souple, et il faut lui conserver ce caractère de souplesse, en 
raison de notre ignorance. Je crois qu’il serait téméraire de fixer, 
à priori, l'expression des frottements latéraux. 

En particulier, nous ne connaissons pas les cavitations, qui entourent 
le culot, et qui modifient profondément la nature des frottements latéraux 
sur certaines zones de la surface de l’obus. Nous ne savons rien au sujet 
de ces zones, sur lesquelles les frottements latéraux sont presque inexis- 
tants. Il est dangereux d’imposer à priori telle ou telle forme à des 
phénomènes trop peu connus, à des phénomènes dont la nature peut 
varier fortement, d’un point à l’autre de la paroi. 

En bref, pour le moment, n’affirmons rien, au sujet de My. C’est 
prudent. 

Dans les travaux anciens, on assiste à un déploiement extraordinaire 
de quadratures, au sujet du calcul de la poussée R. Tous ces calculs ne 
nous ont rien appris, parce qu’ils reposent sur des hypothèses concernant 
les actions de l’air sur le projectile. Et nous ne savons absolument rien 
de la valeur de ces hypothèses. Si l’on fait maintenant des calculs du 
même genre, au sujet de l'expression de My, il est possible que cela 
représente encore un effort inutile. 

Considérons le troisième moment M. Il est bien vraisemblable que 
nous pouvons le négliger, après les perturbations initiales, Nous faisons 


i) Ernest ESCLANGON, Sur le mouvement des projectiles autour de leur centre de 
gravité. Mémorial de l'Artillerie Française, 1927. 
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plusieurs approximations de simplification ; celle-ci n’est pas la plus dan- 
_ gereuse. D’ailteurs, si nous connaissions Mk, il serait, en général, facile 
d’en tenir compte. | 

Je dois faire une remarque, au sujet de l'effet de vide à l'arrière du 
projectile, effet qui a sa répercussion sur chacun des trois moments. La 
. cavitation arrière est représentée schématiquement par une poche W 
_ (hachurée dans la figure 1). Sans préciser davantage, nous pouvons dire 

que cette poche a une action plus forte, si la vitesse v est plus grande. Sur 

une trajectoire de plein fouet, nous savons que, pour la partie utile de la 
trajectoire, v, P, Q sont décroissants lorsque ¢ croît. Donc, sur la partie 
utile de la trajectoire, l’effet de la poche W est toujours décroissant. Pour 
une même valeur de à, la valeur de + sera moindre, si v est plus grand. Si 
. cette notion intuitive est juste, le rapport y : d et le coefficient k seront 
croissants, avec ¢, sur la partie utile de la trajectoire. 
Et cela aurait pour effet de régulariser les variations de P et de Q, sur 
.” la trajectoire. 

Il est possible d’admettre que, pour un bon projectile, le coefficient k et 
la vitesse v ont, sur une trajectoire, des variations de sens contraires. 
M. Burzio, en se plaçant à un autre point de vue, a proposé cette hypothèse, 
qui parait érès judicieuse, et que l’on peut adopter provisoirement (’). 

Si les paramètres k et / ont des variations faibles et de sens contraires, 
cela permettra de supposer que le produit k/ est à peu près constant. Dans 
certaines questions, fort heureusement, les valeurs numériques de k et de 
l n’interviennent pas. 
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2. LES ÉQUATIONS DU MOUVEMENT GYROSCOPIQUE. 


La figure 2 représente une trajectoire. Le sens d’orientation étant direct, 
l’axe O’Z est dirigé vers l'arrière. Je désigne par S le sommet, par D le 
point de wilesse minimum, par E le point de vitesse maximum, par C le 
point de chute. Si la dérivation était nulle, le point C’ serait sur l’axe O’X. 

Je désigne par M le point où la fonction Q de Mayevski est minimum. 
Dans le cas du tir de plein fouet, les points M, D, E sont au delà du point 
de chute. Dans le cas du tir courbe, les points D et M peuvent être très 
voisins du sommet. 

Je suppose que le point D existe ; c’est ce qui a lieu, en général. Cepen- 
dant M. Esclangon a fait remarquer qu’il existe certaines trajectoires, sur 
lesquelles v est toujours décroissant (?). 


(1) F. Burzto, Recherches sur la rotation des projectiles, Mémorial de VArtillerie 
Francaise, 1927, page 820. Le us 

(2) ERNEST ESCLANGON, L’Acoustique des canons et des projectiles, Mémorial de 
l'Artillerie Française, 1925. page 56. 


ur une bonne approx! 
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Faut-il prolonger la trajectoire indéfiniment ? On le fait parfois, et cela 
paraît naturel, mais, en même temps, cela nous met dans l’embarras. 


eggs 


RiGeros 


Prolonger la trajectoire à U’infini, c’est supprimer la terre. Que dire 
alors de l'accélération de la pesanteur ? Que dire de la densité de l’air ? Je 
ne crois pas que ce prolongement puisse nous éclairer beaucoup. 
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FIG. a: 


sée R, par rapport à G. L’observateur placé sur Gx voit GT tourner dans 
le sens direct, pour venir sur GA, à étant essentiellement positif (ou nul). 
Soit un point O, sur GT, tel que l’on ait : 
GO = 1. 


En ce point O, nous prenons un plan OHN perpendiculaire sur GT. Le 
LVIU, I 4 


D”: tete a 


rase us PE proj ojectile Bi ve le isin ag référence ( 
oot a. En suivant le mouvement de a, nous voyons le mouvement g 
copique (Fig.5). Gs 

Au point G, nous prenons un plan Gay, perpendiculaire sur la tangente; f 
Gz, est parallèle à OH, et Gy, est parallèle à ON. L’axe Gz, est confondu — : 
avec la tangente, en sens inverse du mouvement de G. Ce triédre Ga,y,z, 
nous servira pour définir les angles d’Euler (Fig. 4). Le demi-axe Gz est KE 
confondu avec l’axe de figure, dans le sens du culot. Le triédre GUVz est 
lié au projectile ; c’est le tièdre lié. Je dirai que le triédre G xyz est le 
trièdre libre. Nous passons de GUVz à Gayz, en posant : = 0. 

J’ai établi les formules fondamentales suivantes (°) : 


ve p = d'+ 8’ cos y + n/ cos 8 sin y, 
: (43) { q = y’ sin d —n sin 8 sin d+ cos d (— 8 sin y+ n'cos 8 cos y), 
p = p'+w’ cos d—n' sin 8 cos d—sin d (—6' sin y+n/ cos 8 cos y). 


On suppose, en général, que la trajectoire plane est une approximation 
suffisante de la trajectoire vraie. Je dirai que cela constitue l’approxima- 
tion TP. Dans ces conditions, nous remplacerons 8 par + et 8’ par tT’ dans 
les formules ci-dessus. 

On suppose, en général, que à est assez pelit, ce qui signifie que l’on 
considère exclusivement des tirs assez réguliers. Nous pouvons, dans ces 
conditions, substituer à à sin à, et À à cos à. Je dirai que cela constitue 
Papproximation SC. 

Les formules (13) prennent alors la forme suivante : 


{ p = 0 —J cos y+ n' cost. sin y, 
(14) { q=w'd+ J sin y +n cost. cos y— dn sint, 
p= p'+ y’—dJ sin y— 1 sin T— dn' cost. cos w. 


Posons : p=p'+r. | 
Les composantes de la rotation instantanée du trièdre lié sont : p, q, p 
Les composantes de la rotation instantanée du trièdre libre sont: p, q, r. 
Je m’occupe toujours de la rayure à gauche, et, par suite, la dérivation _ 
normale est à gauche du plan de tir, pour l’observateur placé derrièrele 
canon. 
La dérivation est toujours négligée dans la première approximation. 


(1) Sur le mouvement pendulaire des projectiles stables, Mémorial de l'Artillerie 
Française,1933. Voir aussi mon Mémoire: Les équations du mouvement gyroscopique 


+ a Dérivation des projectiles tournants, Ann. de la Société scientifique de Bruxelles, 
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pose : n =0 (Fig. 3). Je dirai que cela constitue Papproximation DN. 2 
1 $ Les équations (14) prennent alors la forme suivante : 


L p=d —Jcosy, 
8) qg=yd+Jsiny, 

pt p= 9'-+ y’— dJ siny, 
et nous avons : r = y'— bd] sin y. 


L’approximation DN est acceptable pour un tir de plein fouet, avec un 
_ bon projectile. Elle ne le serait pas avec un projectile défectueux. Dans le 
cas du tir courbe, il faudrait voir si cette approximation ne modifie pas 
trop les caractères du phénomène gyroscopique. 

Enfin, je fais souvent une approximation, en substituant w’ à r. Je dirai 
que cela constitue l’approximation A. Il faudra vérifier, avec soin, la vali- EL 
. -dité de cette approximation A. a 
7 Considérons les dérivées par rapport au temps, w” et 7’. Substituons 

wy” ar’. Cela constitue approximation A’. Il n’est pas certain, à priori, 
; que l’on puisse faire l’approximation A’. Cela devra être examiné de près. 
J 
; 


J’emploie toujours les équations du théorème du moment cinétique 
sous la forme suivante : 
‘ 


dp SOS 
dia q(vp —r)= B M , 


4 dq 4 ae Ÿ 
(16) ar BWP") = My, (v = 5) 
> Frs A M:. 


C’est la forme adoptée par M. Sugot, suggérée par M. Saintillan (’). 

Je dirai que les équations (13) et (16) constituent le système GY. 

Si nous négligeons la dérivation, les équations (15) et (16) constituent 
le système GYS. 


3. LES PERTURBATIONS INITIALES. 


Il faut signaler les perturbations initiales, qui ont été ignorées autrefois, 
et qui, actuellement, ne peuvent être omises (*). Nous avons, après le 
départ, une période pendant laquelle l'écart, partant de la valeur zéro, a 
une suite de maxima et de minima très rapprochés. Cest la période 
d’instabilité initiale, ou d’irrégularité initiale. 


(1) G. Sucor, Balistique Extérieure Théorique, Mémorial de l’Artillerie Française, 
1927. | 

(2) R. H. FowLer, E. G. GaLLoP, etc. Aérodynamique d’un projectile tournant, trad. 
Mémorial de l Artillerie Française, I, fasc. 2 et 3. 
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de prends les conditions les plus favorables, c’est-à-dire que Je suppose 


le projectile bon et le départ normal. (à siyup i 

Nous pouvons avoir une régularisation rapide de y'’, avec un écart pas 
trop voisin de zéro. Nous avons alors un mouvement du point a à peu près 
circulaire et uniforme, de centre O. Je dirai que cela constitue la presta- 
bilisation. À 

Je me garde bien d'affirmer que la prestabilisation a toujours lieu et 
qu’elle a toujours lieu à une petite distance de la bouche, par exemple 
‘après une ou deux secondes. Mais si l’on n’a pas une régularisation de ce 
genre, il paraît difficile de faire la théorie du mouvement gyroscopique. 
Il est possible que ce mouvement soit désordonné sur toute la trajectoire. 

Les gaz qui sortent de la bouche, en jets violents, sont très rapidement 
dilués dans l’air, de sorte que, si le projectile est bon, les perturbations 
initiales doivent être rapidement amorties, ce qui correspond à une régu- 
larisation rapide des phénomènes sonores. 

_ D’une manière générale, quand on étudie le mouvement gyroscopique, 
on suppose le projectile assez bon, les perturbations initiales assez faibles, 
Pamortissement des phénomènes sonores assez rapide, la dispersion faible. 

En un mot, nous nous occupons d’un tir assez satisfaisant. Sinon, il 
serait absurde, par exemple, de dire que M: est négligeable, que / est à 
peu près constant, que k varie peu et varie régulièrement. 

Je donne ici une simple indication sur la prestabilisation, cette question 
devant être reprise ultérieurement. Cette notion a un fondement physique. 
En supposant que la prestabilisation a lieu à une petite distance de la 
bouche, je suppose que les phénomènes sonores sont devenus réguliers à 
une petite distance de la bouche. 

Il faut consulter, à ce sujet, les études de M. E. Esclangon, que j'ai 
citées antérieurement. 


4. L’ÉQUATION DIFFÉRENTIELLE DE LA PRECESSION EULÉRIENNE. 


Il faut étudier maintenant la première équation du théorème du 
moment cinétique, que l’on peut écrire sous la forme d’une équation 
différentielle par rapport à y. Il s’agit ici de la précession eulérienne, et 
Je précise, parce que nous aurons aussi à considérer la précession 
balistique, qui est autre chose. 


Lorsque la composante p est variable, nous pouvons poser : 
(17) EG) = vo(t), 
E est la vitesse réduite 


Nous pouvons, sans inconvénient, supposer que p est constant, après la 
période d’irrégularité initiale, et peut-être même pendant cette période. 


2 réduite est ane w, oe la “formule (). 
rrions avoir, aussi bien, après la période | 


. 2 P= Q— = const. 


tant assez voisin deQ. hy 
Dans ce > cas, la vitesse réduite serait : 


we 4 


Ww hm 


730 


il faudrait Eatitoes Q, à Q dans l’expression de Pp. 
~ Padopterai, en général, la relation (18), mais cela n’a rien d’essentiel. 
Dans ces BOO NODS, les trois équations du moment cinétique se réduisent 


à deux : 


PU 


0) a eae pe ee ar 
en a Pur = EM. 


Et si j’emploie l’approximation A, ces deux nÉUAÈRTS prennent. la 
forme suivante : 


pay Sn RM 
: : 
(23): - st re — plu —y)=< 5 M. | 


Il faut voir ce que natin cette approximation A. 
Nous avons, d’après ce qui précède : 


(24) 


Je suppose v assez grand et 5 assez pelit ; il en résulte que J est très 
petit par rapport à 1. : 

Supposons que la trajectoire soit de plein fouet. L’angle de départ esi 
inférieur à 40° ou 12°, et la vitesse initiale est supérieure à 500 ms ou 
600 ms. La valeur absolue de t ne dépasse pas 15° ou 18°, par exemple. 
Nous pouvons admettre que l’on a, au point de chute : 


v > 250 ms. 
Supposons : v=—9250 et cost —1. Nous aurons approximativement : 


sin y =i 2 . 


_=0( 11 et ety 


42 I ‘ +, hits 4 ya ] 
# Les civconstasceseont tn oe différentes pour une 
Si nous avons, au sommet : v—100, nous sy its | 


PRIT 1=0 (6) 


Sil’ sty est de 15°, nous avons approximativement : 


Cela donne : 


C’est encore un nombre petit par rapport a 1. 

Mais, dans le cas du tir courbe, tout devient plus difficile ; en particu- 
lier, nous ne savons pas bien si nous pouvons négliger la dérivation. 

ll est prudent d’étudier d’abord le tir de plein fouet et de voir ensuite 
si les résultats subsistent quand on passe au tir courbe. 

Donc je considère un tir de plein fouet et j’emploie les équations GYS. 
En outre, j’admets que la relation (18) est vérifiée, et je dirai que ceci 
constitue AD DTO TUE B. 

Revenons à la relation (24) et considérons le rapport : 

W—Y ~ pz 
ae a w— yl 
Ce rapport a pour valeur 4, silona. 


ù =, (0 ou bien sin w = 0, 
Ge rapport est très voisin de À, si l’on a : 
y= 0, 
ou bien si l’on a : 
(25) ly] = (uw). 


Getle condition (25) joue un rôle très important dans la Théorie 
Élémentaire. 


Supposons maintenant : wy’ >0. 


Ve 


te 


Je suppose d’abord que y' est un très grand nombre ; je suppose, par. 


exemple. que l’on a : 


w = B(y’). = 
I] est clair que, dans ce cas, le rapport ci-dessus est très voisin de 4. 
Je suppose maintenant que l’on a : | w — w' | —0(1). 


pu ae : W — y 
Le rapport He sera négligeable, et, par suite, le rapport er 
sera très voisin de 1. | 
Je suppose toujours: wy’ >0, et je suppose que la relation ci-dessous 
est vérifiée : 


(26) y'~w, 


aE ora Pa 
Dans ce cas, nous ne savons pas, à priori, si le rapport , est très 
— y 


W 


~ petit par rapport à 4. Il faudra examiner chaque cas particulier. 


Mais il faut remarquer que j’emploie souvent l’équation (20) pour 


calculer y’. Si nous pouvons savoir indirectement que y’ est voisin de w, 
nous n’avons pas besoin de faire usage de l’équation (22). 

Je suppose donc que la relation (26) n’est pas vérifiée. Nous pouvons 
substituer l’équation (22) à équation (20). 


En posant, pour simplifier la notation: w' =z, je puis mettre 


l'équation (22) sous la forme suivante : 
(27) vd — wa d + wPd — wJ sin y — à” + J'cos y —0. 

Il faut faire usage de la dérivée, par rapport à ¢, de l’équation (20) ; 
nous verrons que c’est une relation fondamentale. Cette dérivée peut 
être écrite sous la forme ci-dessous : 


d’p ; 
(28) - TE ae oD) — QT = KPO). 


Si nous prenions la dérivée de l’équation approchée (22), nous aurions 
léquation : 
d’p 
2 mn 
Et cette équation serait obtenue, aussi bien, en prenant la dérivée de 
l’équation (27), c’est-à-dire en écrivant : 
(30) | xd — waxd + wP!’—wicosp.c¢—d”—J'siny.xv| _ 0 
| + 2xx'd — wzx'o + wP'd — wJ' sin y + J cos w f 
Posons : D = x? — wæ + wP ; ® est connu, si æ (ou wp’) est connu. 
L’équation (30) prend la forme suivante : 
[| Do + ad (Qe — w) + w Pd UF 
| —wicosy.x—J'siny (x + w) + J” cos y — 0” ( 


+ q(w—2)— qa'= w(Pèy. 


(30) 


des ene Bi sth JRA ini ue 


on is ae avons substitué — qa’ a.—q(a' +S). Done, en dériv 
: À à t, les deux membres de l'équation Re Le nous 
tons le terme : 


cl qS. | 
Pour faire un calcul correct, il faut ajouter ce terme dans le premier : 
membre de l’équation (29), ou bien il faut ajouter le terme : 


via ice 


dans le second membre de l’équation (29). 
Par suite, l’équation (28), provenant de la dérivation de pt ra ur 
peut être mise sous la forme ci-dessous : AE 


(32) D 6 + w’d (2y'— w) + w P'o + qS 
—wJcosw.w’—J'sinw. (w’+ w) + J” cos y — 0” | 
L’équation (32) me permettra d’obtenir l’expression générale de My. 
il faut étudier le terme S et, dans ce but, il faut substituer a l'équa- 


tion (24) l’équation complète (Avis). ‘ 
Je néglige toujours la dérivation. Nous avons : 


(24bis) r = y cos d — sin d .Jsinw. 
Prenons la dérivée par rapport à { ; nous avons l’équation : 
(Au) r= ( y" cos À — cos à . d'J sin y 
— y" sin à . à — sin à (J' sin y + J cos y . w'). 


L’écart d étant assez petit, nous pouvons écrire : 


(31) sy ts 
et nous avons : 
(33) S = — qd — d (J' sin y + J cos w . y’). 


La relation (31) représente une approximation, qui est d’autant plus 
satisfaisante que à est plus voisin de zéro. 


5. EXPRESSION GÉNÉRALE DU DEUXIÈME MOMENT. 


La trajectoire est donnée ; nous connaissons donc, à l’instant 4, les 
valeursidé "Pia i gs 


Supposons le point a donné; c’est-à-dire supposons et y connus. 


‘ 
L 
| 


… - 
A Le : 
ee” Ç - ER 
tito: > 


oF sons, en outre, y’ connu et UT à eee Qe Dans'ces FD ; 
conditions, les valeurs de ® et de qS sont connues. = 
Il faut remarquer que la dérivée première 0’ et que la Ain troisième: 
ym ne figurent pas dans l’équation (27). Je considère actuellement le 
mouvement syroscopique le plus général, et non le mouvement gyrosco- 
pique élémentaire, qui est un mouvement gyroscopique particulier. Dans 
ces conditions, nous devons regarder d' et 0” comme étant des paramètres 
indépendants. 
.  Mettons l'équation (32) sous la forme ci-dessous de 1" à 
BH), Ayu) + Ad + B— =O, 1 
enposant: er las iis dhe 
4 ae Ge | 
er De | w P'o — wJ cos y . y'— J'sinvw. ET + w) + eos Ww 
| — qù(J' sin y + J cos y . y’). 
Les valeurs de A et de B sont connues, si 5, w et w’ sont connus. 
Posons : W= Ad’ + B— 0”, 
C’est une fonction linéaire en 0’ et en 0’”, si nous regardons 0, w, w 


comme étant des constantes données. Il y a deny cas a distinguer, relati- 
: vement a la valeur de w’. Je suppose que l’on ait : 
Ë 
+ 
j 


D (35) ye wee 


Dans ce cas, le produit w’’d est déterminé par l’équation (34), et ce 
produit est une fonction linéaire de à’ et de à”. 
Je suppose maintenant que l’on ait : 


_ (36) 2 y — w — 0. 
Si l’on a, pour cette valeur de y’: 


D (57) W— 0, 
alors w’’d est indéterminé; ce produit peut prendre une valeur finie 


quelconque. Cette question devra être discutée. Si, au contraire, la : 
relation (36) étant vérifiée, on a, pour cette valeur de y’: 


(38) W 40, 
~ alors la valeur de y’’d sera infinie. Nous verrons que ce fait est dépourvu 
de sens, au point de vue de l'interprétation physique. 


Voyons maintenant quelle est l'expression générale de My. 
Je simplifierai la notation en posant : 


M; = BX, et M, = B Y. 
Nous avons, en faisant l’approximation À : 
Op (Er 


at 7 zi ee sehen QUI wn ds de 
“teat Hitter dann 20 ted ja æ fet) RR 


1e di se à con be osent 


wi ap others 


ds prete bé win (ge upiqon2ov ea HS 
(40) Y = y'd + 9@y’ — + wl cosy + d'a | se 


En faisant usage de l’équation (34), nous avons : 

ut d’(2w’ — w) + wJ cosy + J'siny. ; 
È D’après ce qui précède, Y peut être bien déterminé ; c’est le cas général. 
Il peut exister une indétermination, mais seulement pour une valeur 
assez grande de w’. Cela sera examiné. 

Enfin, supposons vérifiées les relations (36) et (38), alors la PS: 
de Y serait infinie. Je ne vois pas la possibilité de l’interprétation physique 
de ce fait. 

En faisant usage de la deuxième équation du moment cinétique, et des 

expressions de p, q, p, que j’ai établies, j'obtiens l'expression générale du 
deuxième moment, sans connaître le détail des phénomènes physiques. 
Il est dangereux de vouloir fixer l'expression du deuxième moment 
d’après des hypothèses sur la nature des phénomènes physiques, 
hypothèses peu sûres, en général. _ 

L'interprétation physique de la formule (41) serait probablement fort 
difficile. On a beaucoup étudié leffet d’un obstacle fixe sur l’écoulement 
d’un fluide (*). Mais il y a ici wn élément nouveau ; il faut considérer un 
obstacle tournant très rapidement autour de son axe de figure, et cette 
rotation doit modifier profondément les phénomènes. 

L'interprétation physique est plus facile lorsque le mouvement gyrosco- 
pique est un mouvement particulier, le mouvement élémentaire. 


(M) Y=— 


6. ETUDE D'UN MINIMUM PRESQUE NUL DE L’ECART. 


Dans une étude précédente (*), j’ai examiné un minimum presque nul de 
l'écart, dans un mouvement d’oscillation. Je vais étudier la même ques- 
tion, dans un mouvement de révolution. Les circonstances ne sont pas 
identiques. 

Dans la figure 6, la courbe A, A À, représente un mouvement d’oscilla- 
tion. Au point A, l'écart est minimum et l’on a : , 


d =), 0” > 0, sin w = + 1, cos p= 0. 


(') A. Focx, Introduction à la Mécanique des fluides. Collection Armand Colin, 1932. 
(?) Théorie du mouvement gyroscopique des projectiles. Le mouvement gyrosco pique 


élémentaire. Essai d'interprétation physique. Annales de la Société scientifique de 
Bruxelles, 1937. 


20 — usd + UP + ul = 0, 


one 


Fic. 6. 


Je fais d’abord une remarque, au sujet de l’approximation A. Pour avoir L 
l'équation complète (la dérivation étant toujours négligée), il faut ajouter 
au premier membre de (43) le terme : 

qd = (yo — J) Jo. 

Il faut donc substituer (— wd + Jd’) y’ a — wdy’. 

L’altération est insignifiante, en général, et surtout lorsque l'écart est 
petit. 

De même il faut substituer (wP — J*)d à wPd. 


L’altération est insignifiante, en raison des valeurs numériques de 
w, Pe Je 


di , ju 
Par exemple, si Q reste supérieur à 10, ou à 5, nous voyons que w fe 


¥ 7 


4 
est trés grand par rapport a J. 


Nous pouvons donc conserver |’équation (43). Nous voyons qu’elle n’a 
pas exactement la même forme que l'équation (42). 


* : premiers termes, dans l'équation (43), est voisine de zéro, comme (x - 


a6) RARE su nag (=f) watt sinber alsa 
À jui bi Lo fa | ty i My wh nojttil 


Nous voyons bien, en effet, que, dans ce cas, la_ somme des 


Pe 

Pour ne pas alourdir mon exposé, je laisse de côté cette D A > 
particulière ; elle sera examinée, si. c’est nécessaire. Et j’emploie léqua; 
tion (43), pour les mouvements de révolution, et Fees (42), ss les 
mouvements d’oscillation. 

Je vais définir maintenant le minimum presque fon de l’écart. C’est une 
valeur de à, très voisine de zéro. Mais l'écart est toujours petit, si le pro- 
jectile a une bonne tenue. Par exemple, si l’on a un angle de 10°, on aura 


approximativement : d — . . 
1 


Si langle est de 5°, on a approximativement : d = ay 

I] faut donc préciser un peu le sens du terme « presque nul », dans cette 
théorie. Et n’oublions pas que, dans les approximations de la Balistique, 
les mots ont souvent un sens un peu élastique ; cela est imposé par la 


nature des choses. Il s’agit, par exemple, de savoir si un nombre € est 
négligeable par rapport à l’unité, ce que j'écris sous cette forme : 


(44) e=6 (1). 


Prenons : e — 100 ? la condition est très bien remplie. 


, nous pouvons admettre que la condition est à peu 


Prenons : €= = , la condition est assez bien remplie. 
1 
Prenons : € = = 


prés remplie, mais ri n’est pas très satisfaisant. 


Si lon a : e = + , on peut dire que la condition (44) n’est pas remplie. 


Certaines approximations sont trés bonnes: d’autres sont seulement 
acceptables. 
J’écris maintenant les équations de Mayevski : 


(MA) ra SA 
wdo= Pd — J sin w. 


Considérons la seconde. J’ai établi qu’elle a une valeur, comme appro- 
ximation, si l’on a : 


(25) dose (4 


W 


; 


[+ ile 


D A te 


érons un minimum de lé écart, représenté nae le eh d'u iF 
mor vement d’oscillation. Nous avons : sin w=+ 1. 
LA seconde équation (MA) prend la forme ci-dessous : 


 ®) y= Pe PE 

F ? Je prends une valeur de l’écart assez petite, par exemple : 
iy : à à | ee. a nf 

BE i, À 
. Nous aurons : A Zz : E 
(47) : | y= P— uw] <0. 


Cette inégalité permet de voir que le passage de la révolution à l’oscil- 
lation ne saurait avoir lieu tout près de l’origine de la trajectoire. 
En effet, la relation (47) peut être mise sous cette forme : 


Q<w. 


Et nous pouvons avoir : w = 150, ou 200. 

Au contraire, la valeur de Q, dans la région de l’origine de la trajectoire, 
pour un tir de plein fouet, est beaucoup plus grande. 

Nous pouvons écrire : 

(48) lw |=wJ—P<wJ, 
ou bien : — | 

La valeur numérique de J est trés inférieure à 1. 

Par suite, la seconde équation (MA) est valable, par un minimum de 
l'écart défini par (46). 

Considérons maintenant un minimum de l'écart, représenté par B 
(Fig. 6), dans un mouvement de révolulion, et supposons à défini par la 
relation (46). 

Nous avons : sin y= —1. 

La seconde équation (MA) prend la forme suivante : 


(49) Y=P+S=P+u >0. 


Nous pouvons écrire : 


Be UD 0). 


aaa i. mtg 
1 Pee Sh Qu AE ay ane 
sie WD ot A i) ale neler 
Ww mt 
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ns ce cas, la seconde équation. (MA) sera Bor ave tS 
Et, a fortiori, elle sera valable si l’on a : Tie 


Es 


un minimum a ù défini par (46). ee 
_ Par suite, lorsque je dirai que à est presque nul, cela pie dire 


~ 


que 0 est beaucoup plus petit que _ 


Il est inutile d’insister ; la suite fera comprendre nettement le sens du 
terme adopté. | 

Il faut maintenant reprendre les équations (MA), dont l’interprétation 
a toujours été superficielle, et même défectueuse. J’ai établi qu'il faut 
comprendre les équations de Mayevski de la manière suivante (*). 

La première équation (MA) est équivalente à celle-ci : 


(50) p=0, 


En particulier, ona: p=p’=0. 
La seconde équation (MA) est une relation approchée : 


(51) q= Po, 


qui a un sens bien clair, si la condition (25) est remplie, et qui n’a aucun 
sens dans le cas contraire. 

Dans l'étude du mouvement gyroscopique, nous ne pouvons séparer les 
deux équations (MA). C’est leur ensemble qui joue, car ces deux équations 
proviennent de la décomposition artificielle d’une seule équation, l’'équa- 
tion fondamentale (20), qui prend la forme (22), si l’on applique 
V'approximation A. 

J'ai employé l’équation (51) ; maintenant je vais faire usage de l’équa- 
tion (50). 

Nous avons, en général, la dérivation étant négligée, l’expression 
suivante : 


(52) | “P= y"—Veosy+Jsiny.y’. 


(') Voir ma Note: Le moment d’impulsion initiale et Pinclinaison des rayures. Les 


périodes d’instabilité virtuelle des projectiles. Comptes rendus des séances de l'Académie 
des Sciences, 26 octobre 1936. 


oo a eee = : > Fi 
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Si l'expression de 
aa Gp ip 
6 x : dt #0, 
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er. 
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la dérivée d” est telle que l’on ait: 


Mur an" 


_Pensemble des équations (MA) n’est pas valable, et cela, quelle que soit 
la valeur de y’. Ceci nous donne immédiatement un résultat. Considérons 
un mouvement d’oscillation, et un minimum de l’écart, représenté par le 
point A. Il est inutile de supposer que à est presque nul; la valeur de d 
est quelconque. Je fais hypothèse suivante : 


4 H : = wJ » 


_ Si le mouvement gyroscopique est régi par les équations de Mayevski, 
_ nous aurons : p’=0, et, par suite: —- 


EG) ‘ dise Joie olen!) a 
EE Nous aurions donc la relation : 
(53) = [ys | ==... 


Il y a contradiction, car la condition (25) n’est pas remplie. 
Je fais maintenant l’hypothèse suivante : 


ai ”=(w+h)J, bo 0". 
Si le mouvement gyroscopique est régi par les équations (MA), nous 

_aurons la relation : 
D (56) lwl/=wt+h>w. 

[Il y a encore contradiction. 

Considérons un mouvement de révolution et un minimum de l'écart, 
représenté par le point B. Si le mouvement gyroscopique est régi par les 
équations (MA), nous avons la relation : 


(ory lear ue 0. 
Faisons l’hypothèse H ; il y a contradiction. Faisons l'hypothèse H'; 
il y a encore contradiction. | | 
Je conclus. Pour un minimum de à, situé sur l’axe H'OH, quelle que 
soit la valeur de à, si l’on à, en ce point : 
(58) a Sw) « 


il est certain que le mouvement gyroscopique ne peut être régi, à cet 
instant, par les équations (MA). 


SE D eh 


Pour achever cette discussion, revenons aux équations (42) et (43). 
Les solutions de l’équation (42) sont données par la formule : 


EUX 
(59) 2e w + \/ut— sup ote. 


(est la solution négative qui nous intéresse, car le mouvement est une 
oscillation. 

Nous faisons maintenant Vhypothése Beant 
(60) > | 0<d<ul, | 

Il est clair que, si la valeur de à est très voisine de zéro, nous aurons 
une solution négative, dont la valeur absolue sera grande. Par suite, la — 
condition (25) ne sera pas remplie, et le mouvement gyroscopique ne . 
sera pas régi par les équations (MA). 

Les solutions de l’équation (43) sont données par la formule : 


a 


D pee 


(61) : 2œ— w + \/u' — 4 wP — 4 


Nous faisons la même hypothèse, exprimée par sr (60). 

Il est clair que, si la valeur de à est très voisine de zéro, les racines 
auront des valeurs imaginaires. Par suite, dans ces conditions, le 
mouvement gyroscopique est impossible. 

La conclusion est simple. Si nous avons un minimum de l’écart, situé 
sur l’axe H’OH, et si la valeur de à est très voisine de zéro (je dis: 
presque nulle), dans ces conditions, a cet instant, le mouvement gyrosco- 
pique du projectile ne saurait étre régi par les équations de Mayevski. 

J’ai donc trouvé, tout simplement et naturellement, des circonstances 
telles que les équations (MA) soient certainement inapplicables. 
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Sur la définition générale de l'entropie 
et les bases de la thermodynamique chimique 


PAR 


]. — INTRODUCTION. — REMARQUES SUR LE THEOREME DE CLAUSIUS 


Soit dQ une quantité de chaleur élémentaire absorbée par un systéme 


à la température T. La somme 
AQ 
js 


étendue à toutes les quantités de chaleur élémentaires absorbées par un 
système qui subit une transformation a été appelée par Clausius la valeur 
de la transformation. Conformément à l’usage qui tend à prévaloir 
aujourd’hui, nous l’appellerons la chaleur réduile. Étudiant au point de 
vue du principe de Carnot les modifications réelles des systèmes fermés, 
dans le temps croissant, Clausius sut établir le théorème fondamental 
qui porte son nom : lorsqu'un système parcourt un cycle, la chaleur 
réduite absorbée au cours de ce cycle ne saurait être positive. Cependant 
la démonstration de Clausius manque de généralité, car elle implique 
cette hypothèse restrictive qu’à chaque instant le système est caractérisé 
par une température uniforme. 

Même sous cette forme restreinte le théorème de Clausius se montra 
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tellement fécond que trés peu d’auteurs se souciérent de le génér 


Parmi eux il convient de citer H. Poincaré ('), qui consacre un lo ; 
chapitre de sa thermodynamique à « l'extension du théorème de Clausius» — 
et rappelle les idées déjà émises à ce sujet par J. Bertrand (?). a 

Envisageant donc le cas où la température varie, non seulement dun 
instant à l’autre au même point, mais aussi d’un point à l’autre au même — 
instant, Poincaré conclut que le théorème de Clausius reste vrai, si la 
somme est étendue d’abord à tous les systèmes élémentaires, dont il 
suppose que le système total peut être considéré comme composé, et 
ensuite, pour chaque système, à tous les intervalles de temps dont se 
compose la durée du cycle. La question se pose de savoir si dQ doit alors 
représenter la chaleur totale qu’un système élémentaire reçoit pendant 
l'intervalle de temps dt, ou seulement la chaleur cédée au système 
élémentaire par des sources de chaleur extérieures au système total. 
Poincaré accepte cette ambiguïté en démontrant que l’une ou l’autre de 
ces deux interprétations laisse subsister le théorème de Clausius. Mais la 
démonstration indiquée par Poincaré est basée sur la décomposition du 
système total en systèmes élémentaires infinitésimaux, auxquels le théorème 
de Clausius, sous sa forme restreinte, serait directement applicable. Or, une 
telle décomposition n’est pas toujours possible, comme il est facile de s’en 
assurer en songeant à un mélange de plusieurs gaz dans lequel la concen- 
tration varie d’un point à l’autre. La diffusion plus ou moins active, dont 
un tel système est nécessairement le siège, empêche de considérer les 
éléments de volume ou de masse comme des systèmes partiels auxquels 
le théorème de Clausius serait directement applicable et pour lesquels les 
raisonnements de Poincaré garderaient leur clarté et leur rigueur. 

Il est toutefois possible d’éviter la décomposition du système étudié en 
systèmes élémentaires, et de démontrer d’une manière parfaitement 
générale (*) que, dans l’évaluation de la chaleur réduite qui forme l’objet 
du théorème de Clausius généralisé, la somme doit être étendue à tous 
les intervalles de temps et, pour chaque intervalle, à tous les éléments de 
surface du système. L’élément de l'intégrale | T ‘dQ est donc la quantité 


de chaleur qui traverse, vers l’intérieur du système, l'élément de surface 
do pendant l'intervalle de temps dt, rapportée à la température instanta- 

née, de cet élément de surface. Le théorème de Clausius pour un cycle : 
quelconque peut donc s’exprimer comme suit : : 


ty n—1 d’Q 
i qu oN do dt < 0. | 


ty 


(1) H. Poincar, Thermodynamique (rédaction de J. Blondin) :2™° édit., p. 219 à 258 x 
Paris (1923). 7 ’ 

(@) J. BERTRAND, Thermodynamique, p. 265 à 273, Paris (1887). 

(5) Ann. Soc. scient. Brux., t. 45 (1 part.), p. 273 à 279 (1926). 
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Il est bien entendu d’ailleurs, conformément à la remarque judicieuse 
. formulée par Poincaré, qu'il est prudent de limiter les raisonnements 
thermodynamiques basés sur le théorème de Clausius à la chaleur trans- 
mise par voie de conduction thermique, à l’exclusion de la chaleur 
rayonnée. 

- Au début de cette introduction, nous avons tenu à rappeler en détail la 
signification du théorème de Clausius généralisé, non seulement parce que 
l'application de ce théorème jouera un rôle essentiel dans la suite de ce 
travail, mais aussi parce que ces développements préliminaires permettent 
de mieux saisir les données du problème que nous nous proposons 
- d'aborder. 
À Lorsqu'il est possible de faire passer d’une manière réversible un fi 
* certain système A d’un état (1) en un état (2), tous les auteurs s’accordent 
_ _à définir la variation d’entropie S, —S, en Videntifiant avec la chaleur 
_ réduite absorbée dans la transformation réversible (1) — (2). Il résulte 
» en effet du théorème de Clausius (grâce à un raisonnement classique) que, 

_ dans toutes les modifications réversibles qui peuvent amener le système 
de l’état (1) en l’état (2), la chaleur réduite présente la même valeur. On 
peut dès lors se demander comment il convient de définir l’entropie du 
système A, lorsque ce dernier ne se trouve pas dans un état d'équilibre. 
En effet, une modification réversible n’est définie que comme une suite 

= continue d'états dont chacun est infiniment voisin d’un état d'équilibre. 

= Un état très différent d’un état d'équilibre ne pourra donc en aucun cas 
être considéré comme le terme d’une transformation réversible. Par 
conséquent, l’entropie du système exige une définition plus explicite 
lorsque le système n’est pas en équilibre. Comme exemple d’un système 
très éloigné de l'équilibre nous pouvons nous figurer la masse gazeuse qui 
s'échappe de la bouche d’un canon immédiatement après le départ 
de l’obus. Dans cette masse la température et la pression varient d’un 
point à l’autre, il s’y propage des ondes, il s’y produit des réactions 
chimiques avec des vitesses localement très différentes et la compose 
tion chimique de ‘chaque élément de volume s’y modifie au surplus 
par l’effet de la diffusion. 

Dans la troisième partie de ce travail, nous formalerons la définition de 
l’entropie qui paraît convenir pour l’étude du cas général, et nous en 
examinerons quelques conséquences. Toutefois, afin de pouvoir donner 
à cette définition de l’entropie toute la précision désirable, nous aurons 
à faire usage de certaines propriétés thermodynamiques bien connues des 
systèmes homogènes, propriétés dont la présentation sous une forme 
appropriée, avec certains commentaires, fera l’objet de la deuxième partie. 
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- at les sons ouverts. Les premiers peuvent se inc 
d’une surface fermée qui se déforme et se dilate ou se contra 
ia modification considérée, de telle manière qu'aucune particul 
ne lat traverse. Le contenu d’une surface fermée quelconque, ne satisf 
pas à cette condition, constitué par définition un système ouvert. 

Un système homogène (au sens strict du mot), qu’il soit En vast nie 
fermé, se caractérise par l’uniformité vérifiée à chaque instant par la 
distribution de la température et de la concentration de chaque espèce 
chimique intervenant dans la composition du système. 

Considérons d’abord un système fermé, homogène, qui ne soit le siège 
d'aucune réaction chimique. L’accroissement de son énergie est définie 
; comme la somme de la chaleur absorbée et du travail des forces exté- 
an rieures appliquées au systéme, chaleur et travail étant exprimés dans 
k une unité commune. Parmi les forces extérieures, les seules que nous 
ayons à considérer sont celles dont l’effet est équivalent à une pression 
normale et uniforme exercée sur la surface du système, de sorte que le 
travail peut s'exprimer par () 


2 
G—- {Pa 
if 


en désignant par P la pression, et par dV l’accroissement élémentaire du 
volume. Puisqu’en vertu du premier principe de la thermodynamique 
Pénergie U est une fonction de l’état du système et qu’en l’absence de 
réaction chimique cet état est complétement défini par le volume V et la 
température absolue T, nous pouvons poser 


(4) dU aoe Pave ak dT + rh dV. 


il est essentiel de remarquer que l’hypothèse d’une température et d’une 
composition chimique constamment uniformes exclut toute transformation 
irréversible dans un système dépourvu de mobilité chimique. Ainsi que 
nous y avons insisté dans une étude antérieure (*), toute modification 
thermiquement ou mécaniquement irréversible suppose, en effet, l’exis- 
tence d’un gradient de température ou de pression. 

En vertu de la définition habituelle de l’entropie et de l'équation (1), 
nous pouvons écrire pour un système fermé, homogène, en absence de 


| 


(1) Voir par ex. Ann. Soc. scient. Brux., t. 44 (I part.) p. 508- 510 (1925). 
(?) Ann. Soc. scient. Brux., t. 56 (série B) ip. 15-34 (1936). 
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| réaction chimique, 
6 s = @_ dU + Pav + (8 
DO woe LC) HG) seb 


J D’autre part, comme l’entropie est une fonction de l’état du système, 


-@ as =(Sr) a+ (5) av 
; La comparaison de (2) et de (3) fournit les égalités classiques : 
D as 1 (aU 
FE Gr), = 7 Gr), 
7 b as 
| @æ) ES = 5 Gy). diced 


Supposons maintenant que le système fermé, a considéré soit 
- au contraire le siège d’une réaction chimique. Le théorème de Clausius 
nous oblige alors à écrire 


(5) ds Be 


puisque toute transformation chimique s’effectuant réellement dans un 
systéme constitue un processus irréversible. En effet, le fait méme que 
le système quitte son état instantané E pour un état E’ montre bien que 
l’état E n’était pas un état d'équilibre et en différait d’autant plus que le 
passage de l’état E en l’état E’ s'effectue plus vivement. L’irréversibilité 
des réactions chimiques réelles est d’ailleurs illustrée d’une manière frap- 
pante par le principe de Le Chatelier, que l’on peut déduire très directe- 
ment du principe de Carnot (*) et dont une des conséquences est précisé- 
ment l'impossibilité de la réaction inverse. 
Nous écrirons l’équation (5) sous la forme 


(5bis) Tds = dQ + dQ’ 


ou dQ’ = ET est une grandeur essentiellement positive. 

La grandeur e?, que Duhem représente souvent par dw, a déjà été con- 
sidérée par Clausius sous le nom de transformation non compensée. De là 
le nom de chaleur non compensée que l’on donne ordinairement au produit 
eT. On peut d’ailleurs critiquer cette terminologie, que nous adopterons 
cependant pour nous conformer à l’usage général. En effet, au second 
membre de l’équation (5), le premier terme peut être considéré comme 
représentant cette partie de l’accroissement d’entropie qui correspondrait 
à absorption physiquement réversible de la chaleur dQ à la température 


(1) Ann. Soc. scient. Brux., T. 56 (série B), p. 20-24 (1936). 
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d’entropie dans le cas d’un phénomèt 
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chaleur réellement absorbée par le système. eon RUE. 
Quoi qu’il en soit de cette question de mots, la considération explicite e 
de la chaleur non compensée s’est montrée extrêmement utile dans l'étude 
des phénomènes irréversibles. Déjà Duhem (!) a systématiquement appliqué 
cette méthode dans plusieurs paragraphes remarquables de son Traité de 
Mécanique Chimique, etil a développénotamment l’expression de lachaleur 
non compensée (représentée chez lui par Tdw) en fonction des variables 
dont dépend en général l’état d’un système. Récemment, dans le domaine 
particulier de la thermodynamique chimique, Th. De Donder (*), suivant 
la même voie, a réussi, non seulement à retrouver les théorèmes fonda- 
- k mentaux, mais parfois à les généraliser et à en développer les nombreuses 
; conséquences. 
Cherchons maintenant les expressions de dU et de dS dans le cas d’un 

système ouvert, homogène, dépourvu de mobilité chimique. Soient n,, 

n,, .… Ng, --- les nombres des molécules-grammes des différentes espèces 

chimiques dont le système se compose à l’instant considéré. Nous compa- 

rerons d’abord entre eux deux systèmes fermés, A et A’, ayant la même 
température et le même volume et contenant les mêmes composants chi- 

miques. Ceux-ci seront représentés par les mêmes nombres de molécules- 

grammes, à l’exception du composant (i), dont le système A contiendra n 
molécules-grammes tandis que le système A’ en contiendra n, = ng + dna. 

Si les deux systèmes subissent la même variation de température, T,-T,, et 

la même variation de volume, V,-V,, sans que leur composition chimique 

se modifie par suite d'une réaction, les chaleurs absorbées et les travaux 

des forces extérieures ne seront pas rigoureusement les mêmes. En vertu 
de la définition de la variation d'énergie nous pouvons donc poser : | 


Puisque la différence entre A et A‘ se réduit à un excès de dn, molécules- 
grammes du corps (i) dans le système A’, nous écrirons 
U; x. U; == Dy +p dye D) 


La différentielle sera considérée comme la différence de deux différen- 
tielles : 


t 


he P. Dunem. Traité élémentaire de Mécanique Chimique. Tome J, pp. 81, 82, pp. 85 
86 et surtout pp. 155 à 158. Paris (1897). 

(2) Les importants travaux thermodynamiques de Tu. DE DONDER se trouvent résumés 
dans l'ouvrage : [Affinité (rédact. P. Van Rysselberghe), Paris, 1936. 
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dU, du, 
dni dni 


Chacune de ces deux différentielles représente l’accroissement que subit 
Pénergie du système A, lorsque, sans modifier la température (T, ou T,) 
de ce système, ni son volume (V, ou V,), on y ajoute un nombre de molé- 
cules-grammes du corps (i) égal à dni. 

Les considérations précédentes, relatives aux systèmes fermés A et A’, 
permettent de comprendre comment on se trouve amené trés naturelle- 
ment, en assimilant un système ouvert à un système fermé, différent d’un 
instant à l’autre, à définir l’énergie du système ouvert par celle du système 
fermé qui lui équivaut à l'instant considéré. L’énergie d’un système ouvert 
est done une fonction des variables T, V, n,, n,, … ni, .… La dérivée de 
énergie U (TVn,n, ... nz …) par rapport à n; est, par définition, l'énergie 


(6) 


dn; — 


dni 


# 


d(U, —U,)— 


. moléculaire partielle du composant (i) et sera représentée par 


dU 
(7 Ui= ( Feil 
) ON / rvn: 
où l’indice n’ signifie que tous les nombres de molécules-grammes sont 
fixes à l’exception de n;. En convenant que l'indice n signifiera que toutes 
les variables n,, n,, ..…. ni, --- sont fixes nous écrirons 


ne rate (nan (am 
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On peut montrer d’une maniére analogue que pour un systéme ouvert, 
homogène, dépourvu du mobilité chimique, l’entropie également peut être 
considérée comme fonction de T, V, n,, n,, +. ns. ---. Nous définirons 
Ventropie moléculaire partielle du composant (i) par lidentité 
(Ts) | = (À). 

TVR 


Oni 


et nous écrirons : 
@ = (H) are (H+ 2(38), du 


Mais il est évident que si la composition chimique, en principe variable, 
d’un systéme ouvert est, en fait, maintenue constante, comme. le signifie 
indice n, ce système se comporte comme un système fermé. Dès lors les 
égalités (42) et (45) entraînent les égalités 


(104) Capo) 


(105) A ARAT À 
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A4) TdS—dU—PdV= TE Sd — E Uidi = 
Dans la dernière équation, W; désigne la différence U; — 
| l'énergie libre moléculaire partielle, qui est, par définition, la 
partielle par rapport 4; de la fonction d'état que Helmholtz a int 
sous le nom d'énergie libre du système : W=U—TS. : 
L’équation (11) peut aussi se mettre sous la forme : 


GA) dE =— PAV — SAT + Ei Vids, one 
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_ puisque 
=. dY = dU — TdS — SaT . AU RES 
Nous voici préparés à enfin traiter le cas important d’un système fermé, 
homogène, au sein duquel se déroule une réaction chimique. 
Cette réaction sera représentée par l’équation chimique 


NM, -+ NM, +. NUM; ++ NoM$ +... 


où M,, M,, … Mi, Mg --. sont les masses moléculaires et N,, No, … Na, Na, + 
des nombres entiers appelés coefficients stoechiométriques. En convenant 
de considérer comme négatifs les coefficients relatifs aux corps consommés 
par la réaction, nous écrirons l’équation chimique générale sous la forme 


(12) EZN;:M: = 0 
Le rapport 


_ Mi 
c= À 


sera appelé la concentration du corps (i) dans le système. Cette concentra- 
tion varie au cours du temps, du fait de la réaction chimique. Nous défini- 
rons la vitesse de réaction par l'identité 


(13) wa 1 dû 


En vertu de la convention que nous avons faite quant au signe de Ni, 
la vitesse de réaction sera la même quel que soit le corps (i) considéré, 
pourvu que ce corps participe à la réaction. 

Ces définitions étant posées, nous sommes en mesure d'interpréter de 
la manière suivante l’artifice de raisonnement véritablement génial (') 
grâce auquel J. W. Gibbs réussit à constituer une thermodynamique 
chimique. 


(') et d’ailleurs implicite. 


ee 


Un système homogène té qui est le siège d’une réaction de vitesse 
_W, peut être assimilé, au point de vue thermodynamique, à un système 
héMogelle ouvert, dakourvh de mobilité chimique, o chaque composant 
serait introduit au cours du temps à la vitesse 


(AA) dni _ ww 


dt 
Cette proposition présente un tel degré d’évidence physique que, 


lorsqu’elle est formulée comme elle vient de l’être, nul ne songera sans 


doute à l’appeler un postulat. Elle nous autorise à appliquer au cas du 
système fermé, homogène, en voie de transformation chimique, l’équa- 
tion (11) que nous avons déduite plus haut en pensant seulement aux 
systèmes ouverts, chimiquement inertes. Bien entendu dn; sera mainte- 
nant remplacé par sa valeur tirée de (14). Nous obtenons ainsi 


(45) TdS — dU — PdV= —VWEYiNedt. 


Si le système était le siège de plusieurs réactions exprimées par des 
équations linéairement indépendantes, on aurait 


(16) TdS — dU—PadV——V EE W; ViNadt . 


Bien que considéré comme ouvert, pour tenir compte de la variation de 
composition déterminée par la réaction chimique, le système dont nous 
parlons pour le moment n’en reste pas moins un système réellement 
fermé. Nous pouvons, par conséquent, lui appliquer l’équation (5°) en 
remarquant que, pour tout système fermé 


dQ = dU + PdV. 
En combinant entre elles les équations (5P) et (16) nous pouvons écrire 


(17) a = VE SW VIN. 


Considérant le cas d’une réaction unique, il importe de faire remarquer 


que la grandeur 

4 dQ’ 
(18) N'a gee ul 
est une fonction de l’état instantané du système. 

Il est évident, en effet, que si l’énergie et l’entropie d’un système 
ouvert, où aucune réaction chimique n’a lieu, sont fonctions des variables 
T, V,n,,n,,... Ns, --», Chacune des dérivées partielles ne dépendra que 
des valeurs de ces variables à l'instant considéré, ou, plus brièvement, 
sera une fonction de l’état instantané du système. Or, à chaque instant, 
un système fermé en voie de transformation chimique peut, s selon Vingé- 
nieuse méthode de Gibbs, étre assimilé, au point de vue de son énergie 


ou de pi 
pendant l'intervalle de temps dt. En effet, grâce aux équations (10), 
termes en dT et en dV qui intervenaient dans les expressions (8) et (9) de 
dU et de dS se sont trouvés éliminés en formant la différence (11). Quant 
à l'augmentation de pression dP, elle dépend à la fois de l’état du système 
et des différentielles dT, dV et Wdt. nats 
Le fait que la somme ZX WiN; est une fonction de l’état du système 


n'avait certes pas échappé a J. W. Gibbs. Cet auteur désigne par wi le 


rapport 
Sat 


qu'il appelle le potentiel chimique du composant (i). En se basant sur les 
équations (8) et (9) et sur d’autres, analogues, qu’il est facile d'écrire en 
considérant U et S comme fonctions de n,, n,,... ni. .… et de certaines 
paires de variables autres que V, T, on peut démontrer d’une manière 
très directe (*) que 


Cowan (aby = Ce sh PV) ‘' Ge) (TE SA 
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_La somme Zu; M;dn; qui intervient dans les raisonnements de Gibbs 
s’identifie ainsi avec le second membre de notre équation (45). Or, il paraît 
évident que, tout au moins d’une manière implicite, les potentiels 
chimiques u; sont considérés par Gibbs comme des fonctions de l’état du 
système, c’est-à-dire comme des grandeurs entièrement indépendantes des 
variations de température et de volume qui peuvent se trouver associées 
à la réaction (*). = ; 

Il parait opportun de prévenir ici une confusion. Si, pour caractériser 


(1) Ann. Soc. scient. Brux., t. 56 (série B), pp. 26-28 (1936). 

(?) C'est ce qui explique la remarque suivante de H. Le Chatelier : « C’est là en réalité 
une hypothèse qui revient à attribuer à chaque corps une force chimique définie, comme 
il possède un poids défini. » (J. W. Gress : L'équilibre des systèmes chimiques, trad. Le 
Chatelier, p. 15, Paris, 1899). Il est douteux que Le Chatelier ait exactement compris sur 
ce point l’idée de Gibbs, car la comparaison entre le potentiel chimique et le poids est 
boiteuse à un point de vue essentiel. En parlant de poids, sans le distinguer explicitement 
de la masse et sans rappeler au lecteur que le poids d’un corps déterminé varierait d’une 
planète à l’autre, Le Chatelier semble bien vouloir évoquer l'idée d’une constante carac- 
téristique. Or, cette idée ne s'applique nullement au potentiel chimique, puisque ce 
dernier dépend, non seulement de la définition qualitative du système dont le corps 
considéré fait partie, mais même de l'état instantané de ce système. | 


l'état physique du système, nous avons recours à une paire de variables 
quelconques, que nous appellerons en général x et y, on peut se demander 
si la fonction de l’état du système que nous avons Jusqu'ici représentée ~ 
par Yi peut s’exprimer par 


(20) wel (Eye oe) ) 04e 
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On peut montrer (') qu’il en est effectivement ainsi, si l’on choisit 
comme variables, non seulement T et V, mais encore T et P, S et V ou S 
et P. Mais P. Van Rijsselberghe (?) a prouvé que dans d’autres cas (par 
exemple avec P et V) il n’existe aucune fonction F(xyn,n, ... ni...) satis- 
faisant à l’équation (20). En ce sens, il est exact d’affirmer que le choix de 
certaines variables exclut l’existence d’un potentiel thermodynamique 


correspondant. Mais quelles que soient les variables x et y employées avec 


M, M,... N;.. pour caractériser l’état du système, cet état reste défini 
11) 2 b) 


sans ambiguité. Les énergies libres moléculaires, qui sont des fonctions 
de l’état du système, restent donc également définies, même si on décide 
de les exprimer au moyen de x, y,n,, Ns, ... ni... Comme, d’autre part, 
tout changement d’état physique peut aussi bien se décrire par dx et dy 
(n,, Ny, ... Ni, ... restant constants) que par dT et dV, la somme 


—ZINiVi(zyn, ny... Ni--+) 


sera aussi bien indépendante de dx et de dy que la somme 


LA 
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était indépendante de dV et de dT. Dans l’esprit de la terminologie de 
Gibbs, cette somme devrait s’appeler le potentiel chimique du syslème 
relatif à la réaction, ou, plus simplement, le potentiel de réaction. 
Pour van t Hoff, qui s’inspirait directement de Helmholtz, la diminution 
d’énergie libre associée 4 une réaction constituait la mesure de cette 
tendance à réagir que les chimistes appelaient depuis longtemps « Paffi- 
nité » (°). En effet, ainsi que Helmholtz l'avait souligné, la diminution de 
énergie libre représente le travail maximum qu’un système, maintenu à 
température constante, est capable d’effectuer à l’encontre des forces 
extérieures en passant d’un certain état initial en un certain état final (et 
qu'il n’effectue d’ailleurs que si ce passage est réversible). Ce travail 


() Voir par exemple : Ann. Soc. scuent. Bruc., t. 56 (série B), p. 24 (1956). 

(2) C. R. Acad. Sc., Paris, 2 décembre 1935. . . a 

(3) Pour l'historique de la notion d’affinité, voir O. SACKUR : Die chemische Affinitat 
und ihre Messung. (Sammlung Wissenschaft), Braunschweig (1908), ainsi que : Lehrbuch 
der Thermochemie und Thermodynamik(2™° édit. par Cl. Simson), pp. 234 et suiv., 


Berlin (1928). 
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maximum peut servir à mesurer les réactions mécaniques que le système, 


dans certaines conditions, serait susceptible d’opposer aux forces exte- 
rieures. Cela étant, l’idée de Van tHoff se réduit à identifier ces réactions 
mécaniques avec les forces affinitaires. Si cette identification PROPRES 
légitime, ce ne sera que dans le cas particulier où l’étatinitial et | état final 
ne différent entre eux, ni par le volume, ni par la température, mais uni- 
quement par la composition chimique du système. Dans le cas général, la 
diminution élémentaire de l'énergie libre peut s'exprimer au moyen des 
équations (14%) et (14), dont la combinaison fournit la relation : 


— d¥ = SaT + PaV — (EN:Ÿ:) VW dt. 


Cette égalité fait distinguer dans la diminution élémentaire de l'énergie 
libre une partie due à la variation de température, une partie due à la 
variation de volume et une partie due à la réaction chimique. Il est clair 
que cette dernière partie seule, laquelle s'exprime par 


— a¥ — SdT — PdV = —(EN;V;) VW dé 


peut servir à la mesure de l’affinité. Encore faut-il qu’elle soit rapportée 
à la transformation chimique d’une quantité définie de matière. Le choix 
de cette transformation unitaire peut s’effectuer de plusieurs manières. 

Par exemple, on peut songer à mesurer laffinité par la diminution 
d'énergie libre qui correspond (à volume et température invariables) à la 
consommation ou à la formation de N; molécules-grammes de chaque 
composant. Puisque dni/Ni= VW dt, on parvient ainsi à la définition de 
Vaffinité qui a été récemment adoptée par Th. De Donder (’) : 


— d¥ — SAT — PdV 
VW dt 


Cette affinité, définie par De Donder, se confond, comme on voit, avec 
le potentiel réactionnel qui est à la base de la thermodynamique chimique 
de Gibbs. Elle est le résidu d’une sorte d'épuration mathématique de ce 
concept d’affinité, encore assez grossier, que l’on rencontre dans les écrits 
de Van Hoff. Il y a lieu seulement de noter, à un point de vue purement 
formel, que De Donder a préféré définir l’affinité en l’égalant au premier 
plutôt qu'au second membre de légalité (18). Nous verrons dans la 
troisième partie de cette étude que ce choix, tout en n’offrant qu’un 
intérêt formel dans le cas des systèmes uniformes, est cependant désavan- 
tageux, en ce sens qu'il exclut la généralisation de cette notion d’affinité 
dans le cas des systèmes non uniformes. 

Il existe encore une autre manière possible de définir la quantité de 
matière dont la transformation chimique s’accompagnera de cette dimi- 
nution d'énergie libre qu’on adoptera comme mesure d'affinité. On peut 


= LIN. 


() DE Donner : L’Affinité (rédact. P. Van Rijsselberghe), p. 19, Paris (1936). 


| 
} 


EME re ee 


4 

L 

# notamment considérer la quantité de matière transformée par unité de 
__ temps dans le cours de la réaction réelle. Cela revient à définir l’affinité 

" par la grandeur : . pee 
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Cest cette définition de l’afinité, sous la forme où elle est exprimée par 
le second membre de la dernière équation, qui nous paraît la plus adéquate 
et que nous adopterons. Cependant, pour ne pas trop nous écarter de la 
terminologie déjà consacrée par les nombreuses publications de De Donder, 
nous appellerons affinité thermodynamique l'affinité définie par cet auteur 
(c’est-à-dire le potentiel réactionnel de Gibbs) et nous appellerons affinité 
réelle le produit de l’affinité thermodynamique par la vitesse de réaction et 
par le volume. L’affinité réelle a donc les dimensions d’une puissance. Cest 

~ sous ce nom que cette grandeur est occasionnellement mentionnée dans les 
travaux de De Donder. 

Le principal inconvénient que nous verrions à accepter telle quelle la 
terminologie proposée par De Donder se dégage des réflexions suivantes. 
Si l’affinité, c’est-à-dire la grandeur par laquelle doit se mesurer la ten- 
dance du système à réagir, est une fonction de l’état «statique » du système, 
indépendante de la vitesse de réaction, cette dernière peut être nulle sans 
que l’affinité s’annule en même temps. C’est là au fond un paradoxe qui ne 
peut se soutenir qu’à la condition d'introduire l’idée d’une espèce de 
« frottement » capable d’enrayer la marche de la réaction en dépit des 
« véritables » forces affinitaires. Cette distinction entre les forces chimiques 
actives et certaines résistances passives remonte à Gibbs et a été fort accen- 
tuée par Duhem. Elle est inspirée par l’analogie avec les phénomènes 
mécaniques étudiés à l’échelle macroscopique et elle devait, à ce titre, 
plaire à l'ingénieur qu'était H. Le Chatelier (’). Cependant, au point de vue 
moderne, Vanalogie paraît bien artificielle. La réaction chimique est le 
résultat statistique d’un grand nombre de processus élémentaires, dont 
chacun met en jeu une ou quelques molécules. Or, rien de ce que nous 
savons au sujet de ces processus élémentaires, que leur aspect quantique 


() J W. GiBBs : Equilibre des systèmes chimiques, Traduct. Le Chatelier, pp. VII et 
VII, Paris (1899) : « Bien des chimistes méconnaissent encore le rôle des résistances pas- 
sives et sont amenés à contester les lois les plus certaines de la mécanique chimique ; ils 
se trouvent dans la situation d'ingénieurs qui voudraient appliquer à leurs machines les 
principes de la mécanique rationnelle en faisant systématiquement abstraction du frotte- 
ment ». Qu'il nous soit permis de faire observer à propos de cette remarque de Le Chate- 
lier que ce qui importe sans doute ce n’est pas tant de distinguer formellement les forces 
passives des forces actives, que simplement d’en tenir compte (fût-ce même sans les con- 
sidérer comme formant une catégorie spéciale). C’est précisément ce que nous proposons 
de faire, à la différence de De Donder, lequel distingue les forces de frottement, mais 
uniquement pour exclure leurs effets dans sa définition de l’aflinité. 


| semble à me rer ses re ign 
anata d’üne distinction entre une force tendant à prov 
réaction et une résistance antagoniste. + 

Il vaut donc mieux, en effet, définir comme affinité le une 
qui s’annule seulement avec la vitesse de réaction et qui cesse dès lor 
se présenter comme une caractéristique purement thermodynamique. 
système. 

Ce point de vue peut se justifier aussi bien en un langage un peu ditt 
rent. Tous les auteurs actuels seront d’accord, en principe, pour définir 
l’affinité par le degré d’irréversibilité propre au phénomène chimique con- 
sidéré. Nous trouvons naturel d'adopter comme mesure de cette irréversi- — 
bilité la chaleur non compensée débitée par unité de temps. Ainsi nous 
nous conformons strictement à la méthode originale de Clausius, qui étudie 
les modifications réelles des systèmes dans le temps croissant. De Donder 
a estimé au contraire qu’il convenait de faire abstraction du temps et de 
| rapporter la chaleur non compensée à l'étape franchie par la transforma- 

; tion chimique ('), c’est-à-dire à dni/N;. Cela revient à envisager l’irréver- 
sibilité d’une réaction virtuelle (c’est-à-dire simplement possible au point 
de vue stoechiométrique) au lieu de lirréversibilité de la réaction réelle. 
Au sens de De Donder, l’affinité est l’irréversibilité qu’une réaction mani- 
festerait si, par ailleurs, cette réaction pouvait avoir lieu. Le « par ailleurs» 
exprime précisément ce facteur affinitaire essentiel qui fait défaut dans la 
définition critiquée (?1. 

Un système où l’affinité réelle est nulle, doit nous apparaître comme étant 
véritablement en équilibre, car, pour un tel système, le degré d’irréversi- 
bilité se réduit à zéro (équat. (17)). Il est bien entendu que les états d’équi- 
libre peuvent être relatifs el qu’il importe de marquer la différence entre 
les états d'équilibre relatif (ou états de « faux équilibre » d’après Duhem) 
et l’état d'équilibre absolu. Ge dernier est déterminé par la condition 


ZN = 0 


Quand cette condition est satisfaite, l’irréversibilité dQ’/dt s’annule en 
vertu de (17), mais alors toute réaction chimique se trouve par le fait 


(1) De Donder désigne par dé le rapport dni/Ni, où Ni est considéré comme une gran- t 
deur algébrique qui sera négative toutes les fois qu’il s’agit d’un composant consommé : 
par la réaction. 

@) Nous n’avons développé ici nos remarques sur la définition de I’« affinité » que par 
égard pour le point de vue de De Donder. Quant 4 nous, nous ne verrions aucune diffi- 
culté à nous en tenir simplement à l'exposé original de Gibbs. Nulle part ce dernier ne 
parait avoir songé à établir une relation entre l’affinité et les fonctions thermodynamiques 
dont on peut bien penser qu’il n’ignorait pas les propriétés. C’est que rien n'empêche de 
considérer l'aspect quantitatif de la notion d’affinité comme s’exprimant par la simple 
définition de la vitesse de réaction. quitte d’ailleurs à constater que, pour le moment, 
cette grandeur ne nous est encore connue que d’une manière empirique. 
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; même exclue, puisque toute réaction chimique se poursuivant avec une 
… vitesse diftéreñte de zéro constituerait nécessairement un processus irré- 


 blement » nul, si l’on peut exprimer ainsi le fait que les deux facteurs 
. Wet N.Y; sont nuls séparément. L’état d’ équilibre absolu jouit du privi- 
| lege de la stabilité entendue dans le sens suivant: au voisinage d’un état 

4 . d'équilibre absolu le système ne peut évoluer que vers cet état. En effet, 

q si la vitesse de réaction est positive, le degré d’irréversibilité dQ’/dt est 


- nécessairement positif lui aussi. En vertu de (17) la somme totale XN;W, 


1 


a 
doit donc être négative. Cela revient à dire que la somme partielle relative 


_ aux corps engendrés par la réaction est plus petite en valeur absolue que 
4 Di somme partielle relative aux corps consommés. Par conséquent, la 
. somme totale ZN:WŸ; doit nécessairement se rapprocher de zéro par les 
| progrès même de la réaction; c’est dire que le système évolue vers son 
* état d'équilibre absolu. 
> En considérant comme des états d'équilibre véritables, quoique relatifs, 
. certains états que, dans la thermodynamique de Van t’Hoff, on n’aurait 
peut-être pas considérés comme des états d'équilibre, on évite une difficulté 
qui serait bien embarrassante. Soit un mélange d’hydrogène et d’azote 
> (8H, + N,) à la température ordinaire. Si nous déclarons que ce système 
n’est pas dans un état d'équilibre véritable (sous prétexte que l’affinité 
+ thermodynamique ne s’annule que pour un mélange contenant une fraction 
- prépondérante d'NH,), nous encourons une grave contradiction. Tous les 
traités de chimie physique décrivent, en effet, une suite d’opérations 
 fictives, strictement réversibles (effectuées à l’aide de pistons semi-per- 
méables et d’une chambre de réaction auxiliaire), susceptibles de trans- 
former en ammoniac le mélange stoechiométrique d'hydrogène et d’azote. 
On supprime cette contradiction en appelant l’état du mélange 3H, + N,, 
non pas un état de « faux équilibre » (comme le faisait Duhem) mais un 
- état d’équilibre relatif, parfaitement compatible avec la réversibilité des 
processus dont il est le point de départ ou le terme. 

Il faut évidemment se garder d’attacher une importance excessive à ces 
questions de mots, où des points de vue différents peuvent tour à tour 
paraître décisifs. Il est peut-être plus instructif de rechercher ce que devient 
le degré d’irréversibilité dans le cas général d’un système où la tempéra- 
ture et les concentrations varient d’un point à l’autre. 


Ill. — SYSTÈMES NON UNIFORMES 


L'importance de l’étude des systèmes non uniformes résulte du fait qu’il 
est au fond contradictoire de considérer un système à température uni- 
forme comme absorbant de la chaleur. Celle-ci ne peut, en effet, pénétrer 


versible. En l’état d'équilibre absolu le produit VWEN:W: est donc «dou- 


ee aa 


~ dans le système et s'y propager qu’à la faveur d’un cane de tempera] 
ture dont nous aurons à rechercher | influence sur le degré d’irréversibilite 
des transformations. | Re 
La définition générale de l'accroissement d’entropie pendant | intervalle 
de temps dt s'impose et ne donne lieu à aucune difficulté après ce que 


A 


nous avons dit plus haut sur la définition de l’entropie d'un systém 
ouvert. Considérant comme un petit système ouvert chaque élément | 
volume d’un système fermé, nous sommes amenés à définir la variation 


de l’entropie du système total par 


(21) as= | dsdV 


Dans cette intégrale, s désigne l’entropie par unité de volume du système 
ouvert contenu à l'instant t dans l'élément de volume dV. L'intégration est 
étendue au volume total V occupé par le système. Dans toute la suite ce 
volume V sera considéré comme constant. On supposera en outre que la 
surface délimitant le système total est rigide et qu’à l’intérieur du systèmes 
les vitesses de convection sont suffisamment faibles pour que l'énergiel 
cinétique puisse être négligée. Toutefois, si, dans un système à volume 
variable, les variations de volume ne sont pas très grandes, le travail cor- 
respondant des forces extérieures sera petit à côté des autres variations” 
d'énergie et les raisonnements précédents resteront valables a titre de“ 
première approximation. ; 

Soit C; la concentration locale du composant (i), c’est-à-dire le nombre 
de molécules-grammes du corps (i) par unité de volume à l’intérieur de 
dV. Nous pouvons maintenant appliquer l’équation (9) en supprimant au 
second membre de cette équation le terme en dV et en supprimant l’in- 
dice V dans les autres termes, puisque chaque volume élémentaire dV 
intervenant dans l'intégrale (21) doit être considéré comme constant, de 
même que le volume V. L’équation (21) devient ainsi 


(22) ie { (Se) dT + E sedCs | dV 


ou s; est ’entropie moléculaire partielle du composant (i) à l’intérieur de dV 
et où l'indice c, remplaçant Vindice n, signifie que toutes les concentra- 
tions sont fixes. 

Cherchons maintenant la définition la plus naturelle de la chaleur non 


compensée dans le cas d’un système non uniforme. A cet effet, nous pose- 
rons (') : 


(1) En évitant de considérer la chaleur non compensée. relative à un élément de 


volume. II n’y a pas lieu, en effet, de parler de la chaleur non compensée d’un système 
ouvert. 
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où dq’ dV est la chaleur réduite, relative à la transformation subie par 
l'élément de volume dV. En remplaçant ds par sa valeur, on obtient 
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Mais, conformément à l'équation (108), 


"(m),= Gt), 


# 


| 1 et, puisque les éléments de volume sont fixes, 
(sr) dTdvV = dudV — (ZuidCi) dV ; | ae 
or NS | 
D Par conséquent, 
4 (23) dQ’ = | dudV — Î Zi — Ts;)dCidV gis! | Tag'uV = L 


Au second membre de cette équation, la première intégrale représente 
la chaleur totale dQ absorbée par le système. En effet, le volume de 


chaque élément de volume restant fixe, on a {duav = dU, et, le volume 
total restant constant, dU = dQ. Par conséquent 
4) { dud¥= 40. 
7a v 
a 4 . . 84 . r 
a Afin d’obtenir la valeur de la troisième intégrale au second membre de 
~  Péquation (23), commençons par préciser la signification de dq’’, car la 
chaleur réduite n’a été définie jusqu'ici que dans le cas d’un système 
- fermé. En désignant par k la conductivité thermique en un point de la 


surface du système total, nous pouvons, conformément au théorème de 
Clausius généralisé, exprimer comme suit la chaleur réduite dQ” absorbée 
par le système total 


4 5 k grad Tdi 
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— 
en représentant par n un vecteur unitaire, normal à l’élément de surface 


do. Appliquant le théorème d’Ostrogradski sur la divergence, nous pou- 
vons remplacer l’intégrale étendue à la surface ©, par une intégrale étendue 


au volume V, et écrire 
EVill, I 6 


= Le 
| Whe CPR 
(26) pu i diy (sere) dV 


Cette transformation n’implique aucune bypothése physique sur les 
systèmes élémentaires ouverts, et suppose seulement que la conductivité 
thermique puisse se définir en chaque point. 

Comme, d’autre part, il faut définir dq’ de manière à satisfaire à 
l'identité 


(27) dQ” a= f dq” dV 
nous voyons que 
(28) Ha k ae Tdt 


Mais la divergence d’un vecteur au point P multipliée par l’élément de 
volume qui contient ce point est égale au flux du vecteur par la surface 
qui délimite l’élément de volume dV. L’équation (2%) montre par consé- 
quent que notre définition de la chaleur réduite, relative au cas d’un 
système fermé non uniforme, reste applicable au système ouvert qui est 
à instant t contenu dans l’élément de volume dV. 

La valeur de dq” fournie par l’équation (28) peut maintenant être sub- 


stituée en l’équation (23) qui devient ainsi, en tenant compte également 
de (24) : 


(9) dQ— =| Etui — Tsi) dCiaV + dQ — | T div ee we dV 


Nous pouvons exprimer dQ en étendant à la surface du système Vinté- 
grale du flux du vecteur k grad Tdt (où dt peut être traité comme constant). 
Nous obtenons ainsi 


- 


ay == i (k grad Tdt) n do =| div (k grad Tat) dV 


En développant les divergences des produits d’un vecteur par un ou 
deux scalaires, k grad Tet KT” grad T, on obtient, en désignant par w° le 
laplacien : 


(30) dQ = 


aan : 
(31) ( T div ( ee = av = ( kyTarav + 


kg TdtdV + | grad T grad kdtdV 
Vv ~ Vv 
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{ grad T grad kdtdV + | Tk grad T grad T! dtdV 
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= 83. 


Puisque grad T-* = —T™ grad T, l’équation (29) peut s’écrire, en 


remplaçant d’ailleurs u; — Ts; par wi, à at 


3 1 k v ac; | tS 
on d0' = = CL, LEE fie | 
(32) | Q'= ine grad® T E wi Gy) atv 


? > 
Représentons par le vecteur vi la vitesse moyenne avec laquelle chemi- 
nent, à l'instant t, en un point donné du système, les molécules du 
corps (i), et par W, la vitesse de chaque réaction chimique indépendante. 
Nous aurons 


(33) os = 2 New W, — div (v; Cs) 
L’équation (32) prendra ainsi la forme plus explicite 
(34) 3 = perad'T + X wi div (ug Ce) —2 = Ne yiW, ! dv. 


Le degré d’irréversibilité de la transformation du système total, et aussi 
de la transformation du système ouvert, infinitésimal, contenu en chaque 
élément de volume constant, se présente ainsi comme une somme de 
trois termes. Les deux premiers représentent le degré d’irréversibilité des 
phénomènes physiques, le troisième représente le degré d’irréversibilité 
chimique. 

La définition que De Donder a proposé pour l’affinité (thermodyna- 
mique) en l’égalant à (VW) dQ’/dt ne peut donc pas se généraliser. Si 
nous voulions la maintenir dans le cas général d’un systéme dont la 
température et les concentrations ne sont pas uniformes, nous inclurions 
en l’affinité les termes représentant l’irréversibilité des échanges de 
température internes, des processus de diffusion et d’égalisation de 
pression qui s’accomplissent au sein du système en même temps que la 
réaction chimique. Cependant, méme dans le cas d’une température non 
uniforme, le potentiel de réaction local —ZN;w:; peut être appelé sans 


17 
inconvénient l’affinité thermodynamique en l’élément de volume dV. Nous 
n’avons pas démontré, mais nous considérons comme probable, que |’affi- 
nité réelle correspondante, — 2 Ni yi WV, est toujours positive. 
Etudions d’un peu plus prés les deux termes représentant le degré 
d’irréversibilité physique des processus qui se déroulent dans un système 
non uniforme. 


Le terme 
i (x grad*T) aV 


exprime l’irréversibilité purement thermique. Il est nécessairement posi- 


à ne on Sn chaude du nr sque gr: 
‘est petit du second ordre. Nous nu. 4 un € le pi 
ment simple que Virréversibilité peut alors être négligée, comme ¢ e 
toujours été jusqu’ici, lorsqu’on appliquait aux systèmes à températt 
pratiquement uniforme les équations démontrées pour les era 
température strictement uniforme. | 
Considérons une sphère initialement uniforme au point dé vue deta, ; 
température, de la pression et de la composition chimique. Supposons — 
infiniment lents les processus de diffusion et les réactions chimiques et _ 
supposons nul le coefficient de dilatation thermique. Si cette sphère est — 
brusquement plongée dans un bain de température plus élevée, il s’éta- : 
blira à l’intérieur de la sphère un gradient de température, qui, en chaque 
point, ne dépendra que de la distance r de ce point au centre de la sphère. 
Par conséquent, la différence 


À age i T dq’ dV 
: os 

sera négligeable à côté de dQ, puisque (à cause de la constance de k) 

e a : 
=f div (k grad T) dV = ae : rae 
que 
iat 
et que 


dQ ST D) ay. 


On peut ainsi remplacer i} T dq’ dV par dQ et écrire 
dy = | Tas av —aQ 
Cans 


alors méme que la température n’est pas rigoureusement uniforme. 
Considérons enfin le terme 


[ | E wediv (ui Ci) Lay 


En l’absence de réaction chimique, ce terme représente la totalité de 
l'intégrale 


=[ ( Z wi D) dV =— J. 


| 
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Il exprime la vitesse avec laquelle ’énergie libre totale diminue par le 
fait de l’égalisation des concentrations. En absence de réaction chimique, 
pour un systéme fermé, 


f= 
A Vv 
_ de sorte que, si pour certains éléments de volume dCi/dt est positif, cette 
_ dérivée doit être négative pour d’autres éléments de volume. 
Cela étant, considérons l’intégrale J comme une somme d’autant d’inté- 
_ grales J,, J,, ...J;... qu’il y a de composants. Considérons ensuite chaque 
- intégrale J; comme une somme de deux intégrales J‘,et Ji, la première 
_ étendue à tous les éléments de volume où dC,/dt est positif, la deuxième 
_ A tous les éléments de volume où dC,/dt est négatif. La condition pour 
- que lon ait—J =—XJ;—XJ/>0est que les valeurs de y, figurant dans 
_les intégrales J} soient, en moyenne, plus grandes que les valeurs de w, 
figurant dans les intégrales Jj. Si l’on suppose que les w, sont des fonctions 
croissantes des C,, cela revient à dire que la condition —J >0 exprime 
une égalisation générale des concentrations. 

En particulier, si la température du systéme est uniforme, et si ce 
système se réduit à un corps pur de masse moléculaire M, l’énergie libre 
moléculaire s’exprimera par uM et la concentration par p/M, en désignant 

- par p la densité locale, Par conséquent 


nf u Ss dy. 


Il suffit alors que u soit une fonction croissante de p et p une fonction 
croissante de la pression, pour que la condition —J > 0 soit liée à l’égali- 
sation générale des pressions au cours du temps. Il est bien entendu que, 
si l’énergie cinétique par unité de masse est suffisante dans certaines 
régions du système, le transport de matière pourra s'effectuer dans ces 
régions de manière à accentuer temporairement une différence de pres- | 
sion. Le système est alors parcouru par des ondes mécaniques et l’irréver- 
sibilité du phénoméne n’est plus mesurée que par l’amortissement pro- 
gressif de ces ondes (!). Le maximum d’irréversibilité correspond au Cas 
d’une égalisation de pression apériodique (critique). Il est instructif de 
comparer, au point de vue où nous venons de nous placer, Pirréversibilite 
purement thermique avec l’irréversibilité purement mécanique. Il n'existe 
pas d’ondes thermiques. L’égalisation des températures constitue toujours 
un phénomène apériodique, présentant un maximum d’irréversibilite. 
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(1) En définissant le degré d’irréversibilité par dQ/ dt. 


ee avec Patfinite réelle. | 
2° Dans le cas des systèmes non homogènes, on voit s Ra à 
sibilité des transformations chimiques, considérées comme se 
dans le temps, certains termes qui traduisent l’irréversibilité des tran | 
ports internes de chaleur et de matière. L'irréversibili ité totale ne sid fie 


donc plus avec l’affinité. 
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5 x Caleul de la fonction potentielle et des fréquences 
a fondamentales des molécules C,H, D, (x+ ye) 
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| F _ Résumé : Reprenant des travaux antérieurs, on a calculé une fonction potentieile qua- 
- dratique unique comprenant le nombre maximum de constantes compatibles avec la 
.” symétrie du système, soit quinze pour la molécule d’éthylène et ses six substitués, où des 
atomes d'hydrogène sont remplacés par des atomes de deutérium. Dans les trente cas où 
la comparaison du calcul avec des résultats eypérimentaux obtenus pour les substitués 
est possible, la concordance entre les fréquences de vibration calculées et observées se 
fait toujours à moins de 1,4 °/o près ; dans la majorité des cas, l'écart ne dépasse 0,5 °/o. 
L'influence de Panharmonicité des vibrations est importante; on en a tenu compte par 
une correction des masses des atomes d'hydrogène et de deutérium, pour lesquelles on a 
adopté les valeurs 1,088 et 2,126 au lieu de 1,008 et de 2,014. Cette correction, suggérée 
par la théorie de l'oscillateur anharmonique, paraît bien légitimée par les nombreuses 
concordances entre les résultats du calcul et de l’expérience. Sans cette correction de 
masses, ou par des valeurs nettement différentes de celles adoptées, on n’a pas pu aboutir 
à un accord satisfaisant. 


L L'étude théorique de modes fondamentaux de vibration de la molécule 

de l’éthylène a été entreprise en détail par MM. Manneback et Verleysen (’). 
[ls ont représenté l’énergie potentielle U de la molécule par une fonction 
- quadratique, en coordonnées symétriques. Cette fonction peut être décom- 
posée en une somme de quatre termes U; indépendants les uns des autres : 


QU: Ade + BEE + Ciné + QDLEs + Wakins + WiLins (1) 
où i — 1, 2, 3 et 4 correspondent respectivement aux modes bien connus 
S,, Ay, A, et S,. La coordonnée symétrique £ se rapporte aux déplacements 
des hydrogènes suivant la direction de l’axe C—C, et n à la direction per- 
pendiculaire, Z se rapporte aux déplacements des carbones. 

Les paramètres A, B, C, D, E et F étaient choisis de manière à rendre 
compte des fréquences des vibrations de C,H, observées par Bonner (°). 


() Nature, 138, 367 (1936) ; Ann. Soc. sc. Brus. B, 56, 349 (1936) et I, 57, 31 (1937) ; 
Jour. Chem. phys., 5, 989 (1937) 
(2) Jour. Amer. Chem. Soc., 58, 34 (1936). 


Ils ont calculé avec cette fonction les fréquences de vibration des molé- — | 


cules C,D, et C,H,D, provenant de substitution totale ou partielle des 
hydrogènes par des deutériums. 

Avec cette même fonction, nous avons calculé avec M. Lemaître les 
fréquences de la molécule monosubstituée C,H,D ('), en vue de vérifica- 
tions expérimentales. | 

On connaît, depuis lors, une trentaine de raies Raman de ces molécules 
grâce aux nombreuses expériences du professeur M. de Hemptinne et ses 
collaborateurs (7). On connaît également quelques raies d’infra-rouge de la 
molécule complètement lourde, grâce aux observations de M.Sutherland (*). 

La comparaison des fréquences calculées et des fréquences observées fait 
apparaître des écarts systématiques. Les erreurs vont en croissant à mesure 
que les molécules sont plus riches en deutérium. Les erreurs commises 
sur les fréquences du monosubstitué ne dépassent pas 1,2 % ; celles des 
bisubstitués dépassent déjà 3 % ; tandis que celles de la molécule tétra- 
substituée atteignent presque 6 %, ces erreurs élevées se limitant d’ailleurs 
à quelques fréquences seulement. 

La cause d’erreur est évidente. Les auteurs eux-mêmes ont déjà précisé 
que : € Il est probable que les termes supérieurs d’anharmonicilé négligés, 
parce qu’inconnus dans la fonction. potentielle, aient, sur les fréquences 
calculées, une répercussion plus grande que la présence de certains termes 
de couplage dans la fonction potentielle limitée à ses termes quadra- 
tiques ». (4) 

Une autre cause d'erreur moins importante provient de ce que les fré- 
quences de C,H, qu’ils ont adoptées sont mesurées dans des états physiques 
différents. Les raies deS,, A, et A, sont mesurées à l’état de vapeur, tandis 
que celles de S, sont mesurées dans le liquide. 

Nous allons maintenant montrer un procédé très simple, qui, tout en 
conservant absolument la même forme de fonction potentielle qu’ils ont 
adoptée, permet cependant de tenir compte de l’anharmonicité qu’ils ont 
négligée. 

Il existe des relations très simples entre les fréquences Vu et vp des deux 
molécules analogues C,H, et C,D,. En absence de l’anharmonicité, le rap- 
port des produits des fréquences de ces deux molécules ne dépend que de 
masses de l’hydrogène et du deutérium. On trouve les relations suivantes : 


() Nature, 139, 1111 (1987); Ann. Soc. sc. Brux. I, 57, 120 (1937). 

(2) Nature, 140, 393 (1937). 

(*) Résultats aimablement communiqués par l’auteur. 

(*) Il y a aussi à tenir compte du fait important que, dans les travaux cités, les données 
expérimentales n’étaient pas en nombre suffisant pour déterminer sans ambiguïté toutes 
les 15 constantes de la fonction potentielle, spécialement les deux constantes de couplage 


fs Ste C'est sur la fréquence située vers 1500 à 1600 cm2 que cette indétermivation 
a le plus d'effet. 
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ne #7 est la distance C-C ; J est la distance C- H; a est l'angle HCH: May 
met M sont les masses atomiques de l'hydrogène, du deutérium, et du 
carbone. = 

Mais à cause de Panharmonicité, aucune de ces relations n’est te 
exactement par les fréquences observées. 

Nous précisons ici que les données de S, et S, sont des fréquences 
Raman, mesurées à l’état liquide par M. dé Hemptinne. Les trois fré- 
quences de modes A, et A, de C,D, sont mesurées dans l’infra-rouge à l’état 
_ de vapeur par ahead Les quatre fréquences de C,H,, de modes A, 
et A, sont les anciennes mesures de Bonner, également à l’état de vapeur. 
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a 1621 em—1| 2988 em—1] 3107 em—1 3072 em—! 
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1515 em—1| 2192,3 em—1| 2335 + 5 em—1) 2304 em—1. 


CD, 2251 1077,4 ? 780 
981 
Re 1,998 1,868 1,868 1,666 
Re 1,954 1,827 2 1,624 
Rs 1,954 1.824 1,824 1,629 


(2 … 


Dans le Tableau I, R; représente les valeurs numériques du ra 
produits de fréquences, tel que la théorie l'exige, calculées av 
deuxième membre des équations (2). Re représente le rapport effec 
calculé avec les raies expérimentales ; on voit trés clairement leur discor- 
dance. | + i oa à 

C’est justement pour réconcilier ce fait que M. Lemaitre nous a proposé 
de faire une correction sur les masses de l'hydrogène et du deutérium ; le 
choix de ces masses devrait pouvoir rendre compatibles les quatre équa- 
tions (2). 

_ Il a été démontré par M. Manneback (*) que, dans le cas des oscillateurs 
simples comme CH et CD, les constantes spectroscopiques d’anharmonicité 
| sont inversement proportionnelles aux racines carrées des masses réduites 
| des oscillateurs. Dans le cas présent, pour ne pas modifier la masse du 

carbone, nous adoptons pour les masses corrigées : 


my, = Mu (14 — ms, = m, (4 + ee) 


. Le rapport des masses corrigées est imposé par les raies de modeS,, 
j on a c= 0,080 et m!. = 1,088 m! = 2,126 


au lieu des valeurs 1,008 et 2,014. Nous appellerons désormais ces masses 
« masses spectroscopiques ». Dans la suite, nous les écrivons sans prime. 

Dans le Tableau I, Rs représente les valeurs numériques du second 
membre des équations (2), calculées avec les masses spectroscopiques. On 
voit que l’accord entre R, et Rs est établi ; une erreur d’un nombre d’ondes 
sur les mesures permet d’expliquer les petits écarts. 

Avec ces masses spectroscopiques, nous pouvons prédéterminer immé- 
diatement la raie de G,D, qui manque dans la mesure de Sutherland. 
D’aprés notre calcul, cette raie de mode A, doit avoir une fréquence de 
693 + 2 cm. 

Dans la suite, nous allons employer ces masses spectroscopiques pour | 
construire la fonction potentielle de la molécule dé C,H,. L'introduction 
de ces masses se justifie par la cohérence de l’ensemble des résultats. 

Les modes $,, A, et A, ne dépendent chacun que de trois paramètres | 
réellement indépendants. Nous avons maintenant, pour chacun de ces 
modes, deux fréquences expérimentalement déterminées pour C,H, et 
deux fréquences pour C,D,. Grâce aux masses spectroscopiques, nous 
avons justement trois équations, au lieu de quatre équations incompa- 
tibles, entre ces trois paramètres. Leur détermination n’offre aucune diffi- 
culté. Nous notons seulement que nous avons pris pour les raies À, de 


(‘) Ann. Soc. sc. Brux. B, 55, 237 (1935). 


40 et 1 694 em, Dons la limite d'erreur + 5 em— sur f raie 
35 em, nous avons choisi celle qui pone moins ii sur la Bie 
du “bisubstitué cis, 4567 cm1. 
Le mode S, dépend de six on indépendants, mais les trois 
EC iéñtes a C,H, et les trois autres de C,D, donnent six équations 
4 incompatibles. Grane aux masses spectroscopiques, ces équations se ré- 
| duisent à cing. Il subsiste donc une indétermination. Comme S, est le 

mode le plus important parmi les quatre, il vaut la peine de faire inter- 
venir les fréquences des isomères. Ainsi nous levons l’indétermination et 
_déterminons les six paramètres par la condition que, couplant avec le 
mode S,, nous puissions trouver exactement la raie bien connue 4567 cm1, 
fréquence de la molécule bisubstituée, dans la position trans (note 2). 

Dans le tableau suivant, on donne les valeurs des paramètres dont 
_ chaque mode dépend. Nos donnons aussi, dans le Tableau III, les valeurs 
de ces paramètres transformées en forces de liaison ening ues: Toutes les 
‘notations et transformations ont été définies dans le travail de Manneback 
et Verleysen ; nous n’insistons pas ici. - 


; TABLEAU II (2) 

1 

à unité en 10° dynes-cm. 

4 | S, | Ay Ae | S 

: Re. ae M ELA es 

“4 A 41,256 | 9,537 | 15,872 | 17,718 

À B 7,351 (009,537 | 18483 | : 6,877 

’ C 19,538 | 16,833 | 15,872 | 16,314 
D —5,591 | —9,537 | 9,048 | 6,006 
E 8,515 | 7,611 |—15,872 | —15,298 
F —5,881 | —7,611 | —9,048 | —9,063 


©] 


(1) Pour les modes Aj, A2, S les 6 constantes A... F de la fonction potentielle se 
réduisent à 3 distinctes. 


spa 25 ie pote p= 5,399 oie 
PET + vi = 1,276 [ve 1,121 |; = 0,615 | 
gr = 0,184 | go =—0,190 | Lam 4 =—0378 Pr 
f, =—0,387 | | AUCUN E 


“pasa 0601 unité = 405 dynes-cm. 
(a =—y, | 4 | | LE: 


| | ‘ i 
2» On voit que ces valeurs sont nettement supérieures à celles de Manneh 
ch et Verleysen. Ceci se conçoit facilement, puisque les petits termes d’ordre 
supérieur sont contenus implicitement dans notre fonction potentielle. 
Rappelons que les paramètres fe et ye expriment respectivement un cou- 
plage élastique entre les variations de distance C-C, d’une part, et les 
variations de distance G-H, respectivement d’angle HCH, d’autre part. 

Nous avons vérifié que ces coefficients permettent de retrouver les 
fréquences observées de C,H, et de C,D, que nous avons adoptées à 1 em 
près. La différence est due aux petits écarts entre Re et Rs. 

Ayant ces paramètres, il est facile de calculer les fréquences d’autres 
isomères. En voici les résultats, exprimés en nombres d’ondes par cm1. 
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TABLEAU IV 


| Erreur 


Calculées Observees | Calc.-Obs. 
. pour cent 
C,H3D 3094 3104 —10 0,3 
3029 3016 13 0,4 
2986 2965 21 0,7 
2272 2266 6 0,3 
1600 1600 0 0 
1393 1397 —4 0,3 
1269 1285 —16 AG? 
948 = = es 
799 — ur = 
CoH, D, cis 3056 3046 10 0,3 
Ty 2295 2290 5 0,2 À 
1570 1567 3 otre 
1999 a Ex = 
759 763 —4 0,5 
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CHD; 3045 = 4 ra 
2352 a 1. fine 
2974 2272 2 0,1 
2927 2215 12 0,5 
1548 1545 3 0,2 
1256 ht et ax 
1028 ee be ae 

811 ex As di 
720 an Bes re 
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On voit ainsi que, sans rien changer à la structure de la fonction poten- 
tielle de Manneback et Verleysen, l'introduction de masses spectro- 
scopiques permet déjà de réduire notablement les écarts avec l'expérience. 

__ Deux raies, d’intensité extrêmement faible, de fréquences 863 cm1 
et 996, que les expérimentateurs ont identifiées comme dues respective- 
ment à la molécule bisubstituée cis ou TRANS, et à la molécule trisubsti- 


LEA 7 
% 


: ny arg ee Se 
tuée, ne trouvent pas leurs places dans le Tableau IV. Cependant ilr 
pas oublier qu’en plus des vibrations planes, existent des vibre 
gauches. Il est donc naturel, comme c'est le cas pour la raie infra-rouge 
719,5-724,7 cm—! dans les mesures de Sutherland sur CD, que des . 
fréquences gauches se mélent quelquefois a des raies de vibrations planes. 

Toutefois, nous devons signaler que les mesures de raies Raman du ¥ 
monosubstitué, des trois bisubstitués, et du trisubstitué sont faites dans — 
un même état physique : liquide, température —150°C. I] en est de même 
Le pour les raies Raman des modes S, et S, des molécules C,H, et C,D,- 
Toutes les raies ont été obtenues au moyen d’un même spectrographe par 
M. de Hemptinne. Seules, les raies de modes A, et A, de la molécule G,D, 
sont mesurées à l’élat de vapeur, dans Pinfra-rouge, par Sutherland ; 
celles de C,H,, de mêmes modes A, et A,, ont été mesurées sur la vapeur 
également par Bonner. Manquant de données homogènes, nous avons dû 
mêler les données provenant de conditions physiques différentes. 11 n’est 
donc pas étonnant que les fréquences calculées de mode S de la molécule 
bisubstituée, en position trans, soient beaucoup plus voisines des 
fréquences expérimentales que les autres, puisque le mode S du bisubstitué 
est obtenu par couplage de deux modes, et S,, observés dans les mêmes 
conditions physiques que le mode S lui-même. 

Avec les méthodes de calcul que nous indiquons dans les notes, il serait 
possible de déterminer les fréquences infra-rouges que donneraient les 
molécules C,H, et C,D, à l'état liquide, en s'appuyant sur les données 
Raman du monosubstitué, des bisubstitués, et du trisubstitué. Nous 
pourrions ainsi rendre homogènes toutes les données, et par suite, déter- 
miner la fonction potentielle avec une approximation encore plus satis- 
faisante. 

Nous terminons en remerciant vivement MM. les professeurs Lemaître 
et Manneback, des directives qu’ils nous ont données, et des bons conseils 
qu’ils nous ont prodigués. 


NOTE | 


Les mouvements de vibration des molécules monosubstituées et bisub- “ 
stituées peuvent être considérés comme les mouvements de C,H, soumis _ 
à des perturbations. Pour mieux le montrer, faisons le changement de « 


variables : 
A; BEM Ci 
A; = -——— 2 —— (ta a . 
9M k, B Wal ; G Amu Kv’, 
D; Ei F; 
Dj eus E, =} OK = Et 
V8M Mu V8M Ma + Amn i . 


dans lequel k =1306 cm— est un coefficient de réduction qui nous 


permet d’introduire directement les fréquences mesurées en nombres 
d’ondes dans les équations de mouvement. 

Les équations de mouvement des molécules C,H,, C,H,D et C,H,D,, 
Apres ce changement de variables, peuvent être mises sous la même Re, 


| (A, — v°)2,+ DE + E,n,=0 
(1) 7=1,2,3et4 D, +(B,— v')£, + Fin, SAV EE. 
E,7,+ Beep Gv yee Nv En, 


dans lesquelles À est une constante caractéristique de la molécule consi- 
dérée. 
Quand À =0, les deuxièmes membres n’existent plus; on reconnaît 
alors les quatre systèmes séparés des équations du mouvement de U,H,. 
On peut regarder les deuxiémes membres comme étant des termes de 
perturbation, sans vouloir dire par là que ces termes sont petits. Il suffit 


; = De el 
de donner à À la valeur oe pour pouvoir retrouver les équations 
H 


du mouvement de monosubstitué. 
On reconnait aussi facilement les équations de mouvement des bisubsti- 


2 (Mn — 


tués en y mettant À — 7 Ma), Les séparations en deux systèmes de 


H 

six équations correspondant 4 deux modes de mouvement sont déja 

indiquées dans le travail de Manneback et Verleysen. Insistons seulement 

que les signes 2 ne s’étendent qu’à deux indices seulement dans ce cas-la. 
Dans les deux cas, les fréquences de vibration sont les racines non nulles 


de l’équation 
Gs pay (Der ht SE 


où les À, sont les déterminants formés par les coefficients du premier 
membre des équations (1), et les X,, Y, et Z, en sont certains mineurs. 
Les racines s’obtiennent facilement en suivant la méthode qui a déjà été 
indiquée. 

On voit donc l’avantage de ce changement de variables. Il nous permet 
pratiquement de traiter les trois molécules isomères au moyen des mêmes 
coefficients A,, B,, C,, etc., et mêmes polynomes X,, Y,, Z; et A,, et de ce 
fait, on peut facilement déterminer, en cas de nécessité, les paramètres 
d’un système partiel, correspondant à une valeur particulière de 7, en 
connaissant des valeurs des fréquences de certains substitués. 


Note II 


Les relations entre les fréquences et paramètres de mode S, sont (’) 


A+B+C=iv=a, 
A+ u(B+ C)=riv=a, 
BC —F°4 AC + AB— D? — E?=v2v2=6, 


u2 (BC — F2) + u (AC + AB — D?— Ex; v2? = b, 


ABC + 2 DEF — AF?— BE? — CD? = Ti v2 =c, 


4 My 
ou 4=—- 
m, 


Les équations (1) et (2) donnent 


a, Ha, 
A= du 
H ay 
Les équations (3) et (4) donnent 
Bo — pe Ms Po 
u(1 — p) 


comme 
(B— C+ 4F*=(B+ C)—4(BC— F))=f4f 


(4) 
(2) 
(3) 
(4) 
(5) 


(6) 


(7) 


(8) 


est une quantité connue par (6) et (7), nous pouvons écrire, en introdui- 


sant un paramètre angulaire @ tel que 
B—C=2fcosé, 
F— fsiné®, 
alors 
B=g+ fcos®, 
C=g—feos®. 
Les équations (3) et (5) permettent d’écrire 
D’+ E*= A(B + OC) + (BC —F°) —b =’, 
CD’ + BE’ — 2 DEF = A(BC — F’) —c¢,=j, 
h et j sont des quantités connues. 
En remplaçant les valeurs B, C et F dans (12), on a 


g(D* + E°) — f(D’ — E*) cos 8 — 9 DE f sin 6 =j, 


() Nous supprimons partout l’indice 1 dont il faudrait affecter les A, B, etc. 


(9) 


(10) 


(11) 
(12) 


= ADEE . 


| Gi AD*ES) cost = =. 4D'E? sinto —4n DE sit sin 8 a nt, 


D°E? — nDE sin 0 + wa Mere Lg ci 


# 


cos 8 


7 ePasons-" FE = COS p, ce qui est une quantité connue. 
Ona me 


“4 ht 
. | DE — sin (6 +), 


Le 


en résolvant (41) et (13), ona: 


D—+hcs 5, 


1 AT 


4 


ÆE=then +9. 

Ainsi, sauf la constante A, les cinq autres paramètres peuvent être 
exprimés en fonction d’un seul paramètre. Ce résultat met bien en 
. évidence l’indétermination dont il a été question dans le texte. 

L’équation générale de la Note | permet de déterminer ce paramètre, 
connaissant une fréquence d’un des substitués. 


4 Louvain, Institut de Physique de l'Université. 
| Février 1938. 
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SUR LA TAILLE DES LAMES DE QUARTZ PRPENDICURAL EMS 
A L'AXE OPTIQUE a 


_par M. A. Brot 


Sommaire. — A l'aide d’un conoscope muni d’un plateau tournant et d’un 
réticule dans le plan focal arrière de l'objectif, il est possible de déterminer les 
erreurs de perpendicularité des lames de quartz à faces planes et parallèles. Il 
suffit pour cela de les incliner convenablement. Pour obtenir une précision 
de l’ordre de 1’ il convient de placer le quartz non directement sur le plateau 
mais sur une charnière intermédiaire cs le mouvement est commandé par 
une vis graduée. 


Te) eke probleme. - — La question a été étudiée par M. R. Servant 
dans le numéro de ‘juin-juillet 1937 de la Revue d'Optique. Con- 
sidérant une lame de quartz à faces planes et parallèles taillée presque 
perpendiculairement à l’axe optique, il établit d’abord mathématique- 


ment la propriété suivante suffisamment approchée tant que l'erreur © 


de taille est faible : 
Plaçons la lame dans un conoscope. Les anneaux que l’on observe 


sont circulaires et centrés sur un point correspondant 2 à l’axe optique 


du quartz. 

cue proposition est d’ailleurs quasi-évidente. 

M. Servant décrit ensuite un dispositif dont l’emploi permet de 
déceler des erreurs de taille de l’ordre de 1’. 

Pendant plusieurs années, utilisant la proposition rappelée ci- 
dessus, nots nous sommes servi d’une méthode différente de la 
sienne et qui est assez simple pour qu'il soit utile de la faire connaître. 

20) Remarques préliminaires. — A. Quand on doit utiliser une lame 


tnt vu 


de quartz perpendiculaire à l’axe, il importe évidemment de déter- | 


miner avec quelle précision il convient de la tailler. Cette précision 
dépend naturellement du problème que la lame sert à résoudre. Dans 


< 


le cas des lames pour saccharimètres par exemple une précision de — 
5’ semble suffisante, On sait, en effet, que pour les pièces de quartz « 
utilisées dans ces instruments, ce n’est qu’à partir d’une inclinaison 


de plus de 10’ sur l’axe optique que l’état de polarisation des rayons * 


lumineux commence à devenir pratiquement différent de celui qui 


6 


existe sur cet axe. Cette précision de 5’, qui suffit dans la plupart des — 


a 


Hi je 


moyens très réduits par une méthode que nous allons décrire. Une 


précision de l’ordre de 1’, que l’on peut obtenir également, exige 


par contre des moyens un peu plus élaborés. 

B. Nous utilisons un conoscope. Dans la forme classique de cet 
instrument le cristal à l'étude se trouve placé entre deux pièces 
d'optique pour lesquelles on peut prendre des condensateurs de 
microscope. Il va sans dire que, selon la précision requise dans les 
observations, selon aussi l'épaisseur du cristal à l'étude, on pourra 
être amené à utiliser d’autres systèmes optiques qu’il faudra adapter 
aux circonstances. Remarquons toutefois qu'il est illusoire d’aug- 
menter, dans un but de précision, la distance focale du système 
objectif au-delà de certaines limites : les franges que l’on observe au 


foyer de ce système sont en effet d'autant plus larges, done moins 
_ nettes, que sa distance focale est plus grande. 


C. L'examen att conoscope est possible sur des lames épaisses de 
plus de 5 cm. Il faut naturellement pour cela que la « pureté » optique 
du cristal soit suffisante. Il s’agit ici de l'absence non seulement de 
mâcles de cristaux droits et gauches donnant des spirales, maïs encore 
de ces défauts locaux qui ne font que déformer les anneaux, leur 
donnant par exemple une forme ovale. Peut-être est-il utile de signaler 
ici que parmi les nombreux blocs de quartz que nous avons étudiés, 
les plus purs, au point de vue optique, étaient ceux qui avaient été 
roulés, provenant par conséquent de cristaux assez volumineux, et 
non pas les petits cristaux assez bien formés. 

30) Le dispositif simple. — Considérons un conoscope réglé à 
l'extinction. Dans le plan focal de l'objectif, 1a où se produisent les 
franges, plaçons un réticule gravé de traits parallèles. Il est commode 
que l’on puisse le déplacer dans son plan perpendiculairement aux 
traits : il est possible ainsi d'amener à tout moment en coïncidence 


un trait du micromètre et le bord d’une frange observée. Le quartz 


est placé, comme on sait, entre le condensateur et l'objectif. | 
Supposons d’abord qu'il soit taillé bien perpendiculairement a son 
axe optique. Pratiquement les franges sont des circonférences centrées, 
même si l’axe du quartz ne coincide pas avec l'axe optique du système 
de lentilles, c’est-à-dire si la lame est inclinée sur l'axe optique de ce 
systéme. Si le quartz est porté par un plateau, il nest pas nêces- 
saire (1) que le plan de celui-ci soit perpendiculaire à l’axe optique des 


lentilles. Mais si ce plateau est rotatif, il convient que l’axe de rota- 
tion soit perpendiculaire à la face qui reçoit le quartz avec une préci- 


(1) Pour la seule observation de la figure d’interférence, 


cas de la pratique, est facile à atteindre et rapidement, avec des . 


~ par contre, les pees nent 

ee si l’on soulève convenablement le 

. rotation, réfracté dans la lame, coïncide dans celle-c 
Tr quartz. Si € est l’angle d’inclinaison de l’a Pour fo qu 

perpendiculaire aux faces de la lame, l'angle toute ik 

le quartz vaut n€, n étant l'indice du En suivant l'axe, 


Loire aes lamelles dieginene connue, ce qui eS de ce 
l'angle dont on l’incline. Par ailleurs, il est clair que l’axe du q 
. se trouve dans le plan perpendiculaire a l'e arête du BRIE bis 
= que l’on crée dessous en le soulevant. 10048 
40) Le dispositif plus précis. — L'usage des cales pets devenir 
pénible, voire délicat, si l’on recherche une précision de l’ordre de 1. 

Si l’on veut éviter cet inconvénient, on ne placera pas le quartz 
directement sur le plateau tournant, mais sur un intermédiaire 
consistant en une sorte de charnière dont l’un des côtés peut se soil 

lever sous l’action d’une vis graduée, permettant d'évaluer le cen-. 
tième de mm. On peut réaliser l’ensemble du dispositif de différentes 
manières. On veillera seulement à ce que la lame de quartz puisse 
venir se placer très près du condensateur du conoscope. 
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